Protoplanets Research Articles

Could Uranus and Neptune form by collisions of planetary embryos?

ABSTRACT The origin of Uranus and Neptune remains a challenge for planet formation models. A potential explanation is that the planets formed from a population of a few planetary embryos with masses of a few Earth masses which formed beyond Saturn’s orbit and migrated inwards. These embryos can collide and merge to form Uranus and Neptune. In this work, we revisit this formation scenario and study the outcomes of such collisions using 3D hydrodynamical simulations. We investigate under what conditions the perfect-merging assumption is appropriate, and infer the planets’ final masses, obliquities, and rotation periods, as well as the presence of proto-satellite discs. We find that the total bound mass and obliquities of the planets formed in our simulations generally agree with N-body simulations therefore validating the perfect-merging assumption. The inferred obliquities, however, are typically different from those of Uranus and Neptune, and can be roughly matched only in a few cases. In addition, we find that in most cases, the planets formed in this scenario rotate faster than Uranus and Neptune, close to break-up speed, and have massive discs. We therefore conclude that forming Uranus and Neptune in this scenario is challenging, and further research is required. We suggest that future planet formation models should aim to explain the various physical properties of the planets such as their masses, compositions, obliquities, rotation rates, and satellite systems.

Linking planetary embryo formation to planetesimal formation

Context. The growth-timescales of planetary embryos and their formation process are imperative for our understanding on how planetary systems form and develop. They determine the subsequent growth mechanisms during the life stages of a circumstellar disk. Aims. We quantify the timescales and spatial distribution of planetary embryos through collisional growth and fragmentation of dynamically forming 100 km sized planetesimals. In our study, the formation timescales of viscous disk evolution and planetesimal formation are linked to the formation of planetary embryos in the terrestrial planet zone. Methods. We connected a one-dimensional model for viscous gas evolution, dust and pebble dynamics, and pebble flux-regulated planetesimal formation to the N-body code LIPAD. Our framework enabled us to study the formation, growth, fragmentation, and evolution of planetesimals with an initial size of 100 km in diameter for the first million years of a viscous disk. Results. Our study shows the effect of the planetesimal surface density evolution on the preferential location and timescales of planetary embryo formation. Only the innermost embryos (<2 au) in our study form well within the lifetime of an active pebble flux for any disk studied. Higher planetesimal disk masses and steeper planetesimal surface density profiles result in more massive embryos within a larger area, rather than in a higher number of embryos. A one-dimensional analytically derived model for embryo formation based on the local planetesimal surface density evolution is presented. This model can reproduce the spatial distribution, formation rate, and total number of planetary embryos at a fraction of the computational cost of the N-body simulations. Conclusions. The formation of planetary embryos in the terrestrial planet zone occurs simultaneously with the formation of planetesimals. The local planetesimal surface density evolution and the orbital spacing of planetary embryos in the oligarchic regime are good constraints for modeling planetary embryo formation analytically. Our embryo formation model is a valuable asset in future studies of planet formation.

Oligarchic growth in a fully interacting system

Context. In the oligarchic growth model, protoplanets develop in the final stage of planet formation via collisions between planetesimals and planetary embryos. The majority of planetesimals are accreted by the embryos, while the remnant planetesimals acquire dynamically excited orbits. The efficiency of planet formation can be defined by the mass ratio between formed protoplanets and the initial mass of the embryo-planetesimal belt. Aims. In numerical simulations of the oligarchic growth, the gravitational interactions between planetesimals are usually neglected due to computational difficulties. In this way, computations require fewer resources. We investigated the effect of this simplification by modeling the planet formation efficiency in a belt of embryos with self-interacting or non-self-interacting planetesimals. Methods. We used our own graphics processing unit-based direct N-body integrator for the simulations. We compared 2D models using different initial embryo numbers, different initial planetesimal numbers, and different total initial belt masses. For a limited number of cases, we compared the 2D and 3D simulations. Results. We found that planet formation efficiency is higher if the planetesimal self-interaction is taken into account in models that contain the commonly used 100 embryos. The observed effect can be explained by the damping of planetesimal eccentricities by their self-gravity. The final numbers of protoplanets are independent of planetesimal self-gravity, while the average mass of the formed protoplanets is larger in the self-interacting models. We also found that the non-self-interacting and self-interacting models qualitatively give the same results above 200 embryos. Our findings show that the higher the initial mass of the embryo-planetesimal belt, the higher the discrepancy between models that use self-interacting or non-self-interacting planetesimals is. The study of 3D models showed quantitatively the same results as the 2D models for low average inclination. We conclude that it is important to include planetesimal self-interaction in both 2D and 3D models in cases where the initial embryo number is less than 200.

