Showing 1–50 of 173 results for author: Nesvorny, D

The Dynamical Origins of the Dark Comets and a Proposed Evolutionary Track

Authors: Aster G. Taylor, Jordan K. Steckloff, Darryl Z. Seligman, Davide Farnocchia, Luke Dones, David Vokrouhlicky, David Nesvorny, Marco Micheli

Abstract: So-called 'dark comets' are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubbl… ▽ More So-called 'dark comets' are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the $ν_6$ resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 17 pages, 7 figures, 12 supplementary figures & pages. Accepted for publication in Icarus

Self-Gravitational Dynamics Within the Inner Oort Cloud

Authors: Konstantin Batygin, David Nesvorný

Abstract: The formation of the Inner Oort Cloud (IOC) - a vast halo of icy bodies residing far beyond Neptune's orbit - is an expected outcome of the solar system's primordial evolution within a stellar cluster. Recent models have shown that the process of early planetesimal capture within the trans-Neptunian region may have been sufficiently high for the cumulative mass of the Cloud to approach several Ear… ▽ More The formation of the Inner Oort Cloud (IOC) - a vast halo of icy bodies residing far beyond Neptune's orbit - is an expected outcome of the solar system's primordial evolution within a stellar cluster. Recent models have shown that the process of early planetesimal capture within the trans-Neptunian region may have been sufficiently high for the cumulative mass of the Cloud to approach several Earth masses. In light of this, here we examine the dynamical evolution of the IOC, driven by its own self-gravity. We show that the collective gravitational potential of the IOC is adequately approximated by the Miyamoto-Nagai model, and use a semi-analytic framework to demonstrate that the resulting secular oscillations are akin to the von Zeipel-Lidov-Kozai resonance. We verify our results with direct $N$-body calculations, and examine the effects of IOC self-gravity on the long-term behavior of the solar system's minor bodies using a detailed simulation. Cumulatively, we find that while the modulation of perihelion distances and inclinations can occur within an observationally relevant range, the associated timescales vastly surpass the age of the sun, indicating that the influence of IOC self-gravity on the architecture of the solar system is negligible. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 14 pages, 6 figures, accepted for publication in Celestial Mechanics and Dynamical Astronomy

EL meteorites do date the giant planet instability

Authors: Chrysa Avdellidou, Marco Delbo, David Nesvorny, Kevin J. Walsh, Alessandro Morbidelli

Abstract: In our recent work, we combined dynamical simulations, meteoritic data and thermal models as well as asteroid observations to argue that the current parent body of the EL meteorites was implanted into the asteroid belt not earlier than 60 Myr after the beginning of the Solar System and that the most likely capture mechanism was the giant planet orbital instability. In the study "The link between A… ▽ More In our recent work, we combined dynamical simulations, meteoritic data and thermal models as well as asteroid observations to argue that the current parent body of the EL meteorites was implanted into the asteroid belt not earlier than 60 Myr after the beginning of the Solar System and that the most likely capture mechanism was the giant planet orbital instability. In the study "The link between Athor and EL meteorites does not constrain the timing of the giant planet instability" that appeared in arXiv, Izidoro and collaborators argue that the implantation of Athor into the asteroid belt does not necessarily require that the giant planet orbital instability occurred at the implantation time. Here we provide further arguments that, in the end, the giant planet instability is still the most likely dynamical process to implant asteroid Athor into the asteroid main belt between 60 and 100 Myr after the beginning of the Solar System. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 4 pages, 1 figure; comment on arXiv:2404.10828

NEOMOD 3: The Debiased Size Distribution of Near Earth Objects

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Carson Fuls, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Marco Delbo

Abstract: Our previous model (NEOMOD2) for the orbital and absolute magnitude distribution of Near Earth Objects (NEOs) was calibrated on the Catalina Sky Survey observations between 2013 and 2022. Here we extend NEOMOD2 to include visible albedo information from the Wide-Field Infrared Survey Explorer. The debiased albedo distribution of NEOs can be approximated by the sum of two Rayleigh distributions wit… ▽ More Our previous model (NEOMOD2) for the orbital and absolute magnitude distribution of Near Earth Objects (NEOs) was calibrated on the Catalina Sky Survey observations between 2013 and 2022. Here we extend NEOMOD2 to include visible albedo information from the Wide-Field Infrared Survey Explorer. The debiased albedo distribution of NEOs can be approximated by the sum of two Rayleigh distributions with the scale parameters p_V,dark=0.03 and p_V,bright=0.17. We find evidence for smaller NEOs having (on average) higher albedos than larger NEOs; this is likely a consequence of the size-dependent sampling of different main belt sources. These inferences and the absolute magnitude distribution from NEOMOD2 are used to construct the debiased size distribution of NEOs. We estimate 830+/-60 NEOs with diameters D>1 km and 20,000+/-2,000 NEOs with D>140 m. The new model, NEOMOD3, is available via the NEOMOD Simulator -- an easy-to-operate code that can be used to generate user-defined samples (orbits, sizes and albedos) from the model. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Icarus, in press

Terrestrial planet formation from a ring: long-term simulations accounting for the giant planet instability

Authors: J. M. Y. Woo, D. Nesvorny, J. Scora, A. Morbidelli

Abstract: The process leading to the formation of the terrestrial planet remains elusive. In a previous publication, we have shown that, if the first generation of planetesimals forms in a ring at about 1 AU and the gas disk's density peaks at the ring location, planetary embryos of a few martian masses can grow and remain in the ring. In this work, we extend our simulations beyond the gas-disk stage, cover… ▽ More The process leading to the formation of the terrestrial planet remains elusive. In a previous publication, we have shown that, if the first generation of planetesimals forms in a ring at about 1 AU and the gas disk's density peaks at the ring location, planetary embryos of a few martian masses can grow and remain in the ring. In this work, we extend our simulations beyond the gas-disk stage, covering 200 Myr and accounting for the phase of giant planet instability, assumed to happen at different times. About half of the simulations form a pair of Venus and Earth analogues and, independently, about 10% form a Mars analogue. We find that the timing of the giant planet instability affects statistically the terrestrial system's excitation state and the timing of the last giant impacts. Hence a late instability (about 60 to 100 Myr after the Solar system's birth) is more consistent with a late Moon-formation time, as suggested by radioactive chronometers. However, the late veneer mass (LVM: mass accreted after the last giant impact) of Earth-sized planets suffering a giant impact after 80 My is usually an order of magnitude lower than the value inferred from geochemistry. In addition, the final angular momentum deficit (AMD) of the terrestrial planets tends to be too high. We tested the effect on the final AMD of the generation of debris during collisions and found that it is too small to change these conclusions. We argue that the best-case scenario is that the Moon-forming event occurred between 50 and 80 My, possibly just following the giant planet instability. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Icarus

Generation of Low-Inclination, Neptune-Crossing TNOs by Planet Nine

Authors: Konstantin Batygin, Alessandro Morbidelli, Michael E. Brown, David Nesvorny

Abstract: The solar system's distant reaches exhibit a wealth of anomalous dynamical structure, hinting at the presence of a yet-undetected, massive trans-Neptunian body - Planet 9. Previous analyses have shown how orbital evolution induced by this object can explain the origins of a broad assortment of exotic orbits, ranging from those characterized by high perihelia to those with extreme inclinations. In… ▽ More The solar system's distant reaches exhibit a wealth of anomalous dynamical structure, hinting at the presence of a yet-undetected, massive trans-Neptunian body - Planet 9. Previous analyses have shown how orbital evolution induced by this object can explain the origins of a broad assortment of exotic orbits, ranging from those characterized by high perihelia to those with extreme inclinations. In this work, we shift the focus toward a more conventional class of TNOs, and consider the observed census of long-period, nearly planar, Neptune-crossing objects as a hitherto-unexplored probe of the Planet 9 hypothesis. To this end, we carry out comprehensive $N-$body simulations that self-consistently model gravitational perturbations from all giant planets, the Galactic tide, as well as passing stars, stemming from initial conditions that account for the primordial giant planet migration and sun's early evolution within a star cluster. Accounting for observational biases, our results reveal that the orbital architecture of this group of objects aligns closely with the predictions of the P9-inclusive model. In stark contrast, the P9-free scenario is statistically rejected at a $\sim5\,σ$ confidence-level. Accordingly, this work introduces a new line of evidence supporting the existence of Planet 9 and further delineates a series of observational predictions poised for near-term resolution. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 11 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Accretion and Uneven Depletion of the Main Asteroid Belt

Authors: Rogerio Deienno, David Nesvorny, Matthew S. Clement, William F. Bottke, Andre Izidoro, Kevin J. Walsh

Abstract: The main asteroid belt (MAB) is known to be primarily composed of objects from two distinct taxonomic classes, generically defined here as S- and C-complex. The former probably originated from the inner solar system (interior to Jupiter's orbit), while the latter probably from the outer solar system. Following this definition, (4) Vesta, a V-type residing in the inner MAB (a < 2.5 au), is the sole… ▽ More The main asteroid belt (MAB) is known to be primarily composed of objects from two distinct taxonomic classes, generically defined here as S- and C-complex. The former probably originated from the inner solar system (interior to Jupiter's orbit), while the latter probably from the outer solar system. Following this definition, (4) Vesta, a V-type residing in the inner MAB (a < 2.5 au), is the sole D > 500 km object akin to S-complex that potentially formed in-situ. This provides a useful constraint on the number of D > 500 km bodies that could have formed, or grown, within the primordial MAB. In this work we numerically simulate the accretion of objects in the MAB region during the time when gas in the protoplanetary disk still existed, while assuming different MAB primordial masses. We then accounted for the depletion of that population happening after gas disk dispersal. In our analysis, we subdivided the MAB into five sub-regions and showed that the depletion factor varies throughout the MAB. This results in uneven radial- and size-dependent depletion of the MAB. We show that the MAB primordial mass has to be $\lesssim$ 2.14$\times$10$^{-3}$ Earth masses. Larger primordial masses would lead to the accretion of tens-to-thousands of S-complex objects with D > 500 km in the MAB. Such large objects would survive depletion even in the outer sub-regions (a > 2.5 au), thus being inconsistent with observations. Our results also indicate that S-complex objects with D > 200-300 km, including (4) Vesta, are likely to be terrestrial planetesimals implanted into the MAB rather than formed in-situ. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: Accepted for publication in PSJ, 23 pages, 7 figures, 4 tables

Testing MOND on small bodies in the remote solar system

Authors: David Vokrouhlický, David Nesvorný, Scott Tremaine

