Search

Low-albedo asteroids preserve a record of the primordial solar system planetesimals and the conditions in which the solar nebula was active. However, the origin and evolution of these asteroids are not well-constrained. Here we measured visible and near-infrared (0.5 - 4.0 microns) spectra of low-albedo asteroids in the mid-outer main belt. We show that numerous large (d > 100 km) and dark (geometric albedo < 0.09) asteroids exterior to the dwarf planet Ceres' orbit share the same spectral features, and presumably compositions, as Ceres. We also developed a thermal evolution model that demonstrates that these Ceres-like asteroids have highly-porous interiors, accreted relatively late at 1.5 - 3.5 Myr after the formation of calcium-aluminum-rich inclusions, and experienced maximum interior temperatures of < 900 K. Ceres-like asteroids are localized in a confined heliocentric region between 3.0 - 3.4 au but were likely implanted from more distant regions of the solar system during the giant planet's dynamical instability.

Co-orbital systems contain two or more bodies sharing the same orbit around a planet or star. The best-known flavors of co-orbital systems are tadpoles (in which two bodies' angular separations oscillate about the L4/L5 Lagrange points $60^\circ$ apart) and horseshoes (with two bodies periodically exchanging orbital energy to trace out a horseshoe shape in a co-rotating frame). Here, we use N-body simulations to explore the parameter space of many-planet horseshoe systems. We show that up to 24 equal-mass, Earth-mass planets can share the same orbit at 1 au, following a complex pattern in which neighboring planets undergo horseshoe oscillations. We explore the dynamics of horseshoe constellations, and show that they can remain stable for billions of years and even persist through their stars' post-main sequence evolution. With sufficient observations, they can be identified through their large-amplitude, correlated transit timing variations. Given their longevity and exotic orbital architectures, horseshoe constellations may represent potential SETI beacons., 10 pages, 10 figures. Published in MNRAS. YouTube playlist with animations of horseshoe constellation systems here: https://www.youtube.com/playlist?list=PLelMZVM3ka3F335LGLxkxrD1ieiLJYQ5N . Blog post here: https://planetplanet.net/2023/04/20/constellations-of-co-orbital-planets/

One of the common approximations in long-term evolution studies of small bodies is the use of circular orbits averaging the actual eccentric ones, facilitating the coupling of processes with very different timescales, such as the orbital changes and the thermal processing. Here we test a number of averaging schemes for elliptic orbits in the context of the long-term evolution of comets, aiming to identify the one that best reproduces the elliptic orbits' heating patterns and the surface and subsurface temperature distributions. We use a simplified thermal evolution model applied on simulated comets both on elliptic and on their equivalent averaged circular orbits, in a range of orbital parameter space relevant to the inner solar system. We find that time averaging schemes are more adequate than spatial averaging ones. Circular orbits created by means of a time average of the equilibrium temperature approximate efficiently the subsurface temperature distributions of elliptic orbits in a large area of the orbital parameter space, rendering them a powerful tool for averaging elliptic orbits., 12 pages, 7 figures, accepted for publications on the AJ

Some systems of close-in "super-Earths" contain five or more planets on non-resonant but compact and nearly coplanar orbits. The Kepler-11 system is an iconic representative of this class of system. It is challenging to explain their origins given that planet-disk interactions are thought to be essential to maintain such a high degree of coplanarity, yet these same interactions invariably cause planets to migrate into chains of mean motion resonances. Here we mine a large dataset of dynamical simulations of super-Earth formation by migration. These simulations match the observed period ratio distribution as long as the vast majority of planet pairs in resonance become dynamically unstable. When instabilities take place resonances are broken during a late phase of giant impacts, and typical surviving systems have planet pairs with significant mutual orbital inclinations. However, a subset of our unstable simulations matches the Kepler-11 system in terms of coplanarity, compactness, planet-multiplicity and non-resonant state. This subset have dynamical instability phases typically much shorter than ordinary systems. Unstable systems may keep a high degree of coplanarity post-instability if planets collide at very low orbital inclinations ($\lesssim1^\circ$) or if collisions promote efficient damping of orbital inclinations. If planetary scattering during the instability takes place at low orbital inclinations ($\text{i}\lesssim1^\circ$), orbital inclinations are barely increased by encounters before planets collide.When planetary scattering pumps orbital inclinations to higher values ($\gtrsim 1^\circ$) planets tend to collide at higher mutual orbital inclinations, but depending on the geometry of collisions mergers' orbital inclinations may be efficiently damped. Each of these formation pathways can produce analogues to the Kepler-11 system., Comment: Accepted by MNRAS, 9 pages, 10 figures, 1 appendix

The inner solar system's modern orbital architecture provides inferences into the epoch of terrestrial planet formation; a ~100 Myr time period of planet growth via collisions with planetesimals and other proto-planets. While classic numerical simulations of this scenario adequately reproduced the correct number of terrestrial worlds, their semi-major axes and approximate formation timescales, they struggled to replicate the Earth-Mars and Venus-Mercury mass ratios. In a series of past independent investigations, we demonstrated that Mars' mass is possibly the result of Jupiter and Saturn's early orbital evolution, while Mercury's diminutive size might be the consequence of a primordial mass deficit in the region. Here, we combine these ideas in a single modeled scenario designed to simultaneously reproduce the formation of all four terrestrial planets and the modern orbits of the giant planets in broad strokes. By evaluating our Mercury analogs' core mass fractions, masses, and orbital offsets from Venus, we favor a scenario where Mercury forms through a series of violent erosive collisions between a number of ~Mercury-mass embryos in the inner part of the terrestrial disk. We also compare cases where the gas giants begin the simulation locked in a compact 3:2 resonant configuration to a more relaxed 2:1 orientation and find the former to be more successful. In 2:1 cases, the entire Mercury-forming region is often depleted due to strong sweeping secular resonances that also tend to overly excite the orbits of Earth and Venus as they grow. While our model is quite successful at replicating Mercury's massive core and dynamically isolated orbit, the planets' low mass remains extremely challenging to match. Finally, we discuss the merits and drawbacks of alternative evolutionary scenarios and initial disk conditions., Comment: 23 pages, 9 figures, 3 tables, accepted for publication in Icarus

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

Planets between 1-4 Earth radii with orbital periods, 13 pages, 8 figures and 1 table. Accepted for publication in MNRAS

Dynamical scenarios of terrestrial planets formation involve strong perturbations of the inner part of the solar system by the giant-planets, leading to enhanced impact velocities and subsequent collisional erosion. We quantitatively estimate the effect of collisional erosion on the resulting composition of Earth, and estimate how it may provide information on the dynamical context of its formation. We simulate and quantify the erosion of Earth's crust in the context of Solar System formation scenarios, including the classical model and Grand Tack scenario that invokes orbital migration of Jupiter during the gaseous disk phase (Walsh et al., 2011; Raymond et al., 2018). We find that collisional erosion of the early crust is unlikely to produce an excess of about 6% of the Sm/Nd ratio in terrestrial rock samples compared to chondrites for most simulations. Only Grand Tack simulations in which the last giant impact on Earth occurred later than 50 million years after the start of Solar System formation can account for such an offset. However, this time frame is consistent with current cosmochemical and dynamical estimates of the Moon forming impact (Chyba, 1991; Walker, 2009; Touboul et al.,2007, 2009, 2015; Pepin and Porcelli, 2006; Norman et al., 2003; Nyquist et al., 2006; Boyet et al.,2015). Such a late fractionation in the Sm/Nd ratio is unlikely to be responsible for a 20-ppm $^{142}$Nd excess in terrestrial rocks due to the half life of the radiogenic system. Additionally, such a large and late fractionation in the Sm/Nd ratio would accordingly induce non-observed anomalies in the $^{143}$Nd/$^{144}$Nd ratio. Considering our results, the Grand Tack model with a late Moon-forming impact cannot be easily reconciled with the Nd isotopic Earth contents., 22 pages, 3 figures, supplementary, to be published

The formation of the solar system's giant planets predated the ultimate epoch of massive impacts that concluded the process of terrestrial planet formation. Following their formation, the giant planets' orbits evolved through an episode of dynamical instability. Several qualities of the solar system have recently been interpreted as evidence of this event transpiring within the first ~100 Myr after the Sun's birth; around the same time as the final assembly of the inner planets. In a series of recent papers we argued that such an early instability could resolve several problems revealed in classic numerical studies of terrestrial planet formation; namely the small masses of Mars and the asteroid belt. In this paper, we revisit the early instability scenario with a large suite of simulations specifically designed to understand the degree to which Earth and Mars' formation are sensitive to the specific evolution of Jupiter and Saturn's orbits. By deriving our initial terrestrial disks directly from recent high-resolution simulations of planetesimal accretion, our results largely confirm our previous findings regarding the instability's efficiency of truncating the terrestrial disk outside of the Earth-forming region in simulations that best replicate the outer solar system. Moreover, our work validates the primordial 2:1 Jupiter-Saturn resonance within the early instability framework as a viable evolutionary path for the solar system. While our simulations elucidate the fragility of the terrestrial system during the epoch of giant planet migration, many realizations yield outstanding solar system analogs when scrutinized against a number of observational constraints. Finally, we highlight the inability of models to form adequate Mercury-analogs and the low eccentricities of Earth and Venus as the most significant outstanding problems for future numerical studies to resolve., 24 pages, 10 figures, 6 tables, accepted for publication in Icarus

Terrestrial planets in the Solar system generally lack surface liquid water. Earth is at odd with this observation and with the idea of the giant Moon-forming impact that should have vaporized any pre-existing water, leaving behind a dry Earth. Given the evidence available, this means that either water was brought back later or the giant impact could not vaporize all the water.We have looked at Venus for answers. Indeed, it is an example of an active planet that may have followed a radically different evolutionary pathway despite the similar mechanisms at work and probably comparable initial conditions. However, due to the lack of present-day plate tectonics, volatile recycling, and any surface liquid oceans, the evolution of Venus has likely been more straightforward than that of the Earth, making it easier to understand and model over its long term evolution.Here, we investigate the long-term evolution of Venus using self-consistent numerical models of global thermochemical mantle convection coupled with both an atmospheric evolution model and a late accretion N-body delivery model. We test implications of wet and dry late accretion compositions, using present-day Venus atmosphere measurements. Atmospheric losses are only able to remove a limited amount of water over the history of the planet. We show that late accretion of wet material exceeds this sink. CO2 and N2 contributions serve as additional constraints.Water-rich asteroids colliding with Venus and releasing their water as vapor cannot explain the composition of Venus atmosphere as we measure it today. It means that the asteroidal material that came to Venus, and thus to Earth, after the giant impact must have been dry (enstatite chondrites), therefore preventing the replenishment of the Earth in water. Because water can obviously be found on our planet today, it means that the water we are now enjoying on Earth has been there since its formation, likely buried deep in the Earth so it could survive the giant impact. This in turn suggests that suggests that planets likely formed with their near-full budget in water, and slowly lost it with time.

An episode of dynamical instability is thought to have sculpted the orbital structure of the outer solar system. When modeling this instability, a key constraint comes from Jupiter's fifth eccentric mode (quantified by its amplitude M55), which is an important driver of the solar system's secular evolution. Starting from commonly-assumed near-circular orbits, the present-day giant planets' architecture lies at the limit of numerically generated systems, and M55 is rarely excited to its true value. Here we perform a dynamical analysis of a large batch of artificially triggered instabilities, and test a variety of configurations for the giant planets' primordial orbits. In addition to more standard setups, and motivated by the results of modern hydrodynamical simulations of the giant planets' evolution within the primordial gaseous disk, we consider the possibility that Jupiter and Saturn emerged from the nebular gas locked in 2:1 resonance with non-zero eccentricities. We show that, in such a scenario, the modern Jupiter-Saturn system represents a typical simulation outcome, and M55 is commonly matched. Furthermore, we show that Uranus and Neptune's final orbits are determined by a combination of the mass in the primordial Kuiper belt and that of an ejected ice giant., 35 pages, 21 figures, 8 tables, accepted for publication in AJ

The standard model of terrestrial planet formation ignores the role of orbital migration of planetary embryos. A new scenario shows how migration may have sculpted the inner Solar System’s orbital architecture, as long as embryos converged towards about 1 au.

Planets intermediate in size between the Earth and Neptune, and orbiting closer to their host stars than Mercury does the Sun, are the most common type of planet revealed by exoplanet surveys over the last quarter century. Results from NASA's Kepler mission have revealed a bimodality in the radius distribution of these objects, with a relative underabundance of planets between 1.5 and 2.0 $R_{\oplus}$. This bimodality suggests that sub-Neptunes are mostly rocky planets that were born with primary atmospheres a few percent by mass accreted from the protoplanetary nebula. Planets above the radius gap were able to retain their atmospheres ("gas-rich super-Earths"), while planets below the radius gap lost their atmospheres and are stripped cores ("true super-Earths"). The mechanism that drives atmospheric loss for these planets remains an outstanding question, with photoevaporation and core-powered mass loss being the prime candidates. As with the mass-loss mechanism, there are two contenders for the origins of the solids in sub-Neptune planets: the migration model involves the growth and migration of embryos from beyond the ice line, while the drift model involves inward-drifting pebbles that coagulate to form planets close-in. Atmospheric studies have the potential to break degeneracies in interior structure models and place additional constraints on the origins of these planets. However, most atmospheric characterization efforts have been confounded by aerosols. Observations with upcoming facilities are expected to finally reveal the atmospheric compositions of these worlds, which are arguably the first fundamentally new type of planetary object identified from the study of exoplanets., Comment: Invited review in press in JGR: Planets special edition on exoplanets

In a recent paper we proposed that the giant planets' primordial orbits may have been eccentric (~0.05), and used a suite of dynamical simulations to show outcomes of the giant planet instability that are consistent with their present-day orbits. In this follow-up investigation, we present more comprehensive simulations incorporating superior particle resolution, longer integration times, and eliminating our prior means of artificially forcing instabilities to occur at specified times by shifting a planets' position in its orbit. While we find that the residual phase of planetary migration only minimally alters the the planets' ultimate eccentricities, our work uncovers several intriguing outcomes in realizations where Jupiter and Saturn are born with extremely large eccentricities (~0.10 and ~0.25, respectively). In successful simulations, the planets' orbits damp through interactions with the planetesimal disk prior to the instability, thus loosely replicating the initial conditions considered in our previous work. Our results therefore suggest an even wider range of plausible evolutionary pathways are capable of replicating Jupiter and Saturn's modern orbital architecture., Comment: 12 pages, 3 figures, 2 tables, accepted for publication in Icarus

Astronomical observations reveal that protoplanetary disks around young stars commonly have ring- and gap-like structures in their dust distributions. These features are associated with pressure bumps trapping dust particles at specific locations, which simulations show are ideal sites for planetesimal formation. Here we show that our Solar System may have formed from rings of planetesimals -- created by pressure bumps -- rather than a continuous disk. We model the gaseous disk phase assuming the existence of pressure bumps near the silicate sublimation line (at $T \sim$1400~K), water snowline (at $T \sim$170~K), and CO-snowline (at $T \sim$30~K). Our simulations show that dust piles up at the bumps and forms up to three rings of planetesimals: a narrow ring near 1~au, a wide ring between $\sim$3-4~au and $\sim$10-20~au, and a distant ring between $\sim$20~au and $\sim$45~au. We use a series of simulations to follow the evolution of the innermost ring and show how it can explain the orbital structure of the inner Solar System and provides a framework to explain the origins of isotopic signatures of Earth, Mars and different classes of meteorites. The central ring contains enough mass to explain the rapid growth of the giant planets' cores. The outermost ring is consistent with dynamical models of Solar System evolution proposing that the early Solar System had a primordial planetesimal disk beyond the current orbit of Uranus., Comment: Published in Nature Astronomy on Dec 30, 2021. Authors' version including Methods and Supplementary Information

In the pebble accretion scenario, the pebbles that form planets drift inward from the outer disk regions, carrying water ice with them. At the water ice line, the water ice on the inward drifting pebbles evaporates and is released into the gas phase, resulting in water-rich gas and dry pebbles that move into the inner disk regions. Large planetary cores can block the inward drifting pebbles by forming a pressure bump outside their orbit in the protoplanetary disk. Depending on the relative position of a growing planetary core relative to the water ice line, water-rich pebbles might be blocked outside or inside the water ice line. Pebbles blocked outside the water ice line do not evaporate and thus do not release their water vapor into the gas phase, resulting in a dry inner disk, while pebbles blocked inside the water ice line release their water vapor into the gas phase, resulting in water vapor diffusing into the inner disk. As a consequence, close-in sub-Neptunes that accrete some gas from the disk should be dry or wet, respectively, if outer gas giants are outside or inside the water ice line, assuming that giant planets form fast, as has been suggested for Jupiter in our Solar System. Alternatively, a sub-Neptune could form outside the water ice line, accreting a large amount of icy pebbles and then migrating inward as a very wet sub-Neptune. We suggest that the water content of inner sub-Neptunes in systems with giant planets that can efficiently block the inward drifting pebbles could constrain the formation conditions of these systems, thus making these sub-Neptunes exciting targets for detailed characterization (e.g., with JWST, ELT, or ARIEL). In addition, the search for giant planets in systems with already characterized sub-Neptunes can be used to constrain the formation conditions of giant planets as well., Comment: Accepted by A&A Letters

The Solar System’s orbital structure is thought to have been sculpted by a dynamical instability among the giant planets[1–4]. Yet the instability trigger and exact timing have proved hard to pin down[5–9]. The giant planets formed within a gas-dominated disk around the young Sun. Motivated by giant exoplanet systems found in mean motion resonance[10], hydrodynamical modeling has shown that while the disk was present the giant planets migrated into a compact orbital configuration, in a chain of resonances[2,11]. Here we use a suite of dynamical simulations to show that the giant planets’ instability was likely triggered by the dispersal of the Sun’s gaseous disk. As the disk evaporated from the inside-out, its inner edge swept successively across and dynamically perturbed each planet’s orbit in turn. Saturn and each ice giants’ orbits were torqued strongly enough to migrate outward. As a given planet migrated outward with the disk’s inner edge the orbital configuration of the exterior system was compressed, triggering dynamical instability. The final orbits of our simulated systems match those of the Solar System for a viable range of astrophysical parameters. Our results demonstrate that the giant planet instability happened as the gaseous disk dissipated, constrained by astronomical observations to be a few to ten million years after the birth of the Solar System [12]. Late-stage terrestrial planet formation would occur mostly after such an early giant planet instability [13,14], thereby avoiding the possibility of de-stabilizing the terrestrial planets [15] and naturally accounting for the small mass of Mars relative to Earth and the mass depletion of the main asteroid belt [16].

The Solar System's orbital structure is thought to have been sculpted by an episode of dynamical instability among the giant planets

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system., Accepted for publication in AJ (30 pages, 13 figures, 4 tables, 2 appendices)

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., Comment: in press in MNRAS

With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c)., Comment: 50 pages, 21 figures. Accepted for publication in Astronomy and Astrophysics