Linking planetary embryo formation to planetesimal formation

Context.The accretion of pebbles onto planetary cores has been widely studied in recent years and is found to be a highly effective mechanism for planetary growth. While most studies assume planetary cores as an initial condition in their simulation, the question of the manner, location, and time at which these cores form is often neglected.Aims.We study the effect of pebble accretion during the formation phase and subsequent evolution of planetary embryos in the early stages of circumstellar disk evolution. In doing so, we aim to quantify the timescales and local dependence of planetary embryo formation based on the solid evolution of the disk.Methods.We connected a one-dimensional two-population model for solid evolution and pebble-flux-regulated planetesimal formation to theN-body code LIPAD. We focused on the growth of planetesimals with an initial size of 100 km in diameter by planetesimal collisions and pebble accretion for the first one million years of a viscously evolving disk. We compared 18 differentN-body simulations in which we varied the total planetesimal mass after one million years, the surface density profile of the planetesimal disk, the radial pebble flux, and the possibility of pebble accretion.Results.Pebble accretion leads to the formation of fewer but substantially more massive embryos. The area of possible embryo formation is weakly affected by the accretion of pebbles, and the innermost embryos tend to form slightly earlier than in simulations in which pebble accretion is neglected.Conclusions.Pebble accretion strongly enhances the formation of super-Earths in the terrestrial planet region, but it does not enhance the formation of embryos at larger distances.

Formation of Venus, Earth and Mars: Constrained by Isotopes

We discuss the current state of knowledge of terrestrial planet formation from the aspects of different planet formation models and isotopic data from 182Hf-182W, U-Pb, lithophile-siderophile elements, 48Ca/44Ca isotope samples from planetary building blocks, 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne isotope ratios in Venus' and Earth's atmospheres, the expected solar 3He abundance in Earth's deep mantle and Earth's D/H sea water ratios that shed light on the accretion time of the early protoplanets. Accretion scenarios that can explain the different isotope ratios, including a Moon-forming event after ca. 50 Myr, support the theory that the bulk of Earth's mass (>80%) most likely accreted within 10-30 Myr. From a combined analysis of the before mentioned isotopes, one finds that proto-Earth accreted 0.5-0.6 MEarth within the first ~4-5 Myr, the approximate lifetime of the protoplanetary disk. For Venus, the available atmospheric noble gas data are too uncertain for constraining the planet's accretion scenario accurately. However, from the available Ar and Ne isotope measurements, one finds that proto-Venus could have grown to 0.85-1.0 MVenus before the disk dissipated. Classical terrestrial planet formation models have struggled to grow large planetary embryos quickly from the tiniest materials within the typical lifetime of protoplanetary disks. Pebble accretion could solve this long-standing time scale controversy. Pebble accretion and streaming instabilities produce large planetesimals that grow into Mars-sized and larger planetary embryos during this early accretion phase. The later stage of accretion can be explained well with the Grand-Tack, annulus or depleted disk models. The relative roles of pebble accretion and planetesimal accretion/giant impacts are poorly understood and should be investigated with N-body simulations that include pebbles and multiple protoplanets.

Losing oceans: The effects of composition on the thermal component of impact-driven atmospheric loss

ABSTRACT The formation of the Solar system’s terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20–30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be efficiently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. We demonstrate that a giant impact during terrestrial planet formation can remove several Earth oceans’ worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an explanation for the observed depletion in Earth’s light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets.

Grain growth for astrophysics with discontinuous Galerkin schemes

ABSTRACTDepending on their sizes, dust grains store more or less charges, catalyse more or less chemical reactions, intercept more or less photons and stick more or less efficiently to form embryos of planets. Hence, the need for an accurate treatment of dust coagulation and fragmentation in numerical modelling. However, existing algorithms for solving the coagulation equation are overdiffusive in the conditions of 3D simulations. We address this challenge by developing a high-order solver based on the discontinuous Galerkin method. This algorithm conserves mass to machine precision and allows to compute accurately the growth of dust grains over several orders of magnitude in size with a very limited number of dust bins.

The eccentricity distribution of giant planets and their relation to super-Earths in the pebble accretion scenario

Observations of the population of cold Jupiter planets (r >1 AU) show that nearly all of these planets orbit their host star on eccentric orbits. For planets up to a few Jupiter masses, eccentric orbits are thought to be the outcome of planet–planet scattering events taking place after gas dispersal. We simulated the growth of planets via pebble and gas accretion as well as the migration of multiple planetary embryos in their gas disc. We then followed the long-term dynamical evolution of our formed planetary system up to 100 Myr after gas disc dispersal. We investigated the importance of the initial number of protoplanetary embryos and different damping rates of eccentricity and inclination during the gas phase for the final configuration of our planetary systems. We constrained our model by comparing the final dynamical structure of our simulated planetary systems to that of observed exoplanet systems. Our results show that the initial number of planetary embryos has only a minor impact on the final orbital eccentricity distribution of the giant planets, as long as the damping of eccentricity and inclination is efficient. If the damping is inefficient (slow), systems with a larger initial number of embryos harbour larger average eccentricities. In addition, for slow damping rates, we observe that scattering events are already common during the gas disc phase and that the giant planets that formed in these simulations match the observed giant planet eccentricity distribution best. These simulations also show that massive giant planets (above Jupiter mass) on eccentric orbits are less likely to host inner super-Earths as they get lost during the scattering phase, while systems with less massive giant planets on nearly circular orbits should harbour systems of inner super-Earths. Finally, our simulations predict that giant planets are not single, on average, but they live in multi-planet systems.

Unified Simulations of Planetary Formation and Atmospheric Evolution. II. Rapid Disk Clearing by Photoevaporation Yields Low-mass Super-Earth Atmospheres

Abstract Super-Earths possess low-mass H2/He atmospheres (typically less than 10% by mass). However, the origins of super-Earth atmospheres have not yet been ascertained. We investigate the role of rapid disk clearing by photoevaporation during the formation of super-Earths and their atmospheres. We perform unified simulations of super-Earth formation and atmospheric evolution in evolving disks that consider both photoevaporative winds and magnetically driven disk winds. For the growth mode of planetary cores, we consider two cases in which planetary embryos grow with and without pebble accretion. Our main findings are summarized as follows. (i) The time span of atmospheric accretion is shortened by rapid disk dissipation due to photoevaporation, which prevents super-Earth cores from accreting massive atmospheres. (ii) Even if planetary cores grow rapidly by embryo accretion in the case without pebble accretion, the onset of runaway gas accretion is delayed because the isolation mass for embryo accretion is small. Together with rapid disk clearing, the accretion of massive atmospheres can be avoided. (iii) After rapid disk clearing, a number of high-eccentricity embryos can remain in outer orbits. Thereafter, such embryos may collide with the super-Earths, leading to efficient impact erosion of accreted atmospheres. Therefore, we find that super-Earths with low-mass H2/He atmospheres are naturally produced by N-body simulations that consider realistic disk evolution.

The aeolian-erosion barrier for the growth of metre-size objects in protoplanetary discs

ABSTRACT Aeolian erosion is a destructive process that can erode small-size planetary objects through their interaction with a gaseous environment. Aeolian erosion operates in a wide range of environments and under various conditions. Aeolian erosion has been extensively explored in the context of geophysics in terrestrial planets. Here we show that aeolian erosion of cobbles, boulders, and small planetesimals in protoplanetary discs can constitute a significant barrier for the early stages of planet formation. We use analytic calculations to show that under the conditions prevailing in protoplanetary discs small bodies ($10\!-\!10^4 \, \rm {m}$) are highly susceptible to gas-drag aeolian erosion. At this size-range aeolian erosion can efficiently erode the planetesimals down to tens-cm size and quench any further growth of such small bodies. It thereby raises potential difficulties for channels suggested to alleviate the metre-size barrier. Nevertheless, the population of ∼decimetre-size cobbles resulting from aeolian erosion might boost the growth of larger (>km size) planetesimals and planetary embryos through increasing the efficiency of pebble-accretion, once/if such large planetesimals and planetary embryos exist in the disc.

Loss and Fractionation of Noble Gas Isotopes and Moderately Volatile Elements from Planetary Embryos and Early Venus, Earth and Mars

Here we discuss the current state of knowledge on how atmospheric escape processes can fractionate noble gas isotopes and moderately volatile rock-forming elements that populate primordial atmospheres, magma ocean related environments, and catastrophically outgassed steam atmospheres. Variations of isotopes and volatile elements in different planetary reservoirs keep information about atmospheric escape, composition and even the source of accreting material. We summarize our knowledge on atmospheric isotope ratios and discuss the latest evidence that proto-Venus and Earth captured small H$_2$-dominated primordial atmospheres that were lost by EUV-driven hydrodynamic escape after the disk dispersed. All relevant thermal and non-thermal atmospheric escape processes that can fractionate various isotopes and volatile elements are discussed. Erosion of early atmospheres, crust and mantle by large planetary impactors are also addressed. Further, we discuss how moderately volatile elements such as the radioactive heat producing element $^{40}$K and other rock-forming elements such as Mg can also be outgassed and lost from magma oceans that originate on large planetary embryos and accreting planets. Outgassed elements escape from planetary embryos with masses that are $\geq$\,M$_{\rm Moon}$ directly, or due to hydrodynamic drag of escaping H atoms originating from primordial- or steam atmospheres at more massive embryos. We discuss how these processes affect the final elemental composition and ratios such as K/U, Fe/Mg of early planets and their building blocks. Finally, we review modeling efforts that constrain the early evolution of Venus, Earth and Mars by reproducing their measured present day atmospheric $^{36}$Ar/$^{38}$Ar, $^{20}$Ne/$^{22}$Ne noble gas isotope ratios and the role of isotopes on the loss of water and its connection to the redox state on early Mars.

A new and simple prescription for planet orbital migration and eccentricity damping by planet–disc interactions based on dynamical friction

ABSTRACT During planet formation, gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet–disc interaction force and the time-scales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet–disc interactions based on dynamical friction, which can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence, it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

Tidal and general relativistic effects in rocky planet formation at the substellar mass limit using N-body simulations

Context. Recent observational results show that very low mass stars and brown dwarfs are able to host close-in rocky planets. Low-mass stars are the most abundant stars in the Galaxy, and the formation efficiency of their planetary systems is relevant in the computation of a global probability of finding Earth-like planets inside habitable zones. Tidal forces and relativistic effects are relevant in the latest dynamical evolution of planets around low-mass stars, and their effect on the planetary formation efficiency still needs to be addressed. Aims. Our goal is to evaluate the impact of tidal forces and relativistic effects on the formation of rocky planets around a star close to the substellar mass limit in terms of the resulting planetary architectures and its distribution according to the corresponding evolving habitable zone. Methods. We performed a set of N-body simulations spanning the first 100 Myr of the evolution of two systems composed of 224 embryos with a total mass 0.25 M⊕ and 74 embryos with a total mass 3 M⊕ around a central object of 0.08 M⊙. For these two scenarios we compared the planetary architectures that result from simulations that are purely gravitational with those from simulations that include the early contraction and spin-up of the central object, the distortions and dissipation tidal terms, and general relativistic effects. Results. We found that including these effects allows the formation and survival of a close-in (r < 0.07 au) population of rocky planets with masses in the range 0.001 < m∕M⊕ < 0.02 in all the simulations of the less massive scenario, and a close-in population with masses m ~ 0.35 M⊕ in just a few of the simulations of the more massive scenario. The surviving close-in bodies suffered more collisions during the integration time of the simulations. These collisions play an important role in their final masses. However, all of these bodies conserved their initial amount of water in mass throughout the integration time. Conclusions. The incorporation of tidal and general relativistic effects allows the formation of an in situ close-in population located in the habitable zone of the system. This means that both effects are relevant during the formation of rocky planets and their early evolution around stars close to the substellar mass limit, in particular when low-mass planetary embryos are involved.

Escape of rock-forming volatile elements and noble gases from planetary embryos

Large planetesimals and planetary embryos ranging from several hundred to a few thousand kilometers can develop magma oceans through mutual collisions, gravitational energy, and the heating of short-lived radioactive elements. During their solidification after the dissipation of the disk a steam atmosphere will be catastrophically outgassed and may be lost efficiently via hydrodynamic escape, as long as it does not condense. The escaping H-atoms that originate from the dissociation of H2O and H2 will drag heavier trace elements like noble gases such as Ne and Ar and outgassed moderately volatile rock-forming elements such as K, Na, Si, Mg, etc. into space. Under consideration of various EUV flux evolution scenarios of young solar-type stars, we apply an upper atmosphere hydrodynamic escape model that includes the dragging of heavier species by escaping H-atoms. We investigate the atmospheric/elemental escape and fractionation from planetary embryos with masses of 1 MMoon, 0.5 MMars, 1 MMars, and 1.5 MMars at different orbital distances between the orbits of Venus and Mars. Our results indicate that the steam atmospheres and the embedded trace elements will be lost efficiently before they condense for masses ≤0.5 MMars and orbital distances up to 1 AU. For heavier embryos of up to 1.5 MMars almost all of the considered steam atmospheres can be lost within ≈12 Myr, which lies within the time frame of the formation of the first Martian protocrust after ≈20 Myr, i.e. for such steam atmospheres to be lost completely a shallow magma ocean must remain below the atmosphere, which might be achieved through frequent impacts onto the planetary embryo. The considered outgassed noble gases and rock-forming elements will be completely dragged away together with the steam atmospheres under the assumption that the trace elements will reach the thermosphere. For embryos with masses ≤MMoon the gravity is too weak for a dense atmosphere to build up for the high magma ocean related surface temperatures and all outgassed elements will escape immediately to space. For all considered planetary masses and orbits the loss rates of Ar and Ne are so high that there will be no fractionation of their isotopes. The studied planetary embryos, even though not isotopically fractionated, will therefore be severely depleted in noble gases and moderately volatile elements. Hydrodynamic escape might then also affect the final composition of terrestrial planets that accrete out of such planetary embryos, such as the volatile content and the Fe/Mg ratio of a planet.

Earth-size planet formation in the habitable zone of circumbinary stars

ABSTRACT This work investigates the possibility of close binary (CB) star systems having Earth-size planets within their habitable zones (HZs). First, we selected all known CB systems with confirmed planets (totaling 22 systems) to calculate the boundaries of their respective HZs. However, only eight systems had all the data necessary for the computation of HZ. Then, we numerically explored the stability within HZs for each one of the eight systems using test particles. From the results, we selected five systems that have stable regions inside HZs, namely Kepler-34,35,38,413, and 453. For these five cases of systems with stable regions in HZ, we perform a series of numerical simulations for planet formation considering discs composed of planetary embryos and planetesimals, with two distinct density profiles, in addition to the stars and host planets of each system. We found that in the case of the Kepler-34 and 453 systems, no Earth-size planet is formed within HZs. Although planets with Earth-like masses were formed in Kepler-453, they were outside HZ. In contrast, for the Kepler-35 and 38 systems, the results showed that potentially habitable planets are formed in all simulations. In the case of the Kepler-413system, in just one simulation, a terrestrial planet was formed within HZ.

Spitzer’s debris disk legacy from main-sequence stars to white dwarfs

The Spitzer Space Telescope enabled the detection and characterization of infrared excess — attributed to thermal emission from circumstellar dust — towards more than 1,000 main-sequence stars and more than 40 white dwarfs. These systems have been dubbed ‘debris disks’ because they are believed to be planetary systems in which the observed dust is generated by collisions among otherwise undetected planetesimals, planetary embryos and planets. With its superb sensitivity and good spatial resolution, Spitzer revolutionized the study of debris disks, increasing the number of known systems by an order of magnitude and providing detailed disk demographics and properties for a substantial sample for the first time. In doing so, Spitzer set the stage for detailed characterization of the dust using Herschel, the Atacama Large Millimeter/submillimeter Array and high-contrast imaging instruments such as the Hubble Space Telescope’s NICMOS and STIS, the Very Large Telescope’s SPHERE, and the Gemini Planet Imager. Spitzer played a vital role in demonstrating that planetary systems are present and active around white dwarfs. Specifically, it showed that the origin of pollution in white dwarfs is circumstellar material. Finally, Spitzer’s longevity enabled time-domain studies of young debris disks that are actively forming terrestrial planets and rejuvenized the study of dusty white dwarfs. Spitzer’s legacy will endure into the next decade with discoveries expected from new instruments on SOFIA, the James Webb Space Telescope and SPHEREx. The Spitzer Space Telescope made huge advances in the study of debris disks around main-sequence stars and white dwarfs, increasing their number by an order of magnitude, and leading the way for the next generation of space-based infrared missions.

Realistic collisional water transport during terrestrial planet formation

Context. According to the latest theoretical and isotopic evidence, Earth’s water content originates mainly from today’s asteroid belt region, or at least from the same precursor material. This suggests that water was transported inwards to Earth, and to similar planets in their habitable zone, via (giant) collisions of planetary embryos and planetesimals during the chaotic final phase of planet formation. Aims. In current dynamical simulations water delivery to terrestrial planets is still studied almost exclusively by assuming oversimplified perfect merging, even though water and other volatiles are particularly prone to collisional transfer and loss. To close this gap we have developed a computational framework to model collisional water transport by direct combination of long-term N-body computations with dedicated 3D smooth particle hydrodynamics (SPH) collision simulations of differentiated, self-gravitating bodies for each event. Methods. Post-collision water inventories are traced self-consistently in the further dynamical evolution, in accretionary or erosive as well as hit-and-run encounters with two large surviving bodies, where besides collisional losses, water transfer between the encountering bodies has to be considered. This hybrid approach enables us for the first time to trace the full dynamical and collisional evolution of a system of approximately 200 bodies throughout the whole late-stage accretion phase (several hundred Myr). As a first application we choose a Solar System-like architecture with already formed giant planets on either circular or eccentric orbits and a debris disk spanning the whole terrestrial planet region (0.5–4 au). Results. Including realistic collision treatment leads to considerably different results than simple perfect merging, with lower mass planets and water inventories reduced regularly by a factor of two or more. Due to a combination of collisional losses and a considerably lengthened accretion phase, final water content, especially with giant planets on circular orbits, is strongly reduced to more Earth-like values, and closer to results with eccentric giant planets. Water delivery to potentially habitable planets is dominated by very few decisive collisions, mostly with embryo-sized or larger objects and only rarely with smaller bodies, at least if embryos have formed throughout the whole disk initially. The high frequency of hit-and-run collisions and the differences to predominantly accretionary encounters, such as generally low water (and mass) transfer efficiencies, are a crucial part of water delivery, and of system-wide evolution in general.

Dynamical evidence for an early giant planet instability

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. Here 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 (inward and outward) 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. About half of the planetary systems that we consider have a self-triggered instability. Of those, the median instability time is ∼4Myr. 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.

Catastrophic Events in Protoplanetary Disks and Their Observational Manifestations

Abstract Observations of protoplanetary disks with high angular resolution using an ALMA interferometer showed that ring-shaped structures are often visible in their images, indicating strong disturbances in the disks. The mechanisms of their formation are vividly discussed in the literature. This article shows that the formation of such structures can be the result of destructive collisions of large bodies (planetesimals and planetary embryos) accompanied by the formation of a large number of dust particles, and the subsequent evolution of a cloud of dust formed in this way.

Physical properties of terrestrial planets and water delivery in the habitable zone using N-body simulations with fragmentation

Aims. The goal of this research is to study how the fragmentation of planetary embryos can affect the physical and dynamical properties of terrestrial planets around solar-type stars. Our study focuses on the formation and evolution of planets and water delivery in the habitable zone (HZ). We distinguish class A and class B HZ planets, which have an accretion seed initially located inside and beyond the snow line, respectively. Methods. We developed an N-body integrator that incorporates fragmentation and hit-and-run collisions, which is called D3 N-body code. From this, we performed 46 numerical simulations of planetary accretion in systems that host two gaseous giants similar to Jupiter and Saturn. We compared two sets of 23 N-body simulations, one of which includes a realistic collisional treatment and the other one models all impacts as perfect mergers. Results. The final masses of the HZ planets formed in runs with fragmentation are about 15–20% lower than those obtained without fragmentation. As for the class A HZ planets, those formed in simulations without fragmentation experience very significant increases in mass with respect to their initial values, while the growth of those produced in runs with fragmentation is less relevant. We remark that the fragments play a secondary role in the masses of the class A HZ planets, providing less than 30% of their final values. In runs without fragmentation, the final fraction of water of the class A HZ planets keeps the initial value since they do not accrete water-rich embryos. In runs with fragmentation, the final fraction of water of such planets strongly depends on the model used to distribute the water after each collision. The class B HZ planets do not show significant differences concerning their final water contents in runs with and without fragmentation. From this, we find that the collisional fragmentation is not a barrier to the survival of water worlds in the HZ.