Showing 1–50 of 200 results for author: Morbidelli, A

The Solar System could have formed in a low-viscosity disc: A dynamical study from giant planet migration to the Nice model

Authors: Philippine Griveaud, Aurélien Crida, Antoine C. Petit, Elena Lega, Alessandro Morbidelli

Abstract: In the context of low-viscosity protoplanetary discs (PPDs), the formation scenarios of the Solar System should be revisited. In particular, the Jupiter-Saturn pair has been shown to lock in the 2:1 mean motion resonance while migrating generally inwards, making the Grand Tack scenario impossible. We explore what resonant chains of multiple giant planets can form in a low-viscosity disc, and wheth… ▽ More In the context of low-viscosity protoplanetary discs (PPDs), the formation scenarios of the Solar System should be revisited. In particular, the Jupiter-Saturn pair has been shown to lock in the 2:1 mean motion resonance while migrating generally inwards, making the Grand Tack scenario impossible. We explore what resonant chains of multiple giant planets can form in a low-viscosity disc, and whether these configurations can evolve into forming the Solar System in the post gas disc phase. We used hydrodynamical simulations to study the migration of the giant planets in a disc with viscosity $α=10^{-4}$. After a transition phase to a gas-less configuration, we studied the stability of the obtained resonant chains through their interactions with a disc of leftover planetesimals by performing N-body simulations using rebound. The gaps open by giant planets are wider and deeper for lower viscosity, reducing the damping effect of the disc and thus weakening resonant chains. Exploring numerous configurations, we found five stable resonant chains of four or five planets. In a thin PPD, the four giant planets revert their migration and migrate outwards. After disc dispersal, under the influence of a belt of planetesimals, some resonant chains undergo an instability phase while others migrate smoothly over a billion years. For three of our resonant chains, about 1% of the final configurations pass the four criteria to fit the Solar System. The most successful runs are obtained for systems formed in a cold PPD with a massive planetesimal disc. This work provides a fully consistent study of the dynamical history of the Solar System's giant planets, from the protoplanetary disc phase up to the giant planet instability. Although building resonant configurations is difficult in low-viscosity discs, we find it possible to reproduce the Solar System from a cold, low-viscosity protoplanetary disc. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 13 Pages, 13 Figures, Accepted for publications in A&A

Forming the Trappist-1 system in two steps during the recession of the disc inner edge

Authors: Gabriele Pichierri, Alessandro Morbidelli, Konstantin Batygin, Ramon Brasser

Abstract: Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred o… ▽ More Trappist-1 hosts 7 planets where period ratios of neighbouring pairs are close to the 8:5, 5:3, 3:2, 3:2, 4:3, and 3:2 ratios in increasing distance from the star. The Laplace angles associated with neighbouring triplets are observed to be librating, proving the resonant nature of the system. This compact, resonant configuration is a manifest sign of disc-driven migration; however, the preferred outcome of such evolution is the establishment of first-order resonances, not the high-order resonances observed in the inner system. Here, we explain the observed orbital configuration in a model that is largely independent on the specific disc migration and orbital circularisation efficiencies. Together with migration, the two key elements of our model are: i) the inner border of the protoplanetary disc receded with time; and ii) the system was initially separated in two sub-systems. Specifically, the inner b, c d and e planets were initially placed in a 3:2 resonance chain and then evolved to the 8:5 -- 5:3 commensurability between planets b, c and d under the effect of the recession of the inner edge of the disc, whereas the outer planets migrated to the inner edge at a later time, establishing the remaining resonances. Our results pivot on the dynamical role of the presently unobservable recession of the inner edge of protoplanetary discs. They also reveal the role of recurring phases of convergent migration followed by resonant repulsion with associated orbital circularisation when resonant chains interact with migration barriers. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: To appear in Nature Astronomy

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

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

Forming Mercury from excited initial conditions

Authors: Jennifer Scora, Diana Valencia, Alessandro Morbidelli, Seth Jacobson

Abstract: Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found… ▽ More Mercury is notoriously difficult to form in solar system simulations, due to its small mass and iron-rich composition. Smooth particle hydrodynamics simulations of collisions have found that a Mercury-like body could be formed by one or multiple giant impacts, but due to the chaotic nature of collisions it is difficult to create a scenario where such impacts will take place. Recent work has found more success forming Mercury analogues by adding additional embryos near Mercury's orbit. In this work, we aim to form Mercury by simulating the formation of the solar system in the presence of the giant planets Jupiter and Saturn. We test out the effect of an inner disk of embryos added on to the commonly-used narrow annulus of initial material. We form Mercury analogues with core-mass fractions (CMF) $> 0.4$ in $\sim 10\%$ of our simulations, and twice that number of Mercury analogues form during the formation process, but are unstable and do not last to the end of the simulations. Mercury analogues form at similar rates for both disks with and without an inner component, and most of our Mercury analogues have lower CMF than that of Mercury, $\sim 0.7$, due to significant accretion of debris material. We suggest that a more in-depth understanding of the fraction of debris mass that is lost to collisional grinding is necessary to understand Mercury's formation, or some additional mechanism is required to stop this debris from accreting. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 15 pages, accepted for publication in ApJ

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

Dynamics of the Beta Pictoris planetary system and its falling evaporating bodies

Authors: H. Beust, J. Milli, A. Morbidelli, S. Lacour, A. -M. Lagrange, G. Chauvin, M. Bonnefoy, J. Wang

Abstract: For decades, the spectral variations of Beta Pictoris have been modelled as the result of the evaporation of exocomets close to the star, termed falling evaporating bodies (FEBs). Resonant perturbations by a giant planet have been proposed to explain the dynamical origin of these stargrazers. The disk is now known to harbour two giant planets, Beta Pic b and c, orbiting the star at 9.9 au and 2.7… ▽ More For decades, the spectral variations of Beta Pictoris have been modelled as the result of the evaporation of exocomets close to the star, termed falling evaporating bodies (FEBs). Resonant perturbations by a giant planet have been proposed to explain the dynamical origin of these stargrazers. The disk is now known to harbour two giant planets, Beta Pic b and c, orbiting the star at 9.9 au and 2.7 au. While the former almost matches the planet formerly suspected, the discovery of the latter complicates the picture. We first question the stability of the two-planet system. Then we investigate the dynamics of a disk of planetesimals orbiting the star with both planets to check the validity of the FEB generation mechanism. Symplectic N-body simulations are used to determine which regions of the planetesimal disk are dynamically stable. Then we focus on regions where disk particles are able to reach high eccentricities thanks to resonant mechanisms. The first result is that the system is dynamically stable. Both planets may temporarily fall in 7:1 mean motion resonance (MMR). Then, simulations reveal that the whole region extending between ~1.5 au and ~25 au is unstable to planetary perturbations. However, a disk below 1.5 au survives, which appears to constitute an active source of FEBs via high-order MMRs with Beta Pic c. Beta Pic b acts as a distant perturber that helps sustain the whole process. These simulations rule out the preceding FEB generation mechanism model, which placed their origin at around 4-5 au. Conversely, FEBs are likely to originate from a region much further in and related to MMRs with Beta Pic c. That mechanism also appears to last longer, as new planetesimals are able to continuously enter the MMRs and evolve towards the FEB state. Subsequently, the physical nature of the FEBs may differ from that previously thought, and presumably may not be icy. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: 14 pages, 10 figures. To appear in Astronomy & Astrophysics

A super-massive Neptune-sized planet

Authors: L. Naponiello, L. Mancini, A. Sozzetti, A. S. Bonomo, A. Morbidelli, J. Dou, L. Zeng, Z. M. Leinhardt, K. Biazzo, P. Cubillos, M. Pinamonti, D. Locci, A. Maggio, M. Damasso, A. F. Lanza, J. J. Lissauer, A. Bignamini, W. Boschin, L. G. Bouma, P. J. Carter, D. R. Ciardi, K. A. Collins, R. Cosentino, I. Crossfield, S. Desidera , et al. (33 additional authors not shown)

Abstract: Neptune-sized planets exhibit a wide range of compositions and densities, depending onf cators related to their formation and evolution history, such as the distance from their host stars and atmospheric escape processes. They can vary from relatively low-density planets with thick hydrogen-helium atmospheres to higher-density planets with a substantial amount of water or a rocky interior with a t… ▽ More Neptune-sized planets exhibit a wide range of compositions and densities, depending onf cators related to their formation and evolution history, such as the distance from their host stars and atmospheric escape processes. They can vary from relatively low-density planets with thick hydrogen-helium atmospheres to higher-density planets with a substantial amount of water or a rocky interior with a thinner atmosphere, such as HD 95338 b, TOI-849 b and TOI-2196 b. The discovery of exoplanets in the hot-Neptune desert, a region close to the host stars with a deficit of Neptune-sized planets, provides insights into the formation and evolution of planetary systems, including the existence of this region itself. Here we show observations of the transiting planet TOI-1853 b, which has a radius of 3.46 +- 0.08 Earth radii and orbits a dwarf star every 1.24 days. This planet has a mass of 73.2 +- 2.7 Earth masses, almost twice that of any other Neptune-sized planet known so far, and a density of 9.7 +- 0.8 grams per cubic centimetre. These values place TOI-1853 b in the middle of the Neptunian desert and imply that heavy elements dominate its mass. The properties of TOI-1853 b present a puzzle for conventional theories of planetary formation and evolution, and could be the result of several proto-planet collisions or the final state of an initially high-eccentricity planet that migrated closer to its parent star. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: Preprint submitted to Nature. Please refer to the published version for the final parameters estimations