The dynamical structure of the Solar System can be explained by a period of orbital instability experienced by the giant planets. While a late instability was originally proposed to explain the Late Heavy Bombardment, recent work favors an early instability. We model the early dynamical evolution of the outer Solar System to self-consistently constrain the most likely timing of the instability. We first simulate the dynamical sculpting of the primordial outer planetesimal disk during the accretion of Uranus and Neptune from migrating planetary embryos during the gas disk phase, and determine the separation between Neptune and the inner edge of the planetesimal disk. We performed simulations with a range of migration histories for Jupiter. We find that, unless Jupiter migrated inwards by 10 AU or more, the instability almost certainly happened within 100 Myr of the start of Solar System formation. There are two distinct possible instability triggers. The first is an instability that is triggered by the planets themselves, with no appreciable influence from the planetesimal disk. Of those, the median instability time is $\sim4$Myr. Among self-stable systems -- where the planets are locked in a resonant chain that remains stable in the absence of a planetesimal's disk-- our self-consistently sculpted planetesimal disks nonetheless trigger a giant planet instability with a median instability time of 37-62 Myr for a reasonable range of migration histories of Jupiter. The simulations that give the latest instability times are those that invoked long-range inward migration of Jupiter from 15 AU or beyond; however these simulations over-excited the inclinations of Kuiper belt objects and are inconsistent with the present-day Solar System. We conclude on dynamical grounds that the giant planet instability is likely to have occurred early in Solar System history., 46 pages, 26 figures, Article reference YICAR_113605, https://authors.elsevier.com/tracking/article/details.do?aid=113605&jid=YICAR&surname=Ribeiro

It is widely known that giant planets have the capacity to open deep gaps in their natal gaseous protoplanetary discs. It is unclear, however, how gas accretion onto growing planets influences the shape and depth of their growing gaps. We performed isothermal hydrodynamical simulations with the Fargo-2D1D code, which assumes planets accreting gas within full discs that range from 0.1 to 260 AU. The gas accretion routine uses a sink cell approach, in which different accretion rates are used to cope with the broad range of gas accretion rates cited in the literature. We find that the planetary gas accretion rate increases for larger disc aspect ratios and greater viscosities. Our main results show that gas accretion has an important impact on the gap-opening mass: we find that when the disc responds slowly to a change in planetary mass (i.e., at low viscosity), the gap-opening mass scales with the planetary accretion rate, with a higher gas accretion rate resulting in a larger gap-opening mass. On the other hand, if the disc response time is short (i.e., at high viscosity), then gas accretion helps the planet carve a deep gap. As a consequence, higher planetary gas accretion rates result in smaller gap-opening masses. Our results have important implications for the derivation of planet masses from disc observations: depending on the planetary gas accretion rate, the derived masses from ALMA observations might be off by up to a factor of two. We discuss the consequences of the change in the gap-opening mass on the evolution of planetary systems based on the example of the grand tack scenario. Planetary gas accretion also impacts stellar gas accretion, where the influence is minimal due to the presence of a gas-accreting planet.

The orbital architecture of the Solar System is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space but a fraction were trapped on stable orbits in the Kuiper belt and Oort cloud. We use a suite of N-body simulations to map out the diversity of planetesimals' dynamical pathways. We focus on two processes: tidal disruption from very close encounters with a giant planet, and loss of surface volatiles from repeated passages close to the Sun. We show that the rate of tidal disruption is more than a factor of two higher for ejected planetesimals than for surviving objects in the Kuiper belt or Oort cloud. Ejected planetesimals are preferentially disrupted by Jupiter and surviving ones by Neptune. Given that the gas giants contracted significantly as they cooled but the ice giants did not, taking into account the thermal evolution of the giant planets decreases the disruption rate of ejected planetesimals. The frequency of volatile loss and extinction is far higher for ejected planetesimals than for surviving ones and is not affected by the giant planets' contraction. Even if all interstellar objects were ejected from Solar System-like systems, our analysis suggests that their physical properties should be more diverse than those of Solar System small bodies as a result of their divergent dynamical histories. This is consistent with the characteristics of the two currently-known interstellar objects., ApJ Letters, in press. 6 pages, 4 figures

It remains contentious whether the meteoritic material delivered to the terrestrial planets after the end of core formation was rich or poor in water and other volatiles. As Venus’s atmosphere has probably experienced less volatile recycling over its history than Earth’s, it may be possible to constrain the volatile delivery to the primitive Venusian atmosphere from the planet’s present-day atmospheric composition. Here we investigate the long-term evolution of Venus using self-consistent numerical simulations of global thermochemical mantle convection coupled with both an atmospheric evolution model and a late accretion N-body delivery model. We found that atmospheric escape is only able to remove a limited amount of water over the history of the planet, and that the late accretion of wet material exceeds this sink and would result in a present-day atmosphere that is too rich in volatiles. A preferentially dry composition of the late accretion impactors is most consistent with measurements of atmospheric H2O, CO2 and N2. Hence, we suggest that the late accreted material delivered to Venus was mostly dry enstatite chondrite, consistent with isotopic data for Earth, with less than 2.5% (by mass) wet carbonaceous chondrites. In this scenario, the majority of Venus’s and Earth’s water would have been delivered during the main accretion phase., SCOPUS: ar.j, info:eu-repo/semantics/published

The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice Model explains many peculiar qualities of the solar system, including the belt's excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice Model's effect on the belt's structure struggle to reproduce the innermost asteroids' orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets' mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn's, Jupiter's eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24-28 arcsec/yr is related to the inclination problem. As Saturn's precession speeds back up, high inclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt's orbital structure is reshaped, leading to markedly improved simulation outcomes., 5 pages, 4 figures, 2 tables, accepted for publication in MRNAS Letters

The discovery of the first interstellar object passing through the Solar System, 1I/2017 U1 ('Oumuamua), provoked intense and continuing interest from the scientific community and the general public. The faintness of 'Oumuamua, together with the limited time window within which observations were possible, constrained the information available on its dynamics and physical state. Here we review our knowledge and find that in all cases the observations are consistent with a purely natural origin for 'Oumuamua. We discuss how the observed characteristics of 'Oumuamua are explained by our extensive knowledge of natural minor bodies in our Solar System and our current knowledge of the evolution of planetary systems. We highlight several areas requiring further investigation., Comment: Published in Nature Astronomy, 01 July 2019. 16 pages, 4 figures, 1 table

In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number of processes, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone at $R\lesssim 10$ AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical simulations to study the interaction between a low-mass planet ($\sim 1.7 \;{\rm M_\oplus}$) on a fixed orbit and an inviscid pebble-rich disc with solid-to-gas ratio $\epsilon\ge 0.5$. For pebbles with Stokes numbers St=0.1, 0.5, multiple dusty vortices are formed through the Rossby Wave Instability at the planet separatrix. Effects due to gas drag then lead to a strong enhancement in the solid-to-gas ratio, which can increase by a factor of $\sim 10^3$ for marginally coupled particles with St=0.5. As in streaming instabilities, pebble clumps reorganize into filaments that may plausibly collapse to form planetesimals. When the planet is allowed to migrate in a MMSN disc, the vortex instability is delayed due to migration but sets in once inward migration stops due a strong positive pebble torque. Again, particle filaments evolving in a gap are formed in the disc while the planet undergoes an episode of outward migration. Our results suggest that vortex instabilities triggered by low-mass planets could play an important role in forming planetesimals in pebble-rich, inviscid discs, and may significantly modify the migration of low-mass planets. They also imply that planetary dust gaps may not necessarily contain planets if these migrated away., Comment: Accepted in MNRAS

Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r, 25 pages, accepted by A&A, see companion papers by Izidoro et al. (2019) and Lambrechts et al. (2019)