Abstract: 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-syste… ▽ More 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 sub-critical acceleration regime. Using the Bekenstein-Milgrom 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 $10^2$-$10^5$ 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 > 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. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 45 pages, 23 figures

Young asteroid families as the primary source of meteorites

Authors: M. Brož, P. Vernazza, M. Marsset, F. E. DeMeo, R. P. Binzel, D. Vokrouhlický, D. Nesvorný

Abstract: Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent brea… ▽ More Understanding the origin of bright shooting stars and their meteorite samples is among the most ancient astronomy-related questions that at larger scales has human consequences [1-3]. As of today, only ${\sim}\,6\%$ of meteorite falls have been firmly linked to their sources (Moon, Mars, and asteroid (4) Vesta [4-6]). Here, we show that ${\sim}\,70\%$ of meteorites originate from three recent breakups of $D > 30\,{\rm km}$ asteroids that occurred 5.8, 7.5 and less than ${\sim}\,40$ million years ago. These breakups, including the well-known Karin family [7], took place in the prominent yet old Koronis and Massalia families and are at the origin of the dominance of H and L ordinary chondrites among meteorite falls. These young families distinguish themselves amidst all main belt asteroids by having a uniquely high abundance of small fragments. Their size-frequency distribution remains steep for a few tens of millions of years, exceeding temporarily the production of metre-sized fragments by the largest old asteroid families (e.g., Flora, Vesta). Supporting evidence includes the existence of associated dust bands [8-10], the cosmic-ray exposure ages of H-chondrite meteorites [11,12], or the distribution of pre-atmospheric orbits of meteorites [13-15]. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 61 pages, 20 figures

Orbital and absolute magnitude distribution of Jupiter Trojans

Authors: David Vokrouhlický, David Nesvorný, Miroslav Brož, William F. Bottke, Rogerio Deienno, Carson D. Fuls, Frank C. Shelly

Abstract: Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. He… ▽ More Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. Here we use decade long observations from the Catalina Sky Survey (station G96) to determine the bias-corrected orbital and magnitude distributions of JTs. We distinguish the background JT population, filling smoothly the long-term stable orbital zone about L4 and L5 points, and collisional families. We find that the cumulative magnitude distribution of JTs (the background population in our case) has a steep slope for $H\leq 9$, followed with a moderately shallow slope till $H\simeq 14.5$, beyond which the distribution becomes even shallower. At $H=15$ we find a local power-law exponent $0.38\pm 0.01$. We confirm the asymmetry between the magnitude limited background populations in L4 and L5 clouds characterized by a ratio $1.45\pm 0.05$ for $H<15$. Our analysis suggests an asymmetry in the inclination distribution of JTs, with the L4 population being tighter and the L5 population being broader. We also provide a new catalog of the synthetic proper elements for JTs with an updated identification of statistically robust families (9 at L4, and 4 at L5). The previously known Ennomos family is found to consist of two, overlapping Deiphobus and Ennomos families. △ Less

Submitted 27 January, 2024; originally announced January 2024.

Comments: 46 pages, 36 figures, 4 tables, accepted for publication in The Astronomical Journal

NEOMOD 2: An Updated Model of Near-Earth Objects from a Decade of Catalina Sky Survey Observations

Authors: David Nesvorny, David Vokrouhlicky, Frank Shelly, Rogerio Deienno, William F. Bottke, Eric Christensen, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Davide Farnocchia, Mikael Granvik

Abstract: Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This re… ▽ More Catalina Sky Survey (CSS) is a major survey of Near-Earth Objects (NEOs). In a recent work, we used CSS observations from 2005-2012 to develop a new population model of NEOs (NEOMOD). CSS's G96 telescope was upgraded in 2016 and detected over 10,000 unique NEOs since then. Here we characterize the NEO detection efficiency of G96 and use G96's NEO detections from 2013-2022 to update NEOMOD. This resolves previous model inconsistencies related to the population of large NEOs. We estimate there are 936+/-29 NEOs with absolute magnitude H<17.75 (diameter D>1 km for the reference albedo p_V=0.14). The slope of the NEO size distribution for H=25-28 is found to be relatively shallow (cumulative index 2.6) and the number of H<28 NEOs (D>9 m) is determined to be (1.20+/-0.04)x10^7. Small NEOs have a different orbital distribution and higher impact probabilities than large NEOs. We estimate 0.034+/-0.002 impacts of H<28 NEOs on the Earth per year, which is near the low end of the impact flux range inferred from atmospheric bolide observations. Relative to a model where all NEOs are delivered directly from the main belt, the population of small NEOs detected by G96 shows an excess of low-eccentricity orbits with a=1--1.6 au that appears to increase with H. We suggest that the population of very small NEOs is boosted by tidal disruption of large NEOs during close encounters to the terrestrial planets. When the effect of tidal disruption is (approximately) accounted for in the model, we estimate 0.06+/-0.01 impacts of H<28 NEOs on the Earth per year, which is more in line with the bolide data. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: Icarus, in press. arXiv admin note: text overlap with arXiv:2306.09521

Origin and Evolution of Jupiter's Trojan Asteroids

Authors: William F. Bottke, Raphael Marschall, David Nesvorný, David Vokrouhlický

Abstract: The origin of the Jupiter Trojan asteroids has long been a mystery. Dynamically, the population, which is considerably smaller than the main asteroid belt, librates around Jupiter's stable L4 and L5 Lagrange points, 60 deg ahead and behind Jupiter. It is thought that these bodies were captured into these orbits early in solar system history, but any capture mechanism must also explain why the Troj… ▽ More The origin of the Jupiter Trojan asteroids has long been a mystery. Dynamically, the population, which is considerably smaller than the main asteroid belt, librates around Jupiter's stable L4 and L5 Lagrange points, 60 deg ahead and behind Jupiter. It is thought that these bodies were captured into these orbits early in solar system history, but any capture mechanism must also explain why the Trojans have an excited inclination distribution, with some objects reaching inclinations of 35 deg. The Trojans themselves, individually and in aggregate, also have spectral and physical properties that appear consistent with many small bodies found in the outer solar system (e.g., irregular satellites, Kuiper belt objects). In this review, we assemble what is known about the Trojans and discuss various models for their origin and collisional evolution. It can be argued that the Trojans are unlikely to be captured planetesimals from the giant planet zone, but instead were once denizens of the primordial Kuiper belt, trapped by the events taking place during a giant planet instability. The Lucy mission to the Trojans is therefore well positioned to not only answer questions about these objects, but also about their place in planet formation and solar system evolution studies. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 30 pages, 4 figures

Journal ref: Space Sci Rev 219, 83 (2023)

A Crater Chronology for the Jupiter's Asteroids

Authors: S. Marchi, D. Nesvorný, D. Vokrouhlický, W. F. Bottke, H. Levison

Abstract: We present a new crater chronology for Jupiter's Trojan asteroids. This tool can be used to interpret the collisional history of the bodies observed by NASA's Lucy mission. The Lucy mission will visit a total of six Trojan asteroids: Eurybates, Polymele, Orus, Leucus, and the near equal mass binary Patroclus-Menoetius. In addition, Eurybates and Polymele each have a small satellite. Here we presen… ▽ More We present a new crater chronology for Jupiter's Trojan asteroids. This tool can be used to interpret the collisional history of the bodies observed by NASA's Lucy mission. The Lucy mission will visit a total of six Trojan asteroids: Eurybates, Polymele, Orus, Leucus, and the near equal mass binary Patroclus-Menoetius. In addition, Eurybates and Polymele each have a small satellite. Here we present a prediction of Trojan cratering based on current models of how the Solar System and the objects themselves evolved. We give particular emphasis to the time lapsed since their implantation into stable regions near Jupiter's Lagrangian L4 and L4 points. We find that cratering on Trojans is generally dominated by mutual collisions, with the exception of a short period of time (~10 Myr) after implantation, in which cometary impacts may have been significant. For adopted crater scaling laws, we find that the overall spatial density of craters on Trojans is significantly lower than that of Main Belt asteroids on surfaces with similar formation ages. We also discuss specific predictions for similar-sized Eurybates and Orus, and the binary system Patroclus-Menoetius. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: Published on AJ

Journal ref: AJ 166 221 (2023)

Isotopic Trichotomy of Main Belt Asteroids from Implantation of Outer Solar System Planetesimals

Authors: David Nesvorny, Nicolas Dauphas, David Vokrouhlicky, Rogerio Deienno, Timo Hopp

Abstract: Recent analyses of samples from asteroid (162173) Ryugu returned by JAXA's Hayabusa2 mission suggest that Ryugu and CI chondrites formed in the same region of the protoplanetary disk, in a reservoir that was isolated from the source regions of other carbonaceous (C-type) asteroids. Here we conduct $N$-body simulations in which CI planetesimals are assumed to have formed in the Uranus/Neptune zone… ▽ More Recent analyses of samples from asteroid (162173) Ryugu returned by JAXA's Hayabusa2 mission suggest that Ryugu and CI chondrites formed in the same region of the protoplanetary disk, in a reservoir that was isolated from the source regions of other carbonaceous (C-type) asteroids. Here we conduct $N$-body simulations in which CI planetesimals are assumed to have formed in the Uranus/Neptune zone at $\sim15$--25 au from the Sun. We show that CI planetesimals are scattered by giant planets toward the asteroid belt where their orbits can be circularized by aerodynamic gas drag. We find that the dynamical implantation of CI asteroids from $\sim15$--25 au is very efficient with $\sim 5$\% of $\sim 100$-km planetesimals reaching stable orbits in the asteroid belt by the end of the protoplanetary gas disk lifetime. The efficiency is reduced when planetesimal ablation is accounted for. The implanted population subsequently evolved by collisions and was depleted by dynamical instabilities. The model can explain why CIs are isotopically distinct from other C-type asteroids which presumably formed at $\sim5$--10 au. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: EPSL, in press

Dynamical feasibility of (3) Juno as a parent body of the H chondrites

Authors: John W. Noonan, Kathryn Volk, David Nesvorný, William F. Bottke

Abstract: We test the hypothesis that (3) Juno is a parent body of the H chondrites with dynamical modeling of an asteroid-family-forming impact and comparison to current observational data. Using a dynamical model that includes the Yarkovsky force on a simulated Juno family and a simplified cosmic ray exposure age model we examine the expected distribution of Juno family members in both the main belt and n… ▽ More We test the hypothesis that (3) Juno is a parent body of the H chondrites with dynamical modeling of an asteroid-family-forming impact and comparison to current observational data. Using a dynamical model that includes the Yarkovsky force on a simulated Juno family and a simplified cosmic ray exposure age model we examine the expected distribution of Juno family members in both the main belt and near-Earth orbits over 300 Myrs and the cosmic ray exposure distribution for fragments exiting the main belt via the 3:1J, 5:2J, and 8:3J mean motion resonances. We find that the smallest modeled ($D<$10 m) family members of (3) Juno cannot be directly responsible for the observed H chondrite flux and that the breakup of larger family members creates an CRE distribution that resembles the measured H chondrite CRE distribution but is still unable to adequately explain the significant number of H chondrites with CRE ages of 6-8 Myrs. A similar model was performed for the asteroid (6) Hebe, another parent body candidate, and produced a CRE age distribution that is inconsistent with the measured H chondrite CRE ages. We also find from our dynamical models that we can expect $<$7 km-scale Juno family members in near-Earth orbits in the present day, consistent with the recent discovery of the shock-darkened H chondrite-like asteroid (52768) 1998 OR$_{2}$. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 16 pages, 9 figures

Journal ref: Volume 408, 15 January 2024, 115838 Volume 408, 15 January 2024, 115838 Icarus Icarus, Volume 408, 15 January 2024, 115838

Debiased population of very young asteroid families

Authors: David Vokrouhlický, David Nesvorný, Miroslav Brož, William F. Bottke

Abstract: We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also… ▽ More We use observations from the Catalina Sky Survey (CSS) to determine the bias-corrected population of small members in four very young families down to sizes equivalent to several hundred meters. Using the most recent catalog of known asteroids, we identified members from four young families for which the population has grown appreciably over recent times. A large fraction of these bodies have also been detected by CSS. We used synthetic populations of asteroids, with their magnitude distribution controlled by a small number of parameters, as a template for the bias-corrected model of these families. Applying the known detection probability of the CSS observations, we could adjust these model parameters to match the observed (biased) populations in the young families. In the case of three families, Datura, Adelaide, and Rampo, we find evidence that the magnitude distribution transitions from steep to shallow slopes near $300$ to $400$ meters. Conversely, the Hobson family population may be represented by a single power-law model. The Lucascavin family has a limited population; no new members have been discovered over the past two decades. We consider a model of parent body rotational fission with the escaping secondary tidally split into two components (thereby providing three members within this family). In support of this idea, we find that no other asteroid with absolute magnitude $H\leq 18.3$ accompanies the known three members in the Lucascavin family. A similar result is found for the archetypal asteroid pair Rheinland--Kurpfalz. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 32 pages, 27 figures, accepted for publication in Astronomy and Astrophysics

Not So Fast Kepler-1513: A Perturbing Planetary Interloper in the Exomoon Corridor

Authors: Daniel A. Yahalomi, David Kipping, David Nesvorný, Paul A. Dalba, Paul Benni, Ceiligh Cacho-Negrete, Karen Collins, Joel T. Earwicker, John Arban Lewis, Kim K. McLeod, Richard P. Schwarz, Gavin Wang

Abstract: Transit Timing Variations (TTVs) can be induced by a range of physical phenomena, including planet-planet interactions, planet-moon interactions, and stellar activity. Recent work has shown that roughly half of moons would induce fast TTVs with a short period in the range of two-to-four orbits of its host planet around the star. An investigation of the Kepler TTV data in this period range identifi… ▽ More Transit Timing Variations (TTVs) can be induced by a range of physical phenomena, including planet-planet interactions, planet-moon interactions, and stellar activity. Recent work has shown that roughly half of moons would induce fast TTVs with a short period in the range of two-to-four orbits of its host planet around the star. An investigation of the Kepler TTV data in this period range identified one primary target of interest, Kepler-1513 b. Kepler-1513 b is a $8.05^{+0.58}_{-0.40}$ $R_\oplus$ planet orbiting a late G-type dwarf at $0.53^{+0.04}_{-0.03}$ AU. Using Kepler photometry, this initial analysis showed that Kepler-1513 b's TTVs were consistent with a moon. Here, we report photometric observations of two additional transits nearly a decade after the last Kepler transit using both ground-based observations and space-based photometry with TESS. These new transit observations introduce a previously undetected long period TTV, in addition to the original short period TTV signal. Using the complete transit dataset, we investigate whether a non-transiting planet, a moon, or stellar activity could induce the observed TTVs. We find that only a non-transiting perturbing planet can reproduce the observed TTVs. We additionally perform transit origami on the Kepler photometry, which independently applies pressure against a moon hypothesis. Specifically, we find that Kepler-1513 b's TTVs are consistent with an exterior non-transiting $\sim$Saturn mass planet, Kepler-1513 c, on a wide orbit, $\sim$5$\%$ outside a 5:1 period ratio with Kepler-1513 b. This example introduces a previously unidentified cause for planetary interlopers in the exomoon corridor, namely an insufficient baseline of observations. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 20 pages, 13 figures. Accepted to MNRAS. Code available at https://github.com/dyahalomi/Kepler1513

Bayesian Inference on the Isotopic Building Blocks of Mars and Earth

Authors: Nicolas Dauphas, Timo Hopp, David Nesvorny

Abstract: Isotopic anomalies provide a means of probing the materials responsible for the formation of terrestrial planets. By analyzing new iron isotopic anomaly data from Martian meteorites and drawing insights from published data for O, Ca, Ti, Cr, Fe, Ni, Sr, Zr, Mo, Ru, and Si, we scrutinize potential changes in the isotopic composition of the material accreted by Mars and Earth during their formation.… ▽ More Isotopic anomalies provide a means of probing the materials responsible for the formation of terrestrial planets. By analyzing new iron isotopic anomaly data from Martian meteorites and drawing insights from published data for O, Ca, Ti, Cr, Fe, Ni, Sr, Zr, Mo, Ru, and Si, we scrutinize potential changes in the isotopic composition of the material accreted by Mars and Earth during their formation. A Principal Component Analysis of isotopic anomalies in meteorites identifies three main clusters (forming the three parts of the isotopic trichotomy): CI, CC=CM+CO+CV+CR, and NC=EH+EL+H+L+LL. Our results suggest that Earth is primarily an isotopic mixture of ~92% E, 6 % CI, and <2% COCV and O. Mars, on the other hand, appears to be a mixture of ~65% E, 33% O, and <2% CI and COCV. We establish that Earth's CI contribution substantially increased during the latter half of its accretion. Mars began accreting a mix of O and E but predominantly accreted E later. Mars' changing isotopic makeup during accretion can be explained if it underwent gas-driven type I migration from its origin near the O - E boundary to a location well within the E region during the first few million years of solar system history. Earth's later increased CI contribution may be attributed to the stochastic impact of an interloper carbonaceous embryo that moved inside the inner solar system region while nebular gas was still present, and subsequently participated in the stage of chaotic growth. The recent findings of Si isotopic anomalies in enstatite chondrites when compared to terrestrial rocks likely stems from insufficient correction for high-temperature equilibrium isotopic fractionation. With appropriate adjustments for this influence, both the silicate Earth and enstatite chondrites exhibit comparable Si isotopic anomalies, reaffirming a genetic link between them. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: In press, Icarus

A race against the clock: Constraining the timing of cometary bombardment relative to Earth's growth

Authors: Sarah Joiret, Sean N. Raymond, Guillaume Avice, Matthew S. Clement, Rogerio Deienno, David Nesvorný

Abstract: Comets are considered a potential source of inner solar system volatiles, but the timing of this delivery relative to that of Earth's accretion is still poorly understood. Measurements of xenon isotopes in comet 67P/Churyumov-Gerasimenko revealed that comets partly contributed to the Earth's atmosphere. However, there is no conclusive evidence of a significant cometary component in the Earth's man… ▽ More Comets are considered a potential source of inner solar system volatiles, but the timing of this delivery relative to that of Earth's accretion is still poorly understood. Measurements of xenon isotopes in comet 67P/Churyumov-Gerasimenko revealed that comets partly contributed to the Earth's atmosphere. However, there is no conclusive evidence of a significant cometary component in the Earth's mantle. These geochemical constraints would favour a contribution of comets mainly occurring after the last stages of Earth's formation. Here, we evaluate whether dynamical simulations satisfy these constraints in the context of an Early Instability model. We perform dynamical simulations of the solar system, calculate the probability of collision between comets and Earth analogs component embryos through time and estimate the total cometary mass accreted in Earth analogs as a function of time. While our results are in excellent agreement with geochemical constraints, we also demonstrate that the contribution of comets on Earth might have been delayed with respect to the timing of the instability, due to a stochastic component of the bombardment. More importantly, we show that it is possible that enough cometary mass has been brought to Earth after it had finished forming so that the xenon constraint is not necessarily in conflict with an Early Instability scenario. However, it appears very likely that a few comets were delivered to Earth early in its accretion history, thus contributing to the mantle's budget. Finally, we compare the delivery of cometary material on Earth to Venus and Mars. These results emphasize the stochastic nature of the cometary bombardment in the inner solar system. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 26 pages, 12 figures

Journal ref: Icarus, 406, 115754 (2023)

The Properties and Origins of Kuiper Belt Object Arrokoth's Large Mounds

Authors: S. A. Stern, O. L. White, Wm. Grundy, B. A. Keeney, J. D. Hofgartner, D. Nesvorny, W. B. McKinnon, D. C. Richardson, J. C. Marohnic, A. J. Verbiscer, S. D. Benecchi, P. M. Schenk, J. M. Moore

Abstract: We report on a study of the mounds that dominate the appearance of Kuiper Belt Object (KBO) (486958) Arrokoth's larger lobe, named Wenu. We compare the geological context of these mounds, measure and intercompare their shapes, sizes/orientations, reflectance, and colors. We find the mounds are broadly self-similar in many respects and interpret them as the original building blocks of Arrokoth. It… ▽ More We report on a study of the mounds that dominate the appearance of Kuiper Belt Object (KBO) (486958) Arrokoth's larger lobe, named Wenu. We compare the geological context of these mounds, measure and intercompare their shapes, sizes/orientations, reflectance, and colors. We find the mounds are broadly self-similar in many respects and interpret them as the original building blocks of Arrokoth. It remains unclear why these building blocks are so similar in size, and this represents a new constrain and challenge for solar system formation models. We then discuss the interpretation of this interpretation. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 24 pages, 8 figures

Radial Distribution of Distant Trans-Neptunian Objects Points to Sun's Formation in a Stellar Cluster

Authors: David Nesvorny, Pedro Bernardinelli, David Vokrouhlicky, Konstantin Batygin

Abstract: The Scattered Disk Objects (SDOs) are a population of trans-Neptunian bodies with semimajor axes $50< a \lesssim 1000$ au and perihelion distances $q \gtrsim 30$ au. The detached SDOs with orbits beyond the reach of Neptune (roughly $q>35$~au) are of special interest here as an important constraint on the early evolution of the outer Solar System. The semimajor axis profile of detached SDOs at 50-… ▽ More The Scattered Disk Objects (SDOs) are a population of trans-Neptunian bodies with semimajor axes $50< a \lesssim 1000$ au and perihelion distances $q \gtrsim 30$ au. The detached SDOs with orbits beyond the reach of Neptune (roughly $q>35$~au) are of special interest here as an important constraint on the early evolution of the outer Solar System. The semimajor axis profile of detached SDOs at 50--500~au, as characterized from the Dark Energy Survey (DES), is radially extended, but previous dynamical models of Neptune's early migration produce a relatively compact profile. This problem is most likely related to Sun's birth environment in a stellar cluster. We perform new dynamical simulations that account for cluster effects and show that the orbital distribution of SDOs can be explained if a particularly close stellar encounter occurred early on (e.g., M dwarf with the mass $\simeq 0.2$ $M_\odot$ approaching the Sun at $\simeq 200$ au). For such an encounter to happen with a reasonably high probability the Sun must have formed in a stellar cluster with $ηT \gtrsim 10^4$ Myr pc$^{-3}$, where $η$ is the stellar number density and $T$ is the Sun's residence time in the cluster. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Shape models and spin states of Jupiter Trojans: Testing the streaming instability formation scenario

Authors: Josef Hanuš, David Vokrouhlický, David Nesvorný, Josef Ďurech, Robert Stephens, Vladimir Benishek, Julian Oey, Petr Pokorný

Abstract: The leading theory for the origin of Jupiter Trojans (JTs) assumes that JTs were captured to their orbits near the Lagrangian points of Jupiter during the early reconfiguration of the giant planets. The natural source region for the majority of JTs would then be the population of planetesimals born in a massive trans-Neptunian disk. If true, JTs represent the most accessible stable population of s… ▽ More The leading theory for the origin of Jupiter Trojans (JTs) assumes that JTs were captured to their orbits near the Lagrangian points of Jupiter during the early reconfiguration of the giant planets. The natural source region for the majority of JTs would then be the population of planetesimals born in a massive trans-Neptunian disk. If true, JTs represent the most accessible stable population of small Solar System bodies that formed in the outer regions of the Solar System. For this work, we compiled photometric datasets for about 1000 JTs and applied the convex inversion technique in order to assess their shapes and spin states. We obtained full solutions for $79$ JTs, and partial solutions for an additional $31$ JTs. We found that the observed distribution of the pole obliquities of JTs is broadly consistent with expectations from the streaming instability, which is the leading mechanism for the formation of planetesimals in the trans-Neptunian disk. The observed JTs' pole distribution has a slightly smaller prograde vs. retrograde asymmetry (excess of obliquities $>130^\circ$) than what is expected from the existing streaming instability simulations. However, this discrepancy can be plausibly reconciled by the effects of the post-formation collisional activity. Our numerical simulations of the post-capture spin evolution indicate that the JTs' pole distribution is not significantly affected by dynamical processes such as the eccentricity excitation in resonances, close encounters with planets, or the effects of nongravitational forces. However, a few JTs exhibit large latitude variations of the rotation pole and may even temporarily transition between prograde- and retrograde-rotating categories. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A56 (2023)

Impact Rates in the Outer Solar System

Authors: David Nesvorny, Luke Dones, Mario De Pra, Maria Womack, Kevin J. Zahnle

Abstract: Previous studies of cometary impacts in the outer Solar System used the spatial distribution of ecliptic comets (ECs) from dynamical models that assumed ECs began on low-inclination orbits (<5 deg) in the Kuiper belt. In reality, the source population of ECs - the trans-Neptunian scattered disk - has orbital inclinations reaching up to ~30 deg. In Nesvorny et al. (2017), we developed a new dynamic… ▽ More Previous studies of cometary impacts in the outer Solar System used the spatial distribution of ecliptic comets (ECs) from dynamical models that assumed ECs began on low-inclination orbits (<5 deg) in the Kuiper belt. In reality, the source population of ECs - the trans-Neptunian scattered disk - has orbital inclinations reaching up to ~30 deg. In Nesvorny et al. (2017), we developed a new dynamical model of ECs by following comets as they evolved from the scattered disk to the inner Solar System. The model was absolutely calibrated from the population of Centaurs and active ECs. Here we use our EC model to determine the steady-state impact flux of cometary/Centaur impactors on Jupiter, Saturn, Uranus, and their moons. Relative to previous work (Zahnle et al. 2003), we find slightly higher impact probabilities on the outer moons and lower impact probabilities on the inner moons. The impact probabilities are smaller when comet disruption is accounted for. The results provide a modern framework for the interpretation of the cratering record in the outer Solar System. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: PSJ, in press

The Collisional Evolution of the Primordial Kuiper Belt, Its Destabilized Population, and the Trojan Asteroids

Authors: William Bottke, David Vokrouhlicky, Raphael Marshall, David Nesvorny, Alessandro Morbidelli, Rogerio Deienno, Simone Marchi, Luke Dones, Harold Levison

Abstract: The tumultuous early era of outer solar system evolution culminated when Neptune migrated across the primordial Kuiper belt (PKB) and triggered a dynamical instability among the giant planets. This event led to the ejection of approximately 99.9\% of the PKB (here called the destabilized population), heavy bombardment of the giant planet satellites, and the capture of Jupiter's Trojans. While this… ▽ More The tumultuous early era of outer solar system evolution culminated when Neptune migrated across the primordial Kuiper belt (PKB) and triggered a dynamical instability among the giant planets. This event led to the ejection of approximately 99.9\% of the PKB (here called the destabilized population), heavy bombardment of the giant planet satellites, and the capture of Jupiter's Trojans. While this scenario has been widely tested using dynamical models, there have been fewer investigations into how the PKB, its destabilized population, and the Trojans experienced collisional evolution. Here we examined this issue for all three populations with the code Boulder. Our constraints included the size-frequency distributions (SFDs) of the Trojan asteroids and craters on the giant planet satellites. Using this combination, we solved for the unknown disruption law affecting bodies in these populations. The weakest ones, from an impact energy per mass perspective, were 20 m in diameter. Overall, collisional evolution produces a power-law-like shape for multikilometer Trojans and a wavy-shaped SFD in the PKB and destabilized populations. The latter can explain (i) the shapes of the ancient and younger crater SFDs observed on the giant planet satellites, (ii) the shapes of the Jupiter family and long-period comet SFDs, which experienced different degrees of collision evolution, and (iii) the present-day impact frequency of superbolides on Jupiter and smaller projectiles on Saturn's rings. Our model results also indicate that many observed comets, most which are smaller than 10 km in diameter, are likely to be gravitational aggregates formed by large-scale collision events. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 95 pages, 19 figures. Accepted for publication in PSJ

NEOMOD: A New Orbital Distribution Model for Near Earth Objects

Authors: David Nesvorny, Rogerio Deienno, William F. Bottke, Robert Jedicke, Shantanu Naidu, Steven R. Chesley, Paul W. Chodas, Mikael Granvik, David Vokrouhlicky, Miroslav Broz, Alessandro Morbidelli, Eric Christensen, Bryce T. Bolin

Abstract: Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected fro… ▽ More Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $ν_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 \times 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $θ\simeq 0^\circ$ for $a<2.5$~au and $θ\simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: AJ

Early Bombardment of the Moon: Connecting the Lunar Crater Record to the Terrestrial Planet Formation

Authors: David Nesvorny, Fernando V. Roig, David Vokrouhlicky, William F. Bottke, Simone Marchi, Alessandro Morbidelli, Rogerio Deienno

Abstract: The lunar crater record features $\sim 50$ basins. The radiometric dating of Apollo samples indicates that the Imbrium basin formed relatively late -- from the planet formation perspective -- some $\simeq 3.9$ Ga. Here we develop a dynamical model for impactors in the inner solar system to provide context for the interpretation of the lunar crater record. The contribution of cometary impactors is… ▽ More The lunar crater record features $\sim 50$ basins. The radiometric dating of Apollo samples indicates that the Imbrium basin formed relatively late -- from the planet formation perspective -- some $\simeq 3.9$ Ga. Here we develop a dynamical model for impactors in the inner solar system to provide context for the interpretation of the lunar crater record. The contribution of cometary impactors is found to be insignificant. Asteroids produced most large impacts on the terrestrial worlds in the last $\simeq 3$ Gyr. The great majority of early impactors were rocky planetesimals left behind at $\sim 0.5$--1.5 au after the terrestrial planet accretion. The population of terrestrial planetesimals was reduced by disruptive collisions in the first $t \sim 20$ Myr after the gas disk dispersal. We estimate that there were $\sim 4 \times 10^5$ diameter $d>10$ km bodies when the Moon formed (total planetesimal mass $\sim 0.015$ $M_{\rm Earth}$ at $t \sim 50$ Myr). The early bombardment of the Moon was intense. To accommodate $\sim 50$ known basins, the lunar basins that formed before $\simeq 4.35$--4.41 Ga must have been erased. The late formation of Imbrium occurs with a $\sim 15$--35\% probability in our model. About 20 $d>10$-km bodies were expected to hit the Earth between 2.5 and 3.5 Ga, which is comparable to the number of known spherule beds in the late Archean. We discuss implications of our model for the lunar/Martian crater chronologies, Late Veneer, and noble gases in the Earth atmosphere. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: Icarus, in press

TOI-2525 b and c: A pair of massive warm giant planets with a strong transit timing variations revealed by TESS

Authors: Trifon Trifonov, Rafael Brahm, Andres Jordan, Christian Hartogh, Thomas Henning, Melissa J. Hobson, Martin Schlecker, Saburo Howard, Finja Reichardt, Nestor Espinoza, Man Hoi Lee, David Nesvorny, Felipe I. Rojas, Khalid Barkaoui, Diana Kossakowski, Gavin Boyle, Stefan Dreizler, Martin Kuerster, Rene Heller, Tristan Guillot, Amaury H. M. J. Triaud, Lyu Abe, Abdelkrim Agabi, Philippe Bendjoya, Nicolas Crouzet , et al. (22 additional authors not shown)

Abstract: TOI-2525 is a K-type star with an estimated mass of M = 0.849$_{-0.033}^{+0.024}$ M$_\odot$ and radius of R = 0.785$_{-0.007}^{+0.007}$ R$_\odot$ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semi-amplitude of a $\sim$6 hours. We performed TTV… ▽ More TOI-2525 is a K-type star with an estimated mass of M = 0.849$_{-0.033}^{+0.024}$ M$_\odot$ and radius of R = 0.785$_{-0.007}^{+0.007}$ R$_\odot$ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semi-amplitude of a $\sim$6 hours. We performed TTV dynamical, and photo-dynamical light curve analysis of the TESS data, combined with radial velocity (RV) measurements from FEROS and PFS, and we confirmed the planetary nature of these companions. The TOI-2525 system consists of a transiting pair of planets comparable to Neptune and Jupiter with estimated dynamical masses of $m_{\rm b}$ = 0.088$_{-0.004}^{+0.005}$ M$_{\rm Jup.}$, and $m_{\rm c}$ = 0.709$_{-0.033}^{+0.034}$ M$_{\rm Jup.}$, radius of $r_b$ = 0.88$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$ and $r_c$ = 0.98$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$, and with orbital periods of $P_{\rm b}$ = 23.288$_{-0.002}^{+0.001}$ days and $P_{\rm c}$ = 49.260$_{-0.001}^{+0.001}$ days for the inner and the outer planet, respectively. The period ratio is close to the 2:1 period commensurability, but the dynamical simulations of the system suggest that it is outside the mean motion resonance (MMR) dynamical configuration. TOI-2525 b is among the lowest density Neptune-mass planets known to date, with an estimated median density of $ρ_{\rm b}$ = 0.174$_{-0.015}^{+0.016}$ g\,cm$^{-3}$. The TOI-2525 system is very similar to the other K-dwarf systems discovered by TESS, TOI-2202 and TOI-216, which are composed of almost identical K-dwarf primary and two warm giant planets near the 2:1 MMR. △ Less

Submitted 11 February, 2023; originally announced February 2023.

Comments: Accepted for publication in AJ

The gateway from Centaurs to Jupiter-family Comets: thermal and dynamical evolution

Authors: Aurélie Guilbert-Lepoutre, Anastasios Gkotsinas, Sean N. Raymond, David Nesvorny

Abstract: It was recently proposed that there exists a "gateway" in the orbital parameter space through which Centaurs transition to Jupiter-family Comets (JFCs). Further studies have implied that the majority of objects that eventually evolve into JFCs should leave the Centaur population through this gateway. This may be naively interpreted as gateway Centaurs being pristine progenitors of JFCs. This is th… ▽ More It was recently proposed that there exists a "gateway" in the orbital parameter space through which Centaurs transition to Jupiter-family Comets (JFCs). Further studies have implied that the majority of objects that eventually evolve into JFCs should leave the Centaur population through this gateway. This may be naively interpreted as gateway Centaurs being pristine progenitors of JFCs. This is the point we want to address in this work. We show that the opposite is true: gateway Centaurs are, on average, more thermally processed than the rest of the population of Centaurs crossing Jupiter's orbit. Using a dynamically-validated JFC population, we find that only $\sim 20\%$ of Centaurs pass through the gateway prior to becoming JFCs, in accordance with previous studies. We show that more than half of JFC dynamical clones entering the gateway for the first time have already been JFCs -they simply avoided the gateway on their first pass into the inner solar system. By coupling a thermal evolution model to the orbital evolution of JFC dynamical clones, we find a higher than 50\% chance that the layer currently contributing to the observed activity of gateway objects has been physically and chemically altered, due to previously sustained thermal processing. We further illustrate this effect by examining dynamical clones that match the present-day orbits of 29P/Schwassmann-Wachmann 1, P/2019 LD2 (ATLAS), and P/2008 CL94 (Lemmon). △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 16 pages, 13 Figures, Accepted for publication on the ApJ

Comets and Planetesimal Formation

Authors: Jacob B. Simon, Jürgen Blum, Til Birnstiel, David Nesvorný

Abstract: In this chapter, we review the processes involved in the formation of planetesimals and comets. We will start with a description of the physics of dust grain growth and how this is mediated by gas-dust interactions in planet-forming disks. We will then delve into the various models of planetesimal formation, describing how these planetesimals form as well as their resulting structure. In doing so,… ▽ More In this chapter, we review the processes involved in the formation of planetesimals and comets. We will start with a description of the physics of dust grain growth and how this is mediated by gas-dust interactions in planet-forming disks. We will then delve into the various models of planetesimal formation, describing how these planetesimals form as well as their resulting structure. In doing so, we focus on and compare two paradigms for planetesimal formation: the gravitational collapse of particle over-densities (which can be produced by a variety of mechanisms) and the growth of particles into planetesimals via collisional and gravitational coagulation. Finally, we compare the predictions from these models with data collected by the Rosetta and New Horizons missions and that obtained via observations of distant Kuiper Belt Objects. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: Planetesimal Formation Review accepted for publication in Comets III

Formation of Lunar Basins from Impacts of Leftover Planetesimals

Authors: David Nesvorny, Fernando V. Roig, David Vokrouhlicky, William F. Bottke, Simone Marchi, Alessandro Morbidelli, Rogerio Deienno

Abstract: The Moon holds important clues to the early evolution of the Solar System. Some 50 impact basins (crater diameter D>300 km) have been recognized on the lunar surface, implying that the early impact flux was much higher than it is now. The basin-forming impactors were suspected to be asteroids released from an inner extension of the main belt (1.8-2.0 au). Here we show that most impactors were inst… ▽ More The Moon holds important clues to the early evolution of the Solar System. Some 50 impact basins (crater diameter D>300 km) have been recognized on the lunar surface, implying that the early impact flux was much higher than it is now. The basin-forming impactors were suspected to be asteroids released from an inner extension of the main belt (1.8-2.0 au). Here we show that most impactors were instead rocky planetesimals left behind at 0.5-1.5 au after the terrestrial planet accretion. The number of basins expected from impacts of leftover planetesimals largely exceeds the number of known lunar basins, suggesting that the first 200 Myr of impacts is not recorded on the lunar surface. The Imbrium basin formation (age 3.92 Gyr; impactor diameter d~100 km) occurs with a 15-35% probability in our model. Imbrium must have formed unusually late to have only two smaller basins (Orientale and Schrodinger) forming afterwards. The model predicts 20 d>10-km impacts on the Earth 2.5-3.5 Gyr ago (Ga), which is comparable to the number of known spherule beds in the late Archean. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: ApJL, in press

The Transition from the Kuiper Belt to the Jupiter-Family (Comets)

Authors: Wesley C. Fraser, Luke Dones, Kathryn Volk, Maria Womack, David Nesvorný

Abstract: Kuiper Belt Objects, or more generally Trans-Neptunian Objects (TNOs), are planetesimals found beyond the orbit of Neptune. Some TNOs evolve onto Neptune-crossing orbits and become Centaurs. Many Centaurs, in turn, reach Jupiter-crossing orbits and become Jupiter-family comets (JFCs). TNOs are the main source of the JFCs. TNOs offer a different window than the JFCs, of more primordial bodies and o… ▽ More Kuiper Belt Objects, or more generally Trans-Neptunian Objects (TNOs), are planetesimals found beyond the orbit of Neptune. Some TNOs evolve onto Neptune-crossing orbits and become Centaurs. Many Centaurs, in turn, reach Jupiter-crossing orbits and become Jupiter-family comets (JFCs). TNOs are the main source of the JFCs. TNOs offer a different window than the JFCs, of more primordial bodies and over a different size and temperature range. It is in that context that this chapter is written. Here we discuss the dynamical pathways taken from the trans-Neptunian region to the JFCs, and the most important properties of TNOs that relate to the JFC population, including considerations of their origins, compositions, morphologies, and size distributions. We relate these properties to the JFCs whenever possible. We reflect on a few key outstanding issues regarding our incomplete knowledge of TNOs as they pertain to the Centaurs and JFC populations. We finish with a short discussion of notable new and upcoming facilities and the impacts they will have regarding these outstanding questions. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 32 pages, 2 tables, 15 figures. Accepted for publication

Implications of Jupiter Inward Gas-Driven Migration for the Inner Solar System

Authors: Rogerio Deienno, Andre Izidoro, Alessandro Morbidelli, David Nesvorny, William F. Bottke

Abstract: The migration history of Jupiter in the sun's natal disk remains poorly constrained. Here we consider how Jupiter's migration affects small-body reservoirs and how this constrains its original orbital distance from the Sun. We study the implications of large-scale and inward radial migration of Jupiter for the inner solar system while considering the effects of collisional evolution of planetesima… ▽ More The migration history of Jupiter in the sun's natal disk remains poorly constrained. Here we consider how Jupiter's migration affects small-body reservoirs and how this constrains its original orbital distance from the Sun. We study the implications of large-scale and inward radial migration of Jupiter for the inner solar system while considering the effects of collisional evolution of planetesimals. We use analytical prescriptions to simulate the growth and migration of Jupiter in the gas disk. We assume the existence of a planetesimal disk inside Jupiter's initial orbit. This planetesimal disk received an initial total mass and size-frequency distribution (SFD). Planetesimals feel the effects of aerodynamic gas drag and collide with one another, mostly while shepherded by the migrating Jupiter. Our main goal is to measure the amount of mass in planetesimals implanted into the main asteroid belt (MAB) and the SFD of the implanted population. We also monitor the amount of dust produced during planetesimal collisions. We find that the SFD of the planetesimal population implanted into the MAB tends to resemble that of the original planetesimal population interior to Jupiter. We also find that unless very little or no mass existed between 5 au and Jupiter's original orbit, it would be difficult to reconcile the current low mass of the MAB with the possibility that Jupiter migrated from distances beyond 15 au. This is because the fraction of the original disk mass that gets implanted into the MAB is very large. Finally, we discuss the implications of our results in terms of dust production to the so-called NC-CC isotopic dichotomy. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: Accepted for publication in The Astrophysical Journal Letters; In press

Implications for the collisional strength of Jupiter Trojans from the Eurybates family

Authors: Raphael Marschall, David Nesvorný, Rogerio Deienno, Ian Wong, Harold F. Levison, William F. Bottke

Abstract: In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated. The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD… ▽ More In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated. The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD) that remains largely unaltered at large sizes (>10 km) but is depleted at small sizes (10 m to 1 km). This results in a turnover in the SFD, the location of which depends on the collisional strength of the material. It is to be expected that the Trojan SFD bends between 1 and 10 km. Based on the SFD of the Eurybates family, we find that the family was likely the result of a catastrophic impact onto a 100 km rubble pile target. This corresponds to objects with a rather low collisional strength (10 times weaker than that of basaltic material studied in Benz & Asphaug1999). Assuming this weak strength, and an initial cumulative slope of the size frequency distribution of 2.1 between diameters of 2 m and 100 km when the Trojans were captured, the existence of Queta, the satellite of Eurybates, implies an upper limit for the family age of 3.7 Gy. Alternatively, we demonstrate that an unconventional collisional strength with a minimum at 20 m is a plausible candidate to ensure the survival of Queta over the age of the Solar System. Finally, we show how different collisional histories change the expected number of craters on the targets of the Lucy mission and that Lucy will be able to differentiate between them. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 19 pages, 13 figures

Trapping (sub-)Neptunes similar to TOI-216b at the inner disk rim: Implications for the disk viscosity and the Neptunian desert

Authors: Ondřej Chrenko, Raúl O. Chametla, David Nesvorný, Mario Flock

Abstract: [Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive co-rotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, lose the support of the co-rotation torque as their co-rotation regions… ▽ More [Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive co-rotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, lose the support of the co-rotation torque as their co-rotation regions become depleted, and the trapping efficiency then becomes uncertain. We study what it takes to trap such gap-opening planets at the inner disk rim. We performed 2D locally isothermal and non-isothermal hydrodynamic simulations of planet migration. A viscosity transition was introduced in the disk to (i) create a density drop and (ii) mimic the viscosity increase as the planet migrated from a dead zone towards a region with active magneto-rotational instability (MRI). We chose TOI-216b as a Neptune-like upper-limit test case, but we also explored different planetary masses, both on fixed and evolving orbits. For planet-to-star mass ratios $q\simeq(4$-$8)\times10^{-5}$, the density drop at the disk rim becomes reshaped due to a gap opening and is often replaced with a small density bump centered on the planet's corotation. Trapping is possible only if the bump retains enough gas mass and if the co-rotation region becomes azimuthally asymmetric, with an island of librating streamlines that accumulate a gas overdensity ahead of the planet. The overdensity exerts a positive torque that can counteract the negative torque of spiral arms. In our model, efficient trapping depends on the $α$ viscosity and its contrast across the viscosity transition. In order to trap TOI-216b, $α_{\mathrm{DZ}}=10^{-3}$ in the dead zone requires $α_{\mathrm{MRI}}\gtrsim5\times10^{-2}$ in the MRI-active zone. If $α_{\mathrm{DZ}}=5\times10^{-4}$, $α_{\mathrm{MRI}}\gtrsim7.5\times10^{-2}$ is needed. △ Less

Submitted 30 August, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A, 21 pages, 16 figures; v2 with minor text corrections

Journal ref: A&A 666, A63 (2022)

SPH simulations of high-speed collisions between asteroids and comets

Authors: J. Rozehnal, M. Brož, D. Nesvorný, K. J. Walsh, D. D. Durda, D. C. Richardson, E. Asphaug

Abstract: We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter… ▽ More We studied impact processes by means of smoothed-particle hydrodynamics (SPH) simulations. The method was applied to modeling formation of main-belt families during the cometary bombardment (either early or late, ${\sim}\,3.85\,{\rm Gy}$ ago). If asteroids were bombarded by comets, as predicted by the Nice model, hundreds of asteroid families (catastrophic disruptions of diameter $D \ge 100\,{\rm km}$ bodies) should have been created, but the observed number is only 20. Therefore we computed a standard set of 125 simulations of collisions between representative $D = 100\,{\rm km}$ asteroids and high-speed icy projectiles (comets). According to our results, the largest remnant mass $M_{\rm lr}$ is similar as in low-speed collisions, due to appropriate scaling with the effective strength $Q_{\rm eff}$, but the largest fragment mass $M_{\rm lf}$ exhibits systematic differences - it is typically smaller for craterings and bigger for super-catastrophic events. This trend does not, however, explain the non-existence of old families. The respective parametric relations can be used in other statistical (Monte-Carlo) models to better understand collisions between asteroidal and cometary populations. △ Less

Submitted 28 May, 2022; originally announced May 2022.

Journal ref: Icarus 383, 115064, 2022

HD 28109 hosts a trio of transiting Neptunian planets including a near-resonant pair, confirmed by ASTEP from Antarctica

Authors: Georgina Dransfield, Amaury H. M. J. Triaud, Tristan Guillot, Djamel Mekarnia, David Nesvorný, Nicolas Crouzet, Lyu Abe, Karim Agabi, Marco Buttu, Juan Cabrera, Davide Gandolfi, Maximilian N. Günther, Florian Rodler, François-Xavier Schmider, Philippe Stee, Olga Suarez, Karen A. Collins, Martín Dévora-Pajares, Steve B. Howell, Elisabeth C. Matthews, Matthew R. Standing, Keivan G. Stassun, Chris Stockdale, Samuel N. Quinn, Carl Ziegler , et al. (6 additional authors not shown)

Abstract: We report on the discovery and characterisation of three planets orbiting the F8 star HD~28109, which sits comfortably in \tess's continuous viewing zone. The two outer planets have periods of $\rm 56.0067 \pm 0.0003~days$ and $\rm 84.2597_{-0.0008}^{+0.0010}~days$, which implies a period ratio very close to that of the first-order 3:2 mean motion resonance, exciting transit timing variations (TTV… ▽ More We report on the discovery and characterisation of three planets orbiting the F8 star HD~28109, which sits comfortably in \tess's continuous viewing zone. The two outer planets have periods of $\rm 56.0067 \pm 0.0003~days$ and $\rm 84.2597_{-0.0008}^{+0.0010}~days$, which implies a period ratio very close to that of the first-order 3:2 mean motion resonance, exciting transit timing variations (TTVs) of up to $\rm 60\,mins$. These two planets were first identified by \tess, and we identified a third planet in the \textcolor{black}{\tess photometry} with a period of $\rm 22.8911 \pm 0.0004~days$. We confirm the planetary nature of all three planetary candidates using ground-based photometry from Hazelwood, ASTEP and LCO, including a full detection of the $\rm \sim9\,h$ transit of HD~28109 c from Antarctica. The radii of the three planets are \textcolor{black}{$\rm R_b=2.199_{-0.10}^{+0.098} ~R_{\oplus}$, $\rm R_c=4.23\pm0.11~ R_{\oplus}$ and $\rm R_d=3.25\pm0.11 ~R_{\oplus}$}; we characterise their masses using TTVs and precise radial velocities from ESPRESSO and HARPS, and find them to be $\rm M_b=18.5_{-7.6}^{+9.1}~M_{\oplus}$, $\rm M_c=7.9_{-3.0}^{+4.2}~M_{\oplus}$ and $\rm M_d=5.7_{-2.1}^{+2.7}~M_{\oplus}$, making planet b a dense, massive planet while c and d are both under-dense. We also demonstrate that the two outer planets are ripe for atmospheric characterisation using transmission spectroscopy, especially given their position in the CVZ of JWST. The data obtained to date are consistent with resonant (librating) and non-resonant (circulating) solutions; additional observations will show whether the pair is actually locked in resonance or just near-resonant. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 19 pages, 14 figures

Asteroid Families: properties, recent advances and future opportunities

Authors: Bojan Novakovic, David Vokrouhlicky, Federica Spoto, David Nesvorny

Abstract: Collisions are one of the key processes shaping planetary systems. Asteroid families are outcomes of such collisions still identifiable across our solar system. The families provide a unique view of catastrophic disruption phenomena and have been in the focus of planetary scientists for more than a century. Most of them are located in the main belt, a ring of asteroids between Mars and Jupiter. He… ▽ More Collisions are one of the key processes shaping planetary systems. Asteroid families are outcomes of such collisions still identifiable across our solar system. The families provide a unique view of catastrophic disruption phenomena and have been in the focus of planetary scientists for more than a century. Most of them are located in the main belt, a ring of asteroids between Mars and Jupiter. Here we review the basic properties of the families, discuss some recent advances, and anticipate future challenges. This review pays more attention to dynamic aspects such as family identification, age determination, and long-term evolution. The text, however, goes beyond that. Especially, we cover the details of young families that see the major advances in the last years, and we anticipate it will develop even faster in the future. We also discuss the relevance of asteroid families for water-ice content in the asteroid belt and our current knowledge on links between families and main-belt comets. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: Review paper to appear in CeMDA's topical collection on "Main Belt Dynamics"

Journal ref: Celest Mech Dyn Astron 134, 34 (2022)

Exogenous delivery of water to Mercury

Authors: Kateryna Frantseva, David Nesvorný, Michael Mueller, Floris F. S. van der Tak, Inge Loes ten Kate, Petr Pokorný

Abstract: Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impa… ▽ More Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impacts with Mercury. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We add comets and asteroids to the simulations as massless test particles, based on their current orbital distributions. Asteroid impactors are assigned a probability of being water-rich (C-class) based on the measured distribution of taxonomic types. For comets, we assume a constant water fraction. For IDPs, we use a dynamical meteoroid model to compute the dust flux on Mercury. Relative to previous work on asteroid and comet impacts (Moses et al. 1999), we leverage 20 years of progress in minor body surveys. Immediate post-impact ejection of impactor material into outer space is taken into account as is the migration efficiency of water across Mercury's surface to the polar cold traps. We find that asteroids deliver $\sim 1 \times 10^{3}$ kg/yr of water to Mercury, comets deliver $\sim 1 \times 10^{3}$ kg/yr and IDPs deliver $\sim 16 \times 10^{3}$ kg/yr within a factor of several. Over a timescale of $\sim 1$ Gyr, this is enough to deliver the minimum amount of water required by the radar and MESSENGER observations. While other sources of water on Mercury are not ruled out by our analysis, we show that they are not required to explain the currently available observational lower limits. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 38 pages, 10 figures, published in Icarus

Thermal processing of Jupiter Family Comets during their chaotic orbital evolution

Authors: Anastasios Gkotsinas, Aurélie Guilbert-Lepoutre, Sean N. Raymond, David Nesvorný

Abstract: Evidence for cometary activity beyond Jupiter and Saturn's orbits -- such as that observed for Centaurs and long period comets -- suggests that the thermal processing of comet nuclei starts long before they enter the inner Solar System, where they are typically observed and monitored. Such observations raise questions as to the depth of unprocessed material, and whether the activity of JFCs can be… ▽ More Evidence for cometary activity beyond Jupiter and Saturn's orbits -- such as that observed for Centaurs and long period comets -- suggests that the thermal processing of comet nuclei starts long before they enter the inner Solar System, where they are typically observed and monitored. Such observations raise questions as to the depth of unprocessed material, and whether the activity of JFCs can be representative of any primitive material. Here we model the coupled thermal and dynamical evolution of Jupiter Family Comets (JFCs), from the moment they leave their outer Solar System reservoirs until their ejection into interstellar space. We apply a thermal evolution model to a sample of simulated JFCs obtained from dynamical simulations (arXiv:1706.07447) that successfully reproduce the orbital distribution of observed JFCs. We show that due to the stochastic nature of comet trajectories toward the inner solar system, all simulated JFCs undergo multiple heating episodes resulting in significant modifications of their initial volatile contents. A statistical analysis constrains the extent of such processing. We suggest that primordial condensed hypervolatile ices should be entirely lost from the layers that contribute to cometary activity observed today. Our results demonstrate that understanding the orbital (and thus, heating) history of JFCs is essential when putting observations in a broader context. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 30 pages, 10 figures, to be published in ApJ

Dynamical Implantation of Blue Binaries in the Cold Classical Kuiper Belt

Authors: David Nesvorny, David Vokrouhlicky, Wesley C. Fraser

Abstract: Colors and binarity provide important constraints on the Kuiper belt formation. The cold classical objects at radial distance r=42-47 au from the Sun are predominantly very red (spectral slope s>17%) and often exist as equal-size binaries (~30% observed binary fraction). This has been taken as evidence for the in-situ formation of cold classicals. Interestingly, a small fraction (~10%) of cold cla… ▽ More Colors and binarity provide important constraints on the Kuiper belt formation. The cold classical objects at radial distance r=42-47 au from the Sun are predominantly very red (spectral slope s>17%) and often exist as equal-size binaries (~30% observed binary fraction). This has been taken as evidence for the in-situ formation of cold classicals. Interestingly, a small fraction (~10%) of cold classicals is less red with s<17%, and these "blue" bodies are often found in wide binaries. Here we study the dynamical implantation of blue binaries from r<42 au. We find that they can be implanted into the cold classical belt from a wide range of initial radial distances, but the survival of the widest blue binaries -- 2001 QW322 and 2003 UN284 -- implies formation at r>30 au. This would be consistent with the hypothesized less-red to very-red transition at 30<r<40 au. For any reasonable choice of parameters (Neptune's migration history, initial disk profile, etc.), however, our model predicts a predominance of blue singles, rather than blue binaries, which contradicts existing observations. We suggest that wide blue binaries formed in situ at r=42-47 au and their color reflects early formation in a protoplanetary gas disk. The predominantly VR colors of cold classicals may be related to the production of methanol and other hydrocarbons during the late stages of the disk, when the temperature at 45 au dropped to 20 K and carbon monoxide was hydrogenated. △ Less

Submitted 7 January, 2022; originally announced January 2022.

Comments: AJ

Col-OSSOS: Probing Ice Line/Color Transitions within the Kuiper Belt's Progenitor Populations

Authors: Laura E. Buchanan, Megan E. Schwamb, Wesley C. Fraser, Michele T. Bannister, Michaël Marsset, Rosemary E. Pike, David Nesvorný, J. J. Kavelaars, Susan D. Benecchi, Matthew J. Lehner, Shiang-Yu Wang, Nuno Peixinho, Kathryn Volk, Mike Alexandersen, Ying-Tung Chen, Brett Gladman, Stephen Gwyn, Jean-Marc Petit

Abstract: Dynamically excited objects within the Kuiper belt show a bimodal distribution in their surface colors, and these differing surface colors may be a tracer of where these objects formed. In this work we explore radial color distributions in the primordial planetesimal disk and implications for the positions of ice line/color transitions within the Kuiper belt's progenitor populations. We combine a… ▽ More Dynamically excited objects within the Kuiper belt show a bimodal distribution in their surface colors, and these differing surface colors may be a tracer of where these objects formed. In this work we explore radial color distributions in the primordial planetesimal disk and implications for the positions of ice line/color transitions within the Kuiper belt's progenitor populations. We combine a full dynamical model of the Kuiper belt's evolution due to Neptune's migration with precise surface colors measured by the Colours of the Outer Solar System Origins Survey in order to examine the true color ratios within the Kuiper belt and the ice lines within the primordial disk. We investigate the position of a dominant, surface color changing ice-line, with two possible surface color layouts within the initial disk; (1) inner neutral surfaces and outer red, and (2) inner red surfaces and outer neutral. We performed simulations with a primordial disk that truncates at 30 au. By radially stepping the color transition out through 0.5 au intervals we show that both disk configurations are consistent with the observed color fraction. For an inner neutral, outer red primordial disk we find that the color transition can be at $28^{+2}_{-3}$ au at a 95% confidence level. For an inner red, outer neutral primordial disk the color transition can be at $27^{+3}_{-3}$ au at a 95% confidence level. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: 18 pages, 8 figures, accepted for publication in PSJ

The Stability Boundary of the Distant Scattered Disk

Authors: Konstantin Batygin, Rosemary A. Mardling, David Nesvorny

Abstract: The scattered disk is a vast population of trans-Neptunian minor bodies that orbit the sun on highly elongated, long-period orbits. The stability of scattered disk objects is primarily controlled by a single parameter - their perihelion distance. While the existence of a perihelion boundary that separates chaotic and regular motion of long-period orbits is well established through numerical experi… ▽ More The scattered disk is a vast population of trans-Neptunian minor bodies that orbit the sun on highly elongated, long-period orbits. The stability of scattered disk objects is primarily controlled by a single parameter - their perihelion distance. While the existence of a perihelion boundary that separates chaotic and regular motion of long-period orbits is well established through numerical experiments, its theoretical basis as well as its semi-major axis dependence remain poorly understood. In this work, we outline an analytical model for the dynamics of distant trans-Neptunian objects and show that the orbital architecture of the scattered disk is shaped by an infinite chain of $2:j$ resonances with Neptune. The widths of these resonances increase as the perihelion distance approaches Neptune's semi-major axis, and their overlap drives chaotic motion. Within the context of this picture, we derive an analytic criterion for instability of long-period orbits, and demonstrate that rapid dynamical chaos ensues when the perihelion drops below a critical value, given by $q_{\rm{crit}}=a_{\rm{N}}\,\big(\ln((24^2/5)\,(m_{\rm{N}}/M_{\odot})\,(a/a_{\rm{N}})^{5/2})\big)^{1/2}$. This expression constitutes a boundary between the "detached" and actively "scattering" sub-populations of distant trans-Neptunian minor bodies. Additionally, we find that within the stochastic layer, the Lyapunov time of scattered disk objects approaches the orbital period, and show that the semi-major axis diffusion coefficient is approximated by $\mathcal{D}_a\sim(8/(5\,π))\,(m_{\rm{N}}/M_{\odot})\,\sqrt{\mathcal{G}\,M_{\odot}\,a_{\rm{N}}}\,\exp\big[-(q/a_{\rm{N}})^2/2\big]$. We confirm our results with numerical simulations and highlight the connections between scattered disk dynamics and the Chirikov Standard Map. Implications of our results for the long-term evolution of the distant solar system are discussed. △ Less

Submitted 30 October, 2021; originally announced November 2021.

Comments: 14 pages, 3 figures, published in ApJ

Journal ref: The Astrophysical Journal, 2021, Volume 920, Number 2

TOI-216: Resonant Constraints on Planet Migration

Authors: David Nesvorny, Ondrej Chrenko, Mario Flock

Abstract: TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_ψ\simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit.… ▽ More TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_ψ\simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit. We carry out $N$-body simulations of planet migration in a protoplanetary gas disk to explain the orbital configuration of TOI-216 planets. We find that TOI-216b's migration must have been halted near its current orbital radius to allow for a convergent migration of the two planets into the resonance. For the inferred damping-to-migration timescale ratio $τ_e/τ_a \simeq 0.02$, overstable librations in the resonance lead to a limit cycle with $A_ψ\simeq 80^\circ$ and $e_{\rm b}<0.1$. The system could have remained in this configuration for the greater part of the protoplanetary disk lifetime. If the gas disk was removed from inside out, this would have reduced the libration amplitude to $A_ψ\simeq 60^\circ$ and boosted $e_{\rm b}$ via the resonant interaction with TOI-216c. Our results suggest a relatively fast inner disk removal ($\sim 10^5$ yr). Another means of explaining the large libration amplitude is stochastic stirring from a (turbulent) gas disk. For that to work, overstable librations would need to be suppressed, $τ_e/τ_a \simeq 0.05$, and very strong turbulent stirring (or some other source of large stochastic forcing) would need to overcome the damping effects of gas. Hydrodynamical simulations can be performed to test these models. △ Less

Submitted 18 October, 2021; originally announced October 2021.

Comments: ApJ

iSyMBA: A Symplectic Massive Bodies Integrator with Planets Interpolation

Authors: Fernando Roig, David Nesvorný, Rogerio Deienno, Matias J. Garcia

Abstract: A planetary instability occurring at time $<100$ My after formation of the giant planets in our solar system can be responsible for some characteristics of the inner solar system. However, the actual influence of the instability on the terrestrial planet formation is not well understood. The simulations of terrestrial planet formation are very CPU-expensive, and this limits the exploration of diff… ▽ More A planetary instability occurring at time $<100$ My after formation of the giant planets in our solar system can be responsible for some characteristics of the inner solar system. However, the actual influence of the instability on the terrestrial planet formation is not well understood. The simulations of terrestrial planet formation are very CPU-expensive, and this limits the exploration of different instability scenarios. To include the effects of the giant planets instability in the simulations of terrestrial planets formation in a feasible way, we approach the problem in two steps. First, we model and record an evolution of the giant planets that replicates the present outer solar system in the end. Then, we use that orbital record, properly interpolated, as the input for a second step to simulate its effects on the terrestrial planet formation. For this second step, we developed iSyMBA, a symplectic massive bodies algorithm, where ``i'' stands for interpolation. iSyMBA is a very useful code to accurately evaluate the effects of planetary instabilities on minor body reservoirs, while accounting for close encounters among massive objects. We provide a detailed description of how iSyMBA was developed and implemented to study terrestrial planet formation. Adapting iSyMBA for other problems that demand interpolation from previous simulations can be done following the method described here. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Terrestrial planet formation by torque-driven convergent migration of planetary embryos

Authors: M. Brož, O. Chrenko, D. Nesvorný, N. Dauphas

Abstract: Massive cores of the giant planets are thought to have formed in a gas disk by accretion of pebble-size particles whose accretional cross-section is enhanced by aerodynamic gas drag [1][2]. A commonly held view is that the terrestrial planet system formed later (30-200 Myr after the dispersal of the gas disk) by giant collisions of tens of roughly Mars-size protoplanets [3]. Here we propose, inste… ▽ More Massive cores of the giant planets are thought to have formed in a gas disk by accretion of pebble-size particles whose accretional cross-section is enhanced by aerodynamic gas drag [1][2]. A commonly held view is that the terrestrial planet system formed later (30-200 Myr after the dispersal of the gas disk) by giant collisions of tens of roughly Mars-size protoplanets [3]. Here we propose, instead, that the terrestrial planets formed earlier by gas-driven convergent migration of protoplanets toward $\sim\!1\,{\rm au}$ (related ref. [4] invoked a different process to concentrate planetesimals). To investigate situations in which convergent migration occurs, we developed a radiation-hydrodynamic model with realistic opacities [5][6] to determine the thermal structure of the gas and pebble disks in the terrestrial planet zone. We find that protoplanets rapidly grow by mutual collisions and pebble accretion, and gain orbital eccentricities by gravitational scattering and the hot-trail effect [7][8]. The orbital structure of the terrestrial planet system is well reproduced in our simulations, including its tight mass concentration at 0.7-1 au and the small sizes of Mercury and Mars. The early-stage protosolar disk temperature exceeds 1500 K inside 0.4 au implying that Mercury grew in a highly reducing environment, next to the evaporation lines of iron and silicates, influencing Mercury's bulk composition [9]. In a late-stage cold gas disk, accretion of icy/hydrated pebbles would contribute to Earth's water budget. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Comments: original submission; final vers. is available at NA website

Journal ref: Nature Astronomy, 5, 898-902, 2021

A pair of warm giant planets near the 2:1 mean motion resonance around the K-dwarf star TOI-2202

Authors: Trifon Trifonov, Rafael Brahm, Nestor Espinoza, Thomas Henning, Andrés Jordán, David Nesvorny, Rebekah I. Dawson, Jack J. Lissauer, Man Hoi Lee, Diana Kossakowski, Felipe I. Rojas, Melissa J. Hobson, Paula Sarkis, Martin Schlecker, Bertram Bitsch, Gaspar Á. Bakos, Mauro Barbieri, Waqas Bhatti, R. Paul Butler, Jeffrey D. Crane, Sangeetha Nandakumar, Matías R. Díaz, Stephen Shectman, Johanna Teske, Pascal Torres , et al. (15 additional authors not shown)

Abstract: TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P=11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hours. Radial velocity follo… ▽ More TOI-2202 b is a transiting warm Jovian-mass planet with an orbital period of P=11.91 days identified from the Full Frame Images data of five different sectors of the TESS mission. Ten TESS transits of TOI-2202 b combined with three follow-up light curves obtained with the CHAT robotic telescope show strong transit timing variations (TTVs) with an amplitude of about 1.2 hours. Radial velocity follow-up with FEROS, HARPS and PFS confirms the planetary nature of the transiting candidate (a$_{\rm b}$ = 0.096 $\pm$ 0.002 au, m$_{\rm b}$ = 0.98 $\pm$ 0.06 M$_{\rm Jup}$), and dynamical analysis of RVs, transit data, and TTVs points to an outer Saturn-mass companion (a$_{\rm c}$ = 0.155 $\pm$ 0.003 au, m$_{\rm c}$= $0.37 \pm 0.10$ M$_{\rm Jup}$) near the 2:1 mean motion resonance. Our stellar modeling indicates that TOI-2202 is an early K-type star with a mass of 0.82 M$_\odot$, a radius of 0.79 R$_\odot$, and solar-like metallicity. The TOI-2202 system is very interesting because of the two warm Jovian-mass planets near the 2:1 MMR, which is a rare configuration, and their formation and dynamical evolution are still not well understood. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Accepted for publication in AJ

The young Hobson family: Possible binary parent body and low-velocity dispersal

Authors: David Vokrouhlický, Miroslav Brož, Bojan Novaković, David Nesvorný

Abstract: Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on th… ▽ More Asteroid families with ages younger than $1$ Myr offer an interesting possibility of studying the outcomes of asteroid disruptions that are little modified by subsequent evolutionary processes. We analyze a very young asteroid family associated with (18777) Hobson in the central part of the main belt. We aim at (i) understanding its peculiar size distribution, and (ii) setting an upper limit on the characteristic dispersal velocity at subkilometer sizes corresponding to the smallest visible Hobson members. We identified the Hobson family using an up-to-date asteroid catalog. A significant increase in the number of its known members allowed us to study their size distribution and compare it with computer simulations of catastrophic disruptions. Backward orbital integrations of the heliocentric orbits allowed us to confirm the previously suggested age of Hobson and helped to estimate limits of the ejection speed. The Hobson family has an unusual size distribution: two nearly equal-size bodies, followed by a population of smaller asteroids, whose distribution takes a characteristic power law. There are two possibilities to explain these data. Either a canonical impact onto a single parent body, requiring fine-tuned impact conditions that have not been studied so far, or an unconventional model for the parent object of the Hobson family, namely a binary with $\simeq 7-9$ km primary and a $\simeq 2.5$ km secondary. In the latter case, the primary was disrupted, leaving behind the largest remnant (18777) Hobson and a suite of subkilometer asteroids. The second largest asteroid, (57738) 2001 UZ160, is the nearly intact satellite of the parent binary. The excellent convergence of nominal orbits of Hobson members sets an upper limit of $\simeq (10-20)$ m s$^{-1}$ for the initial dispersal velocity of the known members, which is consistent with both formation models. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: 17 pages, 12 figures, 1 table, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 654, A75 (2021)

Can a jumping-Jupiter trigger the Moon's formation impact?

Authors: Sandro R. DeSouza, Fernando Roig, David Nesvorný

Abstract: We investigate the possibility that the Moon's formation impact was triggered by an early dynamical instability of the giant planets. We consider the well-studied "jumping Jupiter" hypothesis for the solar system's instability, where Jupiter and Saturn's semi-major axes evolve in step-wise manner from their primordially compact architecture to their present locations. Moreover, we test multiple di… ▽ More We investigate the possibility that the Moon's formation impact was triggered by an early dynamical instability of the giant planets. We consider the well-studied "jumping Jupiter" hypothesis for the solar system's instability, where Jupiter and Saturn's semi-major axes evolve in step-wise manner from their primordially compact architecture to their present locations. Moreover, we test multiple different configurations for the primordial system of terrestrial planets and the Moon-forming projectile, with particular focus on the almost equal masses impact. We find that the instability/migration of the giant planets excites the orbits of the terrestrial planets through dynamical perturbations, thus allowing collisions between them. About 10% of the simulations lead to a collision with the proto-Earth which resulted in a final configuration of the terrestrial system that reproduces, to some extent, its present architecture. Most of these collisions occur in the hit-and-run domain, but about 15% occur in the partial accretion regime, with the right conditions for a Moon-forming impact. In most of the simulations, there is a delay of more than ~20 My between the time of the instability and the Moon-forming impact. This supports the occurrence of an early instability (< 10 My} after dissipation of the gas in the proto-planetary disk), compatible with the time of the Moon-forming impact (30-60 My) inferred from cosmochemical constraints. In general, the final states of the inner solar system in our simulations show an excess of Angular Momentum Deficit, mostly attributed to the over-excitation of Mercury's eccentricity and inclination. △ Less

Submitted 14 July, 2021; v1 submitted 8 July, 2021; originally announced July 2021.

Dark Primitive Asteroids Account for a Large Share of K/Pg-Scale Impacts on the Earth

Authors: David Nesvorny, William F. Bottke, Simone Marchi

Abstract: A dynamical model for large near-Earth asteroids (NEAs) is developed here to understand the occurrence rate and nature of Cretaceous-Paleogene (K/Pg) scale impacts on the Earth. We find that 16--32 (2--4) impacts of diameter $D>5$ km ($D>10$ km) NEAs are expected on the Earth in 1 Gyr, with about a half of impactors being dark primitive asteroids (most of which start with semimajor axis $a>2.5$ au… ▽ More A dynamical model for large near-Earth asteroids (NEAs) is developed here to understand the occurrence rate and nature of Cretaceous-Paleogene (K/Pg) scale impacts on the Earth. We find that 16--32 (2--4) impacts of diameter $D>5$ km ($D>10$ km) NEAs are expected on the Earth in 1 Gyr, with about a half of impactors being dark primitive asteroids (most of which start with semimajor axis $a>2.5$ au). These results explain why the Chicxulub crater, the third largest impact structure found on the Earth (diameter $\simeq180$ km), was produced by impact of a carbonaceous chondrite. They suggest, when combined with previously published results for small ($D \lesssim 1$ km) NEAs, a size-dependent sampling of the main belt. We conclude that the impactor that triggered the K/Pg mass extinction $\simeq 66$ Myr ago was a main belt asteroid that quite likely ($\simeq 60$\% probability) originated beyond 2.5 au. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: Icarus

Col-OSSOS: The Distinct Colour Distribution of Single and Binary Cold Classical KBOs

Authors: Wesley C. Fraser, Susan D. Benecchi, JJ Kavelaars, Michael Marsset, Rosemary Pike, Michele T. Bannister, Megan E. Schwamb, Kathryn Volk, David Nesvorny, Mike Alexandersen, Ying-Tung Chen, Stephen Gwyn, Matthew J. Lehner, Shiang-Yu Wang

Abstract: The cold classical Kuiper Belt Objects (KBOs) possess a high, $\gtrsim30\%$ binary fraction. Widely separated and dynamically fragile, these binary systems have been useful in tracing the origins of KBOs. A new class of binaries was recently identified by their colours. The so-called blue binaries are unanimously members of the less red compositional class, and exhibit a 100% binary fraction. They… ▽ More The cold classical Kuiper Belt Objects (KBOs) possess a high, $\gtrsim30\%$ binary fraction. Widely separated and dynamically fragile, these binary systems have been useful in tracing the origins of KBOs. A new class of binaries was recently identified by their colours. The so-called blue binaries are unanimously members of the less red compositional class, and exhibit a 100% binary fraction. They appear to be push-out survivors, emplaced in the classical region during Neptune's phases of outward migration. The presence of these binary systems implies that the majority of objects that formed near the cold classical region formed as binaries. Here we present new optical colour measurements of cold classical KBOs from the Colours of the Outer Solar System Origins Survey, including colours of a blue binary discovered by the Solar System Origins Legacy Survey -- 2015 RJ277. The increased size of the colours sample has resulted in order-of-magnitude decrease in the probability that the binaries and singles sample share the same colour distribution. From the Anderson-Darling statistic, this probability is only a 0.3%, while it is only 0.002% when utilizing the difference of means statistic. We find a hint that the blue binaries have inflated free inclinations compared to their red counterparts, consistent with the push-out origin for these bodies. △ Less

Submitted 31 March, 2021; originally announced April 2021.

Comments: Accepted for publication in the PSJ. 4 Figures