Journal ref: Nature, published online on the 30th of August 2023 and printed the 12th of October 2023

An improved model of metal/silicate differentiation during Earth's accretion

Authors: K. I. Dale, D. C. Rubie, M. Nakajima, S. Jacobson, G. Nathan, G. J. Golabek, S. Cambioni, A. Morbidelli

Abstract: We improved the algorithm presented in Rubie et al. (2015) to model the chemical evolution of Earth driven by iron/silicate differentiation during the planet's accretion. The pressure at which the equilibration occurs during a giant impact is no longer a free parameter but is determined by the smooth particle hydrodynamic (SPH) simulations of Nakajima et al. (2021). Moreover, impacting planetesima… ▽ More We improved the algorithm presented in Rubie et al. (2015) to model the chemical evolution of Earth driven by iron/silicate differentiation during the planet's accretion. The pressure at which the equilibration occurs during a giant impact is no longer a free parameter but is determined by the smooth particle hydrodynamic (SPH) simulations of Nakajima et al. (2021). Moreover, impacting planetesimals are now assumed to be too small to cause melting and differentiation and thus their materials are stored in the crystalline upper mantle of the growing planet until a hydrostatically relaxed global magma ocean forms in the aftermath of a giant impact, whose depth is also estimated from Nakajima et al. (2021). With these changes, not all dynamical simulations lead to a satisfactory reproduction of the chemical composition of the bulk silicate Earth (BSE). Thus, the latter becomes diagnostic of the success of dynamical models. In the successful cases also the BSE abundances of W and Mo can be reproduced, that were previously hard to fit (Jennings et al., 2021). △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: Accepted in Icarus

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

An inflationary disk phase to explain extended protoplanetary dust disks

Authors: Raphael Marschall, Alessandro Morbidelli

Abstract: Understanding planetesimal formation is an essential first step to understanding planet formation. The distribution of these first solid bodies will drive the locations where planetary embryos can grow. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at l… ▽ More Understanding planetesimal formation is an essential first step to understanding planet formation. The distribution of these first solid bodies will drive the locations where planetary embryos can grow. We seek to understand the parameter space of possible protoplanetary disk formation and evolution models of our Solar System. A good protoplanetary disk scenario for the Solar System must meet at least the following three criteria: 1) an extended dust disk (at least 45 au); 2) formation of planetesimals in at least two distinct locations; and 3) transport of high temperatures condensates (i.e., calcium-aluminium-rich inclusion, CAIs) to the outer disk. We explore a large parameter space to study the effect of the disk viscosity, the timescale of infall of material into the disk, the distance within which material is deposited into the disk, and the fragmentation threshold of dust particles. We find that scenarios with a large initial disk viscosity ($α>0.05$), relatively short infall timescale ($T_{infall}<100-200$ kyr), and a small centrifugal radius ($R_C\sim0.4$~au; the distance within which material falls into the disk) result in disks that satisfy the criteria for a good protoplanetary disk of the Solar System. The large initial viscosity and short infall timescale result in a rapid initial expansion of the disk, which we dub the inflationary phase of the disk. Furthermore, a temperature-dependent fragmentation threshold, which mimics that cold icy particles break more easily, results in larger and more massive disks. This results in more "icy" than "rocky" planetesimals. Such scenarios are also better in line with our Solar System, which has small terrestrial planets and massive giant planet cores. Finally, we find that scenarios with large $R_C$ cannot transport CAIs to the outer disk and do not produce planetesimals at two locations within the disk. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: accepted to be published in A&A

Journal ref: A&A 677, A136 (2023)

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

Identification of a 4.3 billion year old asteroid family and planetesimal population in the Inner Main Belt

Authors: Salvatore Ferrone, Marco Delbo, Chrysa Avdellidou, Rogerio Deienno, Robert Melikyan, Kevin Walsh, Alessandro Morbidelli

Abstract: After performing a reassessment of the known dynamical asteroid families in the inner main belt, we report a newly discovered ancient asteroid family with an estimated age of $4.3\pm1.7$ billion years. Additionally, we report the most comprehensive list of planetesimals, which are asteroids that survived since the planet forming days of the solar system. After performing a reassessment of the known dynamical asteroid families in the inner main belt, we report a newly discovered ancient asteroid family with an estimated age of $4.3\pm1.7$ billion years. Additionally, we report the most comprehensive list of planetesimals, which are asteroids that survived since the planet forming days of the solar system. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 21 pages, 13 figures

Journal ref: A&A 676, A5 (2023)

In-situ enrichment in heavy elements of hot Jupiters

Authors: A. Morbidelli, K. Batygin, E. Lega

Abstract: Context: Radius and mass measurements of short-period giant planets reveal that many of these planets contain a large amount of heavy elements, in sharp contrast with the expectations of the conventional core-accretion model for the origin of giant planets. Aims: The proposed explanations for the heavy-element enrichment of giant planets fall short of explaining the most enriched planets. We loo… ▽ More Context: Radius and mass measurements of short-period giant planets reveal that many of these planets contain a large amount of heavy elements, in sharp contrast with the expectations of the conventional core-accretion model for the origin of giant planets. Aims: The proposed explanations for the heavy-element enrichment of giant planets fall short of explaining the most enriched planets. We look for additional processes that can explain the full envelope of inferred enrichments. Methods: We revisit the dynamics of pebbles and dust in the vicinity of giant planets using analytic estimates. Although our results are derived in the framework of a viscous alpha-disk we also discuss the case of disks driven by angular momentum removal in magnetized winds. Results: When giant planets are far from the star, dust and pebbles are confined in a pressure bump at the outer edge of the planet-induced gap. Instead, when the planets reach the inner part of the disk (r << 2 au), dust penetrates the gap together with the gas. The dust/gas ratio can be enhanced by more than an order of magnitude if radial drift of dust is not impeded farther out by other barriers. Thus, hot planets undergoing runaway gas accretion can swallow a large amount of dust. Conclusions: Whereas the gas accreted by giant planets in the outer disk is very dust-poor, that accreted by hot planets can be extremely dust-rich. Thus, provided that a large fraction of the atmosphere of hot-Jupiters is accreted in situ, a large amount of dust can be accreted as well. We draw a distinction between this process and pebble accretion, which is ineffective at small stellocentric radii, even for super-Earths. Giant planets farther out in the disk are extremely effective barriers against the flow of pebbles and dust across their gap. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: In press in Astronomy & Astrophysics

Journal ref: A&A 675, A75 (2023)

Cold Jupiters and improved masses in 38 Kepler and K2 small planet systems from 3661 HARPS-N radial velocities. No excess of cold Jupiters in small planet systems

Authors: A. S. Bonomo, X. Dumusque, A. Massa, A. Mortier, R. Bongiolatti, L. Malavolta, A. Sozzetti, L. A. Buchhave, M. Damasso, R. D. Haywood, A. Morbidelli, D. W. Latham, E. Molinari, F. Pepe, E. Poretti, S. Udry, L. Affer, W. Boschin, D. Charbonneau, R. Cosentino, M. Cretignier, A. Ghedina, E. Lega, M. López-Morales, M. Margini , et al. (9 additional authors not shown)

Abstract: The exoplanet population characterized by relatively short orbital periods ($P<100$ d) around solar-type stars is dominated by super-Earths and sub-Neptunes. However, these planets are missing in our Solar System and the reason behind this absence is still unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inw… ▽ More The exoplanet population characterized by relatively short orbital periods ($P<100$ d) around solar-type stars is dominated by super-Earths and sub-Neptunes. However, these planets are missing in our Solar System and the reason behind this absence is still unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inward migration of sub-Neptunes from beyond the water iceline; alternatively, Jupiter may have reduced considerably the inward flux of material (pebbles) required to form super-Earths inside that iceline. Both scenarios predict an anti-correlation between the presence of small planets (SPs) and that of cold Jupiters (CJs) in exoplanetary systems. To test that prediction, we homogeneously analyzed the radial-velocity (RV) measurements of 38 Kepler and K2 transiting SP systems gathered over almost 10 years with the HARPS-N spectrograph, as well as publicly available RVs collected with other facilities. We detected five CJs in three systems, two in Kepler-68, two in Kepler-454, and a very eccentric one in K2-312. We derived an occurrence rate of $9.3^{+7.7}_{-2.9}\%$ for CJs with $0.3-13~M_{Jup}$ and 1-10 AU, which is lower but still compatible at $1.3σ$ with that measured from RV surveys for solar-type stars, regardless of the presence or absence of SPs. The sample is not large enough to draw a firm conclusion about the predicted anti-correlation between SPs and CJs; nevertheless, we found no evidence of previous claims of an excess of CJs in SP systems. As an important by-product of our analyses, we homogeneously determined the masses of 64 Kepler and K2 small planets, reaching a precision better than 5, 7.5 and 10$σ$ for 25, 13 and 8 planets, respectively. Finally, we release the 3661 HARPS-N radial velocities used in this work to the scientific community. [Abridged] △ Less

Submitted 6 September, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 21 pages, 10 figures, 10 tables, published in Astronomy and Astrophysics. The updated version of the article takes into account the A&A language editing and guidelines. Tables 1, A.1 and full Table 2 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/677/A33

Journal ref: A&A 677, A33 (2023)

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

Terrestrial planet formation from a ring

Authors: J. M. Y. Woo, A. Morbidelli, S. L. Grimm, J. Stadel, R. Brasser

Abstract: It has been long proposed that, if all the terrestrial planets form within a tiny ring of solid material at around 1 AU, the concentrated mass-distance distribution of the current system can be reproduced. Recent planetesimal formation models also support this idea. In this study, we revisit the ring model by performing a number of high-resolution N-body simulations for 10 Myr of a ring of self-in… ▽ More It has been long proposed that, if all the terrestrial planets form within a tiny ring of solid material at around 1 AU, the concentrated mass-distance distribution of the current system can be reproduced. Recent planetesimal formation models also support this idea. In this study, we revisit the ring model by performing a number of high-resolution N-body simulations for 10 Myr of a ring of self-interacting planetesimals, with various radial distributions of the gas disc. We found that even if all the planetesimals form at ~1 AU in a minimum mass solar nebula-like disc, the system tends to spread radially as accretion proceeds, resulting in a system of planetary embryos lacking mass-concentration at ~1 AU. Modifying the surface density of the gas disc into a concave shape with a peak at ~1 AU helps to maintain mass concentrated at ~1 AU and solve the radial dispersion problem. We further propose that such a disc should be short lived (<= 1 Myr) and with a shallower radial gradient in the innermost region (< 1 AU) than previously proposed to prevent a too-rapid growth of Earth. Future studies should extend to ~100 Myr the most promising simulations and address in a self-consistent manner the evolution of the asteroid belt and its role in the formation of the terrestrial planets. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: Accepted by Icarus

Keck, Gemini, and Palomar 200-inch visible photometry of red and very-red Neptunian Trojans

Authors: B. T. Bolin, C. Fremling, A. Morbidelli, K. S. Noll, J. van Roestel, E. K. Deibert, M. Delbo, G. Gimeno, J. -E. Heo, C. M. Lisse, T. Seccull, H. Suh

Abstract: Neptunian Trojans (NTs), trans-Neptunian objects in 1:1 mean-motion resonance with Neptune, are generally thought to have been captured from the original trans-Neptunian protoplanetary disk into co-orbital resonance with the ice giant during its outward migration. It is possible, therefore, that the colour distribution of NTs is a constraint on the location of any colour transition zones that may… ▽ More Neptunian Trojans (NTs), trans-Neptunian objects in 1:1 mean-motion resonance with Neptune, are generally thought to have been captured from the original trans-Neptunian protoplanetary disk into co-orbital resonance with the ice giant during its outward migration. It is possible, therefore, that the colour distribution of NTs is a constraint on the location of any colour transition zones that may have been present in the disk. In support of this possible test, we obtained $g$, $r$, and $i$-band observations of 18 NTs, more than doubling the sample of NTs with known visible colours to 31 objects. Out of the combined sample, we found $\approx$4 objects with $g$-$i$ colours of $>$1.2 mags placing them in the very red (VR) category as typically defined. We find, without taking observational selection effects into account, that the NT $g$-$i$ colour distribution is statistically distinct from other trans-Neptunian dynamical classes. The optical colours of Jovian Trojans and NTs are shown to be less similar than previously claimed with additional VR NTs. The presence of VR objects among the NTs may suggest that the location of the red to VR colour transition zone in the protoplanetary disk was interior to 30-35 au. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: 8 pages, 3 figures, 3 tables, accepted for publication in MNRAS:Letters

Formation of Rocky Super-Earths From A Narrow Ring of Planetesimals

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belongin… ▽ More The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here, we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work - which demonstrates that planetesimals can form rapidly at discrete locations in the disk - we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 AU. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals, until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. Our results thus indicate that the absence of short-period super-Earths within the solar system can simply be attributed to the comparatively low mass of the primordial planetesimal ring within the protosolar nebula. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 14 pages, 4 figures, accepted for publication in Nature Astronomy

Gas accretion onto Jupiter mass planets in discs with laminar accretion flows

Authors: R. P. Nelson, E. Lega, A. Morbidelli

Abstract: (Abridged) Studies have shown that a Jovian mass planet embedded in a viscous protoplanetary disc (PPD) can accrete gas efficiently through the gap and doubles its mass in $\sim 0.1$ Myr. The planet also migrates inwards on a timescale of $\sim 0.1$ Myr. These timescales are short compared to PPD lifetimes, and raise questions about the origins of cold giant exoplanets. However, PPDs are unlikely… ▽ More (Abridged) Studies have shown that a Jovian mass planet embedded in a viscous protoplanetary disc (PPD) can accrete gas efficiently through the gap and doubles its mass in $\sim 0.1$ Myr. The planet also migrates inwards on a timescale of $\sim 0.1$ Myr. These timescales are short compared to PPD lifetimes, and raise questions about the origins of cold giant exoplanets. However, PPDs are unlikely to be globally turbulent, and instead they may launch magnetised winds such that accretion towards the star occurs in laminar accretion flows located in narrow layers near the surfaces of the disc. The aim of this study is to examine the rate at which gas accretes onto Jovian mass planets that are embedded in layered PPDs. We use 3D hydrodynamical simulations of planets embedded in PPDs, in which a constant radial mass flux towards the star of ${\dot m} = 10^{-8}$ M$_{\odot}$ yr$^{-1}$ is sustained. We consider a classical viscous alpha model, and also models in which an external torque is applied in narrow surface layers to mimic the effects of a magnetised wind. The accreting layers are parameterised by their column densities $Σ_{\rm A}$, and we consider values in the range 0.1 to 10 g cm$^{-2}$. The viscous model gives results in agreement with previous studies. We find the accretion rate onto the planet in the layered models crucially depends on the planet's ability to block the wind-induced mass flow. For $Σ_{\rm A}=10$ g cm$^{-2}$, the planet torque can block the mass flow through the disc, accretion onto the planet is slow, and a mass doubling time of 10 Myr is obtained. For $Σ_{\rm A}=0.1$ g cm$^{-2}$, accretion is fast and the mass doubling time is 0.2 Myr. Although the radial mass flow through the layered disc models is always $10^{-8}$ M$_{\odot}$ yr$^{-1}$, adopting different values of $Σ_{\rm A}$ leads to very different gas accretion rates onto gas giant planets. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 14 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 670, A113 (2023)

Rocky histories: The effect of high excitations on the formation of rocky planets

Authors: Jennifer Scora, Diana Valencia, Alessandro Morbidelli, Seth Jacobson

Abstract: Rocky planets both in and outside of our solar system are observed to have a range of core-mass fractions (CMFs). Imperfect collisions can preferentially strip mantle material from a planet, changing its CMF, and are therefore thought to be the most likely cause of this observed CMF variation. However, previous work that implements these collisions into N-body simulations of planet formation has s… ▽ More Rocky planets both in and outside of our solar system are observed to have a range of core-mass fractions (CMFs). Imperfect collisions can preferentially strip mantle material from a planet, changing its CMF, and are therefore thought to be the most likely cause of this observed CMF variation. However, previous work that implements these collisions into N-body simulations of planet formation has struggled to reliably form high CMF super-Earths. In this work, we specify our initial conditions and simulation parameters to maximize the prevalence of high-energy, CMF-changing collisions in order to form planets with highly diverse CMFs. High-energy collisions have a large $v_{imp}/v_{esc}$ ratio, so we maximize this ratio by starting simulations with high-eccentricity and inclination disks to increase the difference in their orbital velocities, maximizing $v_{imp}$. Additionally, we minimize $v_{esc}$ by starting with small embryos. The final planets undergo more high-energy, debris-producing collisions, and experience significant CMF change over their formation. However, we find that a number of processes work together to average out the CMF of a planet over time, therefore we do not consistently form high-CMF, high mass planets. We do form high-CMF planets below 0.5 $M_{\oplus}$. Additionally, we find in these highly eccentric environments, loss of debris mass due to collisional grinding has a significant effect on final planet masses and CMFs, resulting in smaller planets and a higher average planet CMF. This work highlights the importance of improving measurements of high-density planets to better constrain their CMFs. △ Less

Submitted 2 December, 2022; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 19 pages. This paper has been accepted to ApJ

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

Planetary Exploration Horizon 2061 Report, Chapter 3: From science questions to Solar System exploration

Authors: Véronique Dehant, Michel Blanc, Steve Mackwell, Krista M. Soderlund, Pierre Beck, Emma Bunce, Sébastien Charnoz, Bernard Foing, Valerio Filice, Leigh N. Fletcher, François Forget, Léa Griton, Heidi Hammel, Dennis Höning, Takeshi Imamura, Caitriona Jackman, Yohai Kaspi, Oleg Korablev, Jérémy Leconte, Emmanuel Lellouch, Bernard Marty, Nicolas Mangold, Patrick Michel, Alessandro Morbidelli, Olivier Mousis , et al. (9 additional authors not shown)

Abstract: This chapter of the Planetary Exploration Horizon 2061 Report reviews the way the six key questions about planetary systems, from their origins to the way they work and their habitability, identified in chapter 1, can be addressed by means of solar system exploration, and how one can find partial answers to these six questions by flying to the different provinces to the solar system: terrestrial p… ▽ More This chapter of the Planetary Exploration Horizon 2061 Report reviews the way the six key questions about planetary systems, from their origins to the way they work and their habitability, identified in chapter 1, can be addressed by means of solar system exploration, and how one can find partial answers to these six questions by flying to the different provinces to the solar system: terrestrial planets, giant planets, small bodies, and up to its interface with the local interstellar medium. It derives from this analysis a synthetic description of the most important space observations to be performed at the different solar system objects by future planetary exploration missions. These observation requirements illustrate the diversity of measurement techniques to be used as well as the diversity of destinations where these observations must be made. They constitute the base for the identification of the future planetary missions we need to fly by 2061, which are described in chapter 4. Q1- How well do we understand the diversity of planetary systems objects? Q2- How well do we understand the diversity of planetary system architectures? Q3- What are the origins and formation scenarios for planetary systems? Q4- How do planetary systems work? Q5- Do planetary systems host potential habitats? Q6- Where and how to search for life? △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 107 pages, 37 figures, Horizon 2061 is a science-driven, foresight exercise, for future scientific investigations

MSC Class: 86-02 ACM Class: A.1; H.1.0

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

A Criterion for the Stability of Planets in Chains of Resonances

Authors: Max Goldberg, Konstantin Batygin, Alessandro Morbidelli

Abstract: Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itse… ▽ More Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itself are poorly understood. Previous work postulated that a secondary resonance between the fastest libration frequency and a difference in synodic frequencies destabilizes the system. Here, we use that hypothesis to produce a simple and general criterion for resonant chain stability that depends only on planet orbital periods and masses. We show that the criterion accurately predicts the maximum mass of planets in synthetic resonant chains up to six planets. More complicated resonant chains produced in population synthesis simulations are found to be less stable than expected, although our criterion remains useful and superior to machine learning models. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 8 pages, 4 figures; accepted to Icarus

Stochastic accretion of the Earth

Authors: Paolo A. Sossi, Ingo L. Stotz, Seth A. Jacobson, Alessandro Morbidelli, Hugh St. C. O'Neill

Abstract: Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to… ▽ More Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to that expected from loss by partial vaporisation/condensation. Here we reconcile these observations by showing, using N-body simulations, that Earth accreted stochastically from many precursor bodies whose variable compositions reflect the temperatures at which they formed. Impact-induced atmospheric loss was efficient only when the proto-Earth was small, and elements that accreted thereafter retain near-chondritic isotope ratios. Earth's composition is reproduced when initial temperatures of planetesimal- to embryo-sized bodies are set by disk accretion rates of (1.08 $\pm$ 0.17) $\times$ 10$^{-7}$ solar masses/yr, although they may be perturbed by $^{26}$Al heating on bodies formed at different times. The model implies a heliocentric gradient in composition and rapid planetesimal formation within $\sim$ 1 Myr, in accord with radiometric volatile depletion ages of Earth. △ Less

Submitted 17 July, 2022; originally announced July 2022.

Comments: 13 pages, 4 figures. Nat Astron (2022)

A Self-Consistent Model for Dust-Gas Coupling in Protoplanetary Disks

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: Various physical processes that ensue within protoplanetary disks -- including vertical settling of icy/rocky grains, radial drift of solids, planetesimal formation, as well as planetary accretion itself -- are facilitated by hydrodynamic interactions between H/He gas and high-$Z$ dust. The Stokes number, which quantifies the strength of dust-gas coupling, thus plays a central role in protoplaneta… ▽ More Various physical processes that ensue within protoplanetary disks -- including vertical settling of icy/rocky grains, radial drift of solids, planetesimal formation, as well as planetary accretion itself -- are facilitated by hydrodynamic interactions between H/He gas and high-$Z$ dust. The Stokes number, which quantifies the strength of dust-gas coupling, thus plays a central role in protoplanetary disk evolution, and its poor determination constitutes an important source of uncertainty within the theory of planet formation. In this work, we present a simple model for dust-gas coupling, and demonstrate that for a specified combination of the nebular accretion rate, $\dot{M}$, and turbulence parameter, $α$, the radial profile of the Stokes number can be calculated uniquely. Our model indicates that the Stokes number grows sub-linearly with orbital radius, but increases dramatically across the water-ice line. For fiducial protoplanetary disk parameters of $\dot{M}=10^{-8}\,M_{\odot}/$year and $α=10^{-3}$, our theory yields characteristic values of the Stokes number on the order of $\mathrm{St}\sim10^{-4}$ (corresponding to $\sim$mm-sized silicate dust) in the inner nebula and $\mathrm{St}\sim10^{-1}$ (corresponding to $\sim$few-cm-sized icy grains), in the outer regions of the disk. Accordingly, solids are expected to settle into a thin sub-disk at large stellocentric distances, while remaining vertically well-mixed inside the ice line. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: accepted for publication in Astronomy & Astrophysics

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

Planet Formation Theory in the Era of ALMA and Kepler: from Pebbles to Exoplanets

Authors: Joanna Drazkowska, Bertram Bitsch, Michiel Lambrechts, Gijs D. Mulders, Daniel Harsono, Allona Vazan, Beibei Liu, Chris W. Ormel, Katherine Kretke, Alessandro Morbidelli

Abstract: Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constr… ▽ More Our understanding of planet formation has been rapidly evolving in recent years. The classical planet formation theory, developed when the only known planetary system was our own Solar System, has been revised to account for the observed diversity of the exoplanetary systems. At the same time, the increasing observational capabilities of the young stars and their surrounding disks bring new constraints on the planet formation process. In this chapter, we summarize the new information derived from the exoplanets population and the circumstellar disks observations. We present the new developments in planet formation theory, from dust evolution to the growth of planetary cores by accretion of planetesimals, pebbles, and gas. We review the state-of-the-art models for the formation of diverse planetary systems, including the population synthesis approach which is necessary to compare theoretical model outcomes to the exoplanet population. We emphasize that the planet formation process may not be spatially uniform in the disk and there are preferential locations for the formation of planetesimals and planets. Outside of these locations, a significant fraction of solids is not growing past the pebble-sizes. The reservoir of pebbles plays an important role in the growth of planetary cores in the pebble accretion process. The timescale of the emergence of massive planetary cores is an important aspect of the present models and it is likely that the cores within one disk form at different times. In addition, there is growing evidence that the first planetary cores start forming early, during the circumstellar disk buildup process. △ Less

Submitted 24 May, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 38 pages, 13 figures, review chapter accepted to Protostars and Planets VII, editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura, corrected label in Fig. 1

Dynamical origin of the Dwarf Planet Ceres

Authors: Rafael Ribeiro de Sousa, Alessandro Morbidelli, Rodney Gomes, Ernesto Vieira Neto, Andre Izidoro, Abreuçon Atanasio Alves

Abstract: The Dwarf Planet Ceres revealed the presence of ammonia and other unique properties compared to other asteroids in the main belt which suggests that it was not formed in situ. We model the early dynamical evolution of the outer Solar System to study possible dynamical mechanisms to implant a Ceres-sized planetesimal in the asteroid belt from the trans-Saturnian region. We calculate that the fracti… ▽ More The Dwarf Planet Ceres revealed the presence of ammonia and other unique properties compared to other asteroids in the main belt which suggests that it was not formed in situ. We model the early dynamical evolution of the outer Solar System to study possible dynamical mechanisms to implant a Ceres-sized planetesimal in the asteroid belt from the trans-Saturnian region. We calculate that the fraction of the population of Ceres-sized planetesimals that are captured in the asteroid belt is in the range of 2.8e-5 to 1.2e-3 depending on the initial location in the outer planetesimal disk. The captured bodies have a 70% probability to have a semimajor axis between 2.5 and 3 au, a 33% probability to have an eccentricity smaller than 0.2 and a 45% probability to have an orbital inclination smaller than 10 degrees. Assuming the existence of 3,600 Ceres-size planetesimals in the inner part of the trans-Saturnian disk, consistent with the estimate of Nesvorny & Vokrouhlicky (2016) for the trans-Neptunian disk, our estimated capture probability and a final 80% depletion of the asteroid belt during the subsequent giant planet instability, lead to capture 1 Ceres in the asteroid belt, with a probability of 15%, 34%, and 51% to be located in the inner, middle and outer belt respectively. △ Less

Submitted 18 February, 2022; originally announced February 2022.

The Metal-Silicate Partitioning of Carbon During Earth's Accretion and its Distribution in the Early Solar System

Authors: I. Blanchard, D. C. Rubie, E. S. Jennings, I. A. Franchi, X. Zhao, S. Petitgirard, N. Miyajima, S. A. Jacobson, A. Morbidelli

Abstract: Carbon is an essential element for the existence and evolution of life on Earth. Its abundance in Earth's crust and mantle (the Bulk Silicate Earth, BSE) is surprisingly high given that carbon is strongly siderophile (metal-loving) at low pressures and temperatures, and hence should have segregated almost completely into Earth's core during accretion. Estimates of the concentration of carbon in th… ▽ More Carbon is an essential element for the existence and evolution of life on Earth. Its abundance in Earth's crust and mantle (the Bulk Silicate Earth, BSE) is surprisingly high given that carbon is strongly siderophile (metal-loving) at low pressures and temperatures, and hence should have segregated almost completely into Earth's core during accretion. Estimates of the concentration of carbon in the BSE lie in the range 100-260 ppm and are much higher than expected based on simple models of core-mantle differentiation. Here we show through experiments at the putative conditions of Earth's core formation (49-71 GPa and 3600-4000 K) that carbon is significantly less siderophile at these conditions than at the low pressures ($\leq$13 GPa) and temperatures ($\leq$2500 K) of previous large volume press studies, but at least an order of magnitude more siderophile than proposed recently based on an experimental approach that is similar to ours. Using our new data along with previously published results, we derive a new parameterization of the pressure-temperature dependence of the metal-silicate partitioning of carbon. We apply this parameterization in a model that combines planet formation and core-mantle differentiation that is based on astrophysical N-body accretion simulations. Because differentiated planetesimals were almost completely depleted in carbon due to sublimation at high temperatures, almost all carbon in the BSE was added by the accretion of fully-oxidized carbonaceous chondrite material from the outer solar system. Carbon is added to the mantle continuously throughout accretion and its concentration reaches values within the BSE range (e.g. 140$\pm$40 ppm) at the end of accretion. The corresponding final core and bulk Earth carbon concentrations are 1270$\pm$300 ppm and 495$\pm$125 ppm respectively. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Supplementary material included in pdf

Journal ref: Earth and Planetary Science Letters. 580, 117374 (2022)

Terrestrial planet formation from lost inner solar system material

Authors: Christoph Burkhardt, Fridolin Spitzer, Alessandro Morbidelli, Gerrit Budde, Jan H. Render, Thomas S. Kruijer, Thorsten Kleine

Abstract: Two fundamentally different processes of rocky planet formation exist, but it is unclear which one built the terrestrial planets of the solar system. Either they formed by collisions among planetary embryos from the inner solar system, or by accreting sunward-drifting millimeter-sized 'pebbles' from the outer solar system. We show that the isotopic compositions of Earth and Mars are governed by tw… ▽ More Two fundamentally different processes of rocky planet formation exist, but it is unclear which one built the terrestrial planets of the solar system. Either they formed by collisions among planetary embryos from the inner solar system, or by accreting sunward-drifting millimeter-sized 'pebbles' from the outer solar system. We show that the isotopic compositions of Earth and Mars are governed by two-component mixing among inner solar system materials, including material from the innermost disk unsampled by meteorites, whereas the contribution of outer solar system material is limited to a few percent by mass. This refutes a pebble accretion origin of the terrestrial planets, but is consistent with collisional growth from inner solar system embryos. The low fraction of outer solar system material in Earth and Mars indicates the presence of a persistent dust-drift barrier in the disk, highlighting the specific pathway of rocky planet formation in the solar system. △ Less

Submitted 20 January, 2022; originally announced January 2022.

Comments: combined main test + supplement of published article

Journal ref: Science Advances 7, 52, abj7601 (2021)

Contemporary formation of early solar system planetesimals at two distinct radial locations

Authors: Alessandro Morbidelli, Kevin Baillie, Konstantin Batygin, Sebastien Charnoz, Tristan Guillot, David C. Rubie, Thorsten Kleine

Abstract: The formation of planetesimals is expected to occur via particle-gas instabilities that concentrate dust into self-gravitating clumps. Triggering these instabilities requires the prior pileup of dust in the protoplanetary disk. Until now, this has been successfully modeled exclusively at the disk's snowline, whereas rocky planetesimals in the inner disk were obtained only by assuming either unreal… ▽ More The formation of planetesimals is expected to occur via particle-gas instabilities that concentrate dust into self-gravitating clumps. Triggering these instabilities requires the prior pileup of dust in the protoplanetary disk. Until now, this has been successfully modeled exclusively at the disk's snowline, whereas rocky planetesimals in the inner disk were obtained only by assuming either unrealistically large particle sizes or an enhanced global disk metallicity. However, planetesimal formation solely at the snowline is difficult to reconcile with the early and contemporaneous formation of iron meteorite parent bodies with distinct oxidation states and isotopic compositions, indicating formation at different radial locations in the disk. Here, by modeling the evolution of a disk with ongoing accretion of material from the collapsing molecular cloud, we show that planetesimal formation may have been triggered within the first 0.5 million years by dust pileup at both the snowline (at approximately 5 au) and the silicate sublimation line (at approximately 1 au), provided turbulent diffusion was low. Particle concentration at approximately 1 au is due to the early outward radial motion of gas and is assisted by the sublimation and recondensation of silicates. Our results indicate that, although the planetesimals at the two locations formed about contemporaneously, those at the snowline accreted a large fraction of their mass (approximately 60 percent) from materials delivered to the disk in the first few 10,000 yr, whereas this fraction is only 30 percent for the planetesimals formed at the silicate line. Thus, provided that the isotopic composition of the delivered material changed with time, these two planetesimal populations should have distinct isotopic compositions, consistent with observations. △ Less

Submitted 31 December, 2021; originally announced December 2021.

Comments: Published online in Nature Astronomy on Dec. 22, 2021

Migration of Jupiter mass planets in discs with laminar accretion flows

Authors: E. Lega, A. Morbidelli, R. P. Nelson, X. S. Ramos, A. Crida, W. Bethune, K. Batygin

Abstract: Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal… ▽ More Migration of giant planets in discs with low viscosity has been studied recently. The proportionality between migration speed and the disc's viscosity is broken by the presence of vortices that appear at the edges of the planet-induced gap. Our goal is to investigate vortex-driven migration in low-viscosity discs in the presence of radial advection of gas, as expected from angular momentum removal by magnetised disc winds. We performed three dimensional simulations using the grid-based code FARGOCA. We mimicked the effects of a disc wind by applying a synthetic torque on a surface layer of the disc characterised by a prescribed column density Sigma_A so that it results in a disc accretion rate of 10^-8 Solar masses per year. Discs with this structure are called 'layered' and the layer where the torque is applied is denoted as 'active'. We also consider the case of accretion focussed near the disc midplane to mimic transport properties induced by a large Hall effect or by weak Ohmic diffusion. We observe two migration phases: in the first phase, the migration of the planet is driven by the vortex and is directed inwards. This phase ends when the vortex disappear. Migration depends on the ability of the torque from the planet to block the accretion flow. When the flow is fast and unimpeded (small Sigma_A) migration is very slow. When the accretion flow is completely blocked, migration is faster and the speed is controlled by the rate at which the accretion flow refills the gap behind the migrating planet. The migration speed of a giant planet in a layered protoplanetary disc depends on the thickness of the accreting layer. The lack of large-scale migration apparently experienced by the majority of giant exoplanets can be explained if the accreting layer is sufficiently thin to allow unimpeded accretion through the disc. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 14 pages, 11 figures

Journal ref: A&A 658, A32 (2022)

Planet Formation

Authors: Ravit Helled, Alessandro Morbidelli

Abstract: Our galaxy is full with planets. We now know that planets and planetary systems are diverse and come with different sizes, masses and compositions, as well as various orbital architectures. Although there has been great progress in understanding planet formation in the last couple of decades, both observationally and theoretically, several fundamental questions remain unsolved. This might not be s… ▽ More Our galaxy is full with planets. We now know that planets and planetary systems are diverse and come with different sizes, masses and compositions, as well as various orbital architectures. Although there has been great progress in understanding planet formation in the last couple of decades, both observationally and theoretically, several fundamental questions remain unsolved. This might not be surprising given the complexity of the process that includes various physical and chemical processes, and spans huge ranges of length-scales, masses, and timescales. In addition, planet formation cannot be directly observed but has to be inferred by gluing together different pieces of information into one consistent picture. "How do planets form?" remains a fundamental question in modern astrophysics. In this review we list some of the key open questions in planet formation theory as well as the challenges and upcoming opportunities. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: To appear as a book chapter in "ExoFrontiers: Big questions in exoplanetary science", Ed. N Madhusudhan (Bristol: IOP Publishing Ltd) AAS-IOP ebooks https://iopscience.iop.org/bookListInfo/aas-iop-astronomy

Constraints on planetesimal accretion inferred from particle-size distribution in CO chondrites

Authors: Gabriel A. Pinto, Yves Marrocchi, Alessandro Morbidelli, Sébastien Charnoz, Maria Eugenia Varela, Kevin Soto, Rodrigo Martínez, Felipe Olivares

Abstract: The formation of planetesimals was a key step in the assemblage of planetary bodies, yet many aspects of their formation remain poorly constrained. Notably, the mechanism by which chondrules -- sub-millimetric spheroids that dominate primitive meteorites -- were incorporated into planetesimals remains poorly understood. Here we classify and analyze particle-size distributions in various CO carbona… ▽ More The formation of planetesimals was a key step in the assemblage of planetary bodies, yet many aspects of their formation remain poorly constrained. Notably, the mechanism by which chondrules -- sub-millimetric spheroids that dominate primitive meteorites -- were incorporated into planetesimals remains poorly understood. Here we classify and analyze particle-size distributions in various CO carbonaceous chondrites found in the Atacama Desert. Our results show that the average circle-equivalent diameters of chondrules define a positive trend with the petrographic grade, which reflects the progressive role of thermal metamorphism within the CO parent body. We show that this relationship could not have been established by thermal metamorphism alone but rather by aerodynamic sorting during accretion. By modeling the self-gravitational contraction of clumps of chondrules, we show that (i) the accretion of the CO parent body(ies) would have generated a gradual change of chondrule size with depth in the parent body, with larger chondrules being more centrally concentrated than smaller ones, and (ii) any subsequent growth by pebble accretion would have been insignificant. These findings give substantial support to the view that planetesimals formed via gravitational collapse. △ Less

Submitted 19 August, 2021; v1 submitted 27 July, 2021; originally announced July 2021.

Comments: Accepted in ApJL. 24 pages, 3 figures, 2 tables, 1 appendix

Report number: CRPG contribution #2777

Migration of Jupiter mass planets in low viscosity discs

Authors: E. Lega, R. P. Nelson, A. Morbidelli, W. Kley, W. Béthune, A. Crida, D. Kloster, H. Méheut, T. Rometsch, A. Ziampras

Abstract: Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instabili… ▽ More Type-II migration of giant planets has a speed proportional to the disc's viscosity for values of the alpha viscosity parameter larger than 1.e-4 . At lower viscosities previous studies, based on 2D simulations have shown that migration can be very chaotic and often characterized by phases of fast migration. The reason is that in low-viscosity discs vortices appear due to the Rossby-wave instability at the edges of the gap opened by the planet. Migration is then determined by vortex-planet interactions. Our aim is to study migration in low viscosity 3D discs. We performed numerical simulations using 2D (including self-gravity) and 3D codes. After selecting disc masses for which self-gravity is not important, 3D simulations without self-gravity can be safely used. In our nominal simulation only numerical viscosity is present. We then performed simulations with prescribed viscosity to assess the threshold below which the new migration processes appear. We show that for alpha viscosity <= 1.e-5 two migration modes are possible which differ from classical Type-II migration, in the sense that they are not proportional to the disc's viscosity. The first occurs when the gap opened by the planet is not very deep. This occurs in 3D simulations and/or when a big vortex forms at the outer edge of the planetary gap, diffusing material into the gap. We call this type of migration "vortex-driven migration". This migration is very slow and cannot continue indefinitely, because eventually the vortex dissolves. The second migration mode occurs when the gap is deep so that the planet's eccentricity grows to a value ~0.2 due to inefficient eccentricity damping by corotation resonances. This second, faster migration mode appears to be typical of 2D models in discs with slower damping of temperature's perturbations. △ Less

Submitted 23 December, 2020; originally announced December 2020.

Comments: 18 pages, 13 figures

Journal ref: A&A 646, A166 (2021)

A re-assessment of the Kuiper belt size distribution for sub-kilometer objects, revealing collisional equilibrium at small sizes

Authors: A. Morbidelli, D. Nesvorny, W. F. Bottke, S. Marchi

Abstract: We combine several constraints provided by the crater records on Arrokoth and the worlds of the Pluto system to compute the size-frequency distribution (SFD) of the crater production function for craters with diameter D<10km. For this purpose, we use a Kuiper belt objects (KBO) population model calibrated on telescopic surveys, that describes also the evolution of the KBO population during the ear… ▽ More We combine several constraints provided by the crater records on Arrokoth and the worlds of the Pluto system to compute the size-frequency distribution (SFD) of the crater production function for craters with diameter D<10km. For this purpose, we use a Kuiper belt objects (KBO) population model calibrated on telescopic surveys, that describes also the evolution of the KBO population during the early Solar System. We further calibrate this model using the crater record on Pluto, Charon and Nix. Using this model, we compute the impact probability on Arrokoth, integrated over the age of the Solar System. This probability is then used together with other observational constraints to determine the slope of the crater-production function on Arrokoth. In addition, we use our Kuiper belt model also to compare the impact rates and velocities of KBOs on Arrokoth with those on Charon, integrated over the crater retention ages of their respective surfaces. This allows us to establish a relationship between the spatial density of sub-km craters on Arrokoth and of D~20km craters on Charon. Together, all these considerations suggest the crater production function on these worlds has a cumulative power law slope of -1.5<q<-1.2. Converted into a projectile SFD slope, we find -1.2<q_KBO<-1.0. These values are close to the cumulative slope of main belt asteroids in the 0.2--2~km range, a population in collisional equilibrium. For KBOs, however, this slope appears to extend from ~2km down to objects a few tens of meters in diameter, as inferred from sub-km craters on Arrokoth. From the measurement of the dust density in the Kuiper belt made by the New Horizons mission, we predict that the SFD of the KBOs becomes steep again below ~10-30m. All these considerations strongly indicate that the size distribution of the KBO population is in collisional equilibrium. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: In press in the special issue of Icarus "Pluto System, Kuiper Belt, and Kuiper Belt Objects"

Are the Moon's nearside-farside asymmetries the result of a giant impact?

Authors: Meng-Hua Zhu, Kai Wunnemann, Ross W. K. Potter, Thorsten Kleine, Alessandro Morbidelli

Abstract: The Moon exhibits striking geological asymmetries in elevation, crustal thickness, and composition between its nearside and farside. Although several scenarios have been proposed to explain these asymmetries, their origin remains debated. Recent remote sensing observations suggest that (1) the crust on the farside highlands consists of two layers: a primary anorthositic layer with thickness of ~30… ▽ More The Moon exhibits striking geological asymmetries in elevation, crustal thickness, and composition between its nearside and farside. Although several scenarios have been proposed to explain these asymmetries, their origin remains debated. Recent remote sensing observations suggest that (1) the crust on the farside highlands consists of two layers: a primary anorthositic layer with thickness of ~30-50 km and on top a more mafic-rich layer ~10 km thick; and (2) the nearside exhibits a large area of low-Ca pyroxene that has been interpreted to have an impact origin. These observations support the idea that the lunar nearside-farside asymmetries may be the result of a giant impact. Here, using quantitative numerical modeling, we test the hypothesis that a giant impact on the early Moon can explain the striking differences in elevation, crustal thickness, and composition between the nearside and farside of the Moon. We find that a large impactor, impacting the current nearside with a low velocity, can form a mega-basin and reproduce the characteristics of the crustal asymmetry and structures comparable to those observed on the current Moon, including the nearside lowlands and the farside's mafic-rich layer on top of a primordial anorthositic crust. Our model shows that the excavated deep-seated KREEP (potassium, rare-earth elements, and phosphorus) material, deposited close to the basin rim, slumps back into the basin and covers the entire basin floor; subsequent large impacts can transport the shallow KREEP material to the surface, resulting in its observed distribution. In addition, our model suggests that prior to the asymmetry-forming impact, the Moon may have had an 182W anomaly compared to the immediate post-giant impact Earth's mantle, as predicted if the Moon was created through a giant collision with the proto-Earth. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 45 pages, 12 figures

Journal ref: Journal of Geophysical Research: Planets, 2019

Reconstructing the late accretion history of the Moon

Authors: Meng-Hua Zhu, Natalia Artemieva, Alessandro Morbidelli, Qing-Zhu Yin, Harry Becker, Kai Wunnemann

Abstract: The importance of highly siderophile elements (HSEs) to track planetary late accretion has long been recognized. However, the precise nature of the Moon's accretional history remains enigmatic. There exists a significant mismatch of HSE budgets between the Earth and Moon, with the Earth disproportionally accreted far more HSEs than the Moon did. Several scenarios have been proposed to explain this… ▽ More The importance of highly siderophile elements (HSEs) to track planetary late accretion has long been recognized. However, the precise nature of the Moon's accretional history remains enigmatic. There exists a significant mismatch of HSE budgets between the Earth and Moon, with the Earth disproportionally accreted far more HSEs than the Moon did. Several scenarios have been proposed to explain this conundrum, including the delivery of HSEs to Earth by a few big impactors, the accretion of pebble-sized objects on dynamically cold orbits that enhanced the Earth's gravitational focusing factor, and the "sawtooth model" with much reduced impact flux before ~4.10 Gyr. However, most of these models assume a high impactor retention ratio f (fraction of impactor mass retained on the target) for the Moon. Here, we performed a series of impact simulations to quantify the f-value, followed by a Monte Carlo procedure enacting a monotonically decaying impact flux, to compute the mass accreted into lunar crust and mantle over their histories. We found that the average f-value for the Moon's entire impact history is about 3 times lower than previously estimated. Our results indicate that, to match the HSE budget of lunar crust and mantle, the retention of HSEs should have started ~ 4.35 Gyr ago, when most of lunar magma ocean was solidified. Mass accreted prior to 4.35 Gyr must have lost its HSE to the lunar core, presumably during the lunar mantle crystallization. The combination of a low impactor retention ratio and a late retention of HSEs in the lunar mantle provide a realistic explanation for the apparent deficit of Moon's late accreted mass relative to the Earth. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 3 figures

Journal ref: Nature,2019

Probing the impact of varied migration and gas accretion rates for the formation of giant planets in the pebble accretion scenario

Authors: Nelson Ndugu, Bertram Bitsch, Alessandro Morbidelli, Aurélien Crida, Edward Jurua

Abstract: The final orbital position of growing planets is determined by their migration speed, which is essentially set by the planetary mass. Small mass planets migrate in type I migration, while more massive planets migrate in type II migration, which is thought to depend mostly on the viscous evolution rate of the disc. A planet is most vulnerable to inward migration before it reaches type II migration… ▽ More The final orbital position of growing planets is determined by their migration speed, which is essentially set by the planetary mass. Small mass planets migrate in type I migration, while more massive planets migrate in type II migration, which is thought to depend mostly on the viscous evolution rate of the disc. A planet is most vulnerable to inward migration before it reaches type II migration and can lose a significant fraction of its semi-major axis at this stage. We investigated the influence of different disc viscosities, the dynamical torque and gas accretion from within the horseshoe region as mechanisms for slowing down planet migration. Our study confirms that planets growing in low viscosity environments migrate less, due to the earlier gap opening and slower type II migration rate. We find that taking the gas accretion from the horseshoe region into account allows an earlier gap opening and this results in less inward migration of growing planets. Furthermore, this effect increases the planetary mass compared to simulations that do not take the effect of gas accretion from the horseshoe region. Moreover, combining the effect of the dynamical torque with the effect of gas accretion from the horseshoe region, significantly slows down inward migration. Taking these effects into account could allow the formation of cold Jupiters (a > 1 au) closer to the water ice line region compared to previous simulations that did not take these effects into account. We thus conclude that gas accretion from within the horseshoe region and the dynamical torque play crucial roles in shaping planetary systems. △ Less

Submitted 11 December, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: 13 pages, 8 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

Effects of early intense bombardment on megaregolith evolution and on lunar (and planetary) surface samples

Authors: William K. Hartmann, Alessandro Morbidelli

Abstract: Impact rates in the first 500 Myr of the solar system are critical to an understanding of lunar geological history, but they have been controversial. The widely accepted, post-Apollo paradigm of early lunar impact cratering (ca. 1975-2014) proposed very low or negligible impact cratering in the period from accretion (>4.4 Ga) to about 4.0 Ga ago, followed by a 170-million-year-long spike of catacl… ▽ More Impact rates in the first 500 Myr of the solar system are critical to an understanding of lunar geological history, but they have been controversial. The widely accepted, post-Apollo paradigm of early lunar impact cratering (ca. 1975-2014) proposed very low or negligible impact cratering in the period from accretion (>4.4 Ga) to about 4.0 Ga ago, followed by a 170-million-year-long spike of cataclysmic cratering, during which most prominent multi-ring impact basins formed at age of about 3.9 Ga. More recent dynamical models suggest very early intense impact rates, declining throughout the period from accretion until an age of about 3.0 Ga. These models remove the basin-forming spike. This shift has important consequences on megaregolith evolution and properties of rock samples that can be collected on the lunar surface today. We adopt the Morbidelli et al. (2018) "accretion tail" model of early intense bombardment, declining as a function of time. We find effects differing from the previous models: early crater saturation and supersaturation; disturbance of magma ocean solidification; deep early megaregolith; and erosive destruction of the earliest multi-ring basins, their impact melts, and their ejecta blankets. Our results explain observations such as differences in numbers of early lunar impact melts vs. numbers of early igneous crustal rocks, highland breccias containing impact melts as old as 4.35 Ga, absence of a 170 Myr-long spike in impact melt ages at 3.9 Ga among lunar and asteroidal meteorites, and GRAIL observations of lunar crustal structure. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: in press in Meteoritics and Planetary Science (MAPS)

Evolution of the Earth's Atmosphere during Late Veneer Accretion

Authors: Catriona A. Sinclair, Mark C. Wyatt, Alessandro Morbidelli, David Nesvorny

Abstract: Recent advances in our understanding of the dynamical history of the Solar system have altered the inferred bombardment history of the Earth during accretion of the Late Veneer, after the Moon-forming impact. We investigate how the bombardment by planetesimals left-over from the terrestrial planet region after terrestrial planet formation, as well as asteroids and comets, affects the evolution of… ▽ More Recent advances in our understanding of the dynamical history of the Solar system have altered the inferred bombardment history of the Earth during accretion of the Late Veneer, after the Moon-forming impact. We investigate how the bombardment by planetesimals left-over from the terrestrial planet region after terrestrial planet formation, as well as asteroids and comets, affects the evolution of Earth's early atmosphere. We develop a new statistical code of stochastic bombardment for atmosphere evolution, combining prescriptions for atmosphere loss and volatile delivery derived from hydrodynamic simulations and theory with results from dynamical modelling of realistic populations of impactors. We find that for an initially Earth-like atmosphere impacts cause moderate atmospheric erosion with stochastic delivery of large asteroids giving substantial growth ($\times 10$) in a few $\%$ of cases. The exact change in atmosphere mass is inherently stochastic and dependent on the dynamics of the left-over planetesimals. We also consider the dependence on unknowns including the impactor volatile content, finding that the atmosphere is typically completely stripped by especially dry left-over planetesimals ($<0.02 ~ \%$ volatiles). Remarkably, for a wide range of initial atmosphere masses and compositions, the atmosphere converges towards similar final masses and compositions, i.e. initially low mass atmospheres grow whereas massive atmospheres deplete. While the final properties are sensitive to the assumed impactor properties, the resulting atmosphere mass is close to that of current Earth. The exception to this is that a large initial atmosphere cannot be eroded to the current mass unless the atmosphere was initially primordial in composition. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 30 pages, accepted for publication in MNRAS

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

Authors: Bryce T. Bolin, Christoffer Fremling, Timothy R. Holt, Matthew J. Hankins, Tomás Ahumada, Shreya Anand, Varun Bhalerao, Kevin B. Burdge, Chris M. Copperwheat, Michael Coughlin, Kunal P. Deshmukh, Kishalay De, Mansi M. Kasliwal, Alessandro Morbidelli, Josiah N. Purdum, Robert Quimby, Dennis Bodewits, Chan-Kao Chang, Wing-Huen Ip, Chen-Yen Hsu, Russ R. Laher, Zhong-Yi Lin, Carey M. Lisse, Frank J. Masci, Chow-Choong Ngeow , et al. (20 additional authors not shown)

Abstract: We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of… ▽ More We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 26 pages, 8 figures, 3 Tables, accepted for publication in ApJL

Isotopic evolution of the inner Solar System inferred from molybdenum isotopes in meteorites

Authors: Fridolin Spitzer, Christoph Burkhardt, Gerrit Budde, Thomas S. Kruijer, Alessandro Morbidelli, Thorsten Kleine

Abstract: The fundamentally different isotopic compositions of non-carbonaceous (NC) and carbonaceous (CC) meteorites reveal the presence of two distinct reservoirs in the solar protoplanetary disk that were likely separated by Jupiter. However, the extent of material exchange between these reservoirs, and how this affected the composition of the inner disk are not known. Here we show that NC meteorites dis… ▽ More The fundamentally different isotopic compositions of non-carbonaceous (NC) and carbonaceous (CC) meteorites reveal the presence of two distinct reservoirs in the solar protoplanetary disk that were likely separated by Jupiter. However, the extent of material exchange between these reservoirs, and how this affected the composition of the inner disk are not known. Here we show that NC meteorites display broadly correlated isotopic variations for Mo, Ti, Cr, and Ni, indicating the addition of isotopically distinct material to the inner disk. The added material resembles bulk CC meteorites and Ca-Al-rich inclusions in terms of its enrichment in neutron-rich isotopes, but unlike the latter materials is also enriched in s-process nuclides. The comparison of the isotopic composition of NC meteorites with the accretion ages of their parent bodies reveals that the isotopic variations within the inner disk do not reflect a continuous compositional change through the addition of CC dust, indicating an efficient separation of the NC and CC reservoirs and limited exchange of material between the inner and outer disk. Instead, the isotopic variations among NC meteorites more likely record a rapidly changing composition of the disk during infall from the Sun's parental molecular cloud, where each planetesimal locks the instant composition of the disk when it forms. A corollary of this model is that late-formed planetesimals in the inner disk predominantly accreted from secondary dust that was produced by collisions among pre-existing NC planetesimals. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: Accepted for publication in Astrophysical Journal Letters, June 17, 2020

No evidence for interstellar planetesimals trapped in the Solar System

Authors: A. Morbidelli, K. Batygin, R. Brasser, S. Raymond

Abstract: In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not… ▽ More In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not representative of the past evolution of real bodies. Instead, these simulations are only useful as a means to quantify the short dynamical lifetime of the considered bodies and the fast decay of their population. In light of this fast decay, if the observed bodies were the survivors of populations of objects captured from interstellar space in the early Solar System, these populations should have been implausibly large (e.g. about 10 times the current main asteroid belt population for the retrograde coorbital of Jupiter). More likely, the observed objects are just transient members of a population that is maintained in quasi-steady state by a continuous flux of objects from some parent reservoir in the distant Solar System. We identify in the Halley type comets and the Oort cloud the most likely sources of retrograde coorbitals and highly-inclined Centaurs. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: in press in MNRAS

Formation of Giant Planet Satellites

Authors: Konstantin Batygin, Alessandro Morbidelli

Abstract: Recent analyses have shown that the concluding stages of giant planet formation are accompanied by the development of large-scale meridional flow of gas inside the planetary Hill sphere. This circulation feeds a circumplanetary disk that viscously expels gaseous material back into the parent nebula, maintaining the system in a quasi-steady state. Here we investigate the formation of natural satell… ▽ More Recent analyses have shown that the concluding stages of giant planet formation are accompanied by the development of large-scale meridional flow of gas inside the planetary Hill sphere. This circulation feeds a circumplanetary disk that viscously expels gaseous material back into the parent nebula, maintaining the system in a quasi-steady state. Here we investigate the formation of natural satellites of Jupiter and Saturn within the framework of this newly outlined picture. We begin by considering the long-term evolution of solid material, and demonstrate that the circumplanetary disk can act as a global dust trap, where $s_{\bullet}\sim0.1-10\,$mm grains achieve a hydrodynamical equilibrium, facilitated by a balance between radial updraft and aerodynamic drag. This process leads to a gradual increase in the system's metallicity, and eventually culminates in the gravitational fragmentation of the outer regions of the solid sub-disk into $\mathcal{R}\sim100\,$km satellitesimals. Subsequently, satellite conglomeration ensues via pairwise collisions, but is terminated when disk-driven orbital migration removes the growing objects from the satellitesimal feeding zone. The resulting satellite formation cycle can repeat multiple times, until it is brought to an end by photo-evaporation of the parent nebula. Numerical simulations of the envisioned formation scenario yield satisfactory agreement between our model and the known properties of the Jovian and Saturnian moons. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Comments: 31 pages, 12 figures, accepted for publication in ApJ

The onset of instability in resonant chains

Authors: Gabriele Pichierri, Alessandro Morbidelli

Abstract: There is evidence that most chains of mean motion resonances of type $k$:$k-1$ among exoplanets become unstable once the dissipative action from the gas is removed from the system, particularly for large $N$ (the number of planets) and $k$ (indicating how compact the chain is). We present a novel dynamical mechanism that can explain the origin of these instabilities and thus the dearth of resonant… ▽ More There is evidence that most chains of mean motion resonances of type $k$:$k-1$ among exoplanets become unstable once the dissipative action from the gas is removed from the system, particularly for large $N$ (the number of planets) and $k$ (indicating how compact the chain is). We present a novel dynamical mechanism that can explain the origin of these instabilities and thus the dearth of resonant systems in the exoplanet sample. It relies on the emergence of secondary resonances between a fraction of the synodic frequency $2 π(1/P_1-1/P_2)$ and the libration frequencies in the mean motion resonance. These secondary resonances excite the amplitudes of libration of the mean motion resonances thus leading to an instability. We detail the emergence of these secondary resonances by carrying out an explicit perturbative scheme to second order in the planetary masses and isolating the harmonic terms that are associated with them. Focusing on the case of three planets in the 3:2 -- 3:2 mean motion resonance as an example, a simple but general analytical model of one of these resonances is obtained which describes the initial phase of the activation of one such secondary resonance. The dynamics of the excited system is also briefly described. This scheme shows how one can obtain analytical insight into the emergence of these resonances, and into the dynamics that they trigger. Finally, a generalisation of this dynamical mechanism is obtained for arbitrary $N$ and $k$. This leads to an explanation of previous numerical experiments on the stability of resonant chains, showing why the critical planetary mass allowed for stability decreases with increasing $N$ and $k$. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: Submitted to MNRAS

Planet formation by pebble accretion in ringed disks

Authors: Alessandro Morbidelli

Abstract: Context: Pebble accretion is expected to be the dominant process for the formation of massive solid planets, such as the cores of giant planets and super-Earths. So, far, this process has been studied under the assumption that dust coagulates and drifts throughout the full protoplanetary disk. However, observations show that many disks are structured in rings that may be due to pressure maxima, pr… ▽ More Context: Pebble accretion is expected to be the dominant process for the formation of massive solid planets, such as the cores of giant planets and super-Earths. So, far, this process has been studied under the assumption that dust coagulates and drifts throughout the full protoplanetary disk. However, observations show that many disks are structured in rings that may be due to pressure maxima, preventing the global radial drift of the dust. Aims: We study how the pebble-accretion paradigm changes if the dust is confined in a ring. Results: Planet Type-I migration is stopped in a ring, but not necessarily at its center. If the entropy-driven corotation torque is desaturated, the planet is located in a region with a low density of dust, which severely limits its accretion rate. If instead the planet is near the ring's center, its accretion rate can be similar to the one it would have in a classic (ring-less) disk of equivalent dust density. However, the growth rate of the planet is limited by the diffusion of dust in the ring and the final planet's mass is bounded by the total ring's mass. The DSHARP rings are too far from the star to allow the formation of massive planets within the disk's lifetime. However, a similar ring rescaled to 5 AU could lead to the formation of a planet incorporating the full ring's mass in less than 1/2 My. Conclusions: The existence of rings may not be an obstacle to planet formation by pebble-accretion. However, for accretion to be effective the resting position of the planet has to be relatively near the ring's center and the ring needs to be not too far from the central star. The formation of planets in rings can explain the existence of giant planets with core masses smaller than the so-called pebble isolation mass. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 638, A1 (2020)

Subsolar Al/Si and Mg/Si ratios of non-carbonaceous chondrites reveal planetesimal formation during early condensation in the protoplanetary disk

Authors: A. Morbidelli, G. Libourel, H. Palme, S. A. Jacobson, D. C. Rubie

Abstract: The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the s… ▽ More The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the subsequent formation of residual condensates with strongly sub-solar Al/Si and Mg/Si ratios. The mixing of residual condensates with different amounts of material with solar refractory/Si element ratios explains the Al/Si and Mg/Si values of non-carbonaceous chondrites. To match quantitatively the observed ratios, we find that the first-planetesimals should have accreted when the disk temperature was ~1,330-1,400 K depending on pressure and assuming a solar C/O ratio of the disk. We discuss how this model relates to our current understanding of disk evolution, grain dynamics, and planetesimal formation. We also extend the discussion to moderately volatile elements (e.g., Na), explaining how it may be possible that the depletion of these elements in non-carbonaceous chondrites is correlated with the depletion of refractory elements (e.g., Al). Extending the analysis to Cr, we find evidence for a higher than solar C/O ratio in the protosolar disk's gas when/where condensation from a fractionated gas occurred. Finally, we discuss the possibility that the supra-solar Al/Si and Mg/Si ratios of the Earth are due to the accretion of ~40% of the mass of our planet from the first-generation of refractory-rich planetesimals. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: Accepted for publication in EPSL