Search

Rocky planets may acquire a primordial atmosphere by outgassing of volatiles from their magma ocean. The distribution of O between H$_2$O, CO and CO$_2$ in chemical equilibrium subsequently changes significantly with decreasing temperature. We explore here two chemical models: one where CH$_4$ and NH$_3$ are assumed to be irrevocably destroyed by photolysis, and one where these molecules persist. In the first case, we show that CO cannot co-exist with H$_2$O, since CO oxidizes at low temperatures to form CO$_2$ and H$_2$. In both cases, H escapes from the thermosphere within a few ten million years by absorption of stellar XUV radiation. This escape drives an atmospheric self-oxidation process whereby rocky planet atmospheres become dominated by CO$_2$ and H$_2$O, regardless of their initial oxidation state at outgassing. HCN is considered a potential precursor of prebiotic compounds and RNA. Our oxidizing atmospheres are inefficient at producing HCN by lightning. Instead, we demonstrate that lightning-produced NO, which dissolves as nitrate in the oceans, and interplanetary dust particles may be the main sources of fixed nitrogen to emerging biospheres. Our results highlight the need for origin-of-life scenarios where the first metabolism fixes its C from CO$_2$, rather than from HCN and CO., Comment: Accepted for publication in Astrobiology

We present here results of numerical simulations of the formation and early evolution of rocky planets through pebble accretion, with an with an emphasis on hydrogen envelope longevity and the composition of the outgassed atmosphere. We model planets with a range in mass from 0.1 to 5 Earth masses that orbit between 0.7 and 1.7 AU. The composition of the outgassed atmosphere is calculated with the partial pressure of free oxygen fit to geophysical models of magma ocean self-oxidation. XUV radiation powered photoevaporation is considered as the main driver of atmospheric escape. We model planets that remain below the pebble isolation mass and hence accrete tenuous envelopes only. We consider slow, medium or fast initial stellar rotation for the temporal evolution of the XUV flux. The loss of the envelope is a key event that allows the magma ocean to crystallise and outgas its bulk volatiles. The atmospheric composition of the majority of our simulated planets is dominated by CO$_2$. Our planets accrete a total of 11.6 Earth oceans of water, the majority of which enters the core. The hydrospheres of planets lighter than the Earth reach several times the mass of the Earth's modern oceans, while the hydrospheres of planets ranging from 1 to 3.5 Earth masses are comparable to those of our planet. However, planets of 4-5 Earth masses have smaller hydrospheres due to trapping of volatiles in their massive mantles. Overall, our simulations demonstrate that hydrogen envelopes are easily lost from rocky planets and that this envelope loss triggers the most primordial partitioning of volatiles between the solid mantle and the atmosphere., Comment: Accepted for Astronomy & Astrophysics

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.

During their formation, nascent planetary systems are subject to FU Orionis outbursts that heat a substantial part of the disc. This causes water ice in the affected part of the disc to sublimate as the ice line moves outwards to several to tens of astronomical units. In this paper, we investigate how the subsequent cooling of the disc impacts the particle sizes. We calculate the resulting particle sizes in a disc model with cooling times between 100 and 1,000 years, corresponding to typical FU Orionis outbursts. As the disc cools and the ice line retreats inwards, water vapour forms icy mantles on existing silicate particles. This process is called heterogeneous nucleation. The nucleation rate per surface area of silicate substrate strongly depends on the degree of super-saturation of the water vapour in the gas. Fast cooling results in high super-saturation levels, high nucleation rates, and limited condensation growth because the main ice budget is spent in the nucleation. Slow cooling, on the other hand, leads to rare ice nucleation and efficient growth of ice-nucleated particles by subsequent condensation. We demonstrate that close to the quiescent ice line, pebbles with a size of about centimetres to decimetres form by this process. The largest of these are expected to undergo cracking collisions. However, their Stokes numbers still reach values that are high enough to potentially trigger planetesimal formation by the streaming instability if the background turbulence is weak. Stellar outbursts may thus promote planetesimal formation around the water ice line in protoplanetary discs., Comment: Accepted for publication in Astronomy & Astrophysics. 8 pages, 6 figures

The architecture and masses of planetary systems in the habitable zone could be strongly influenced by outer giant planets, if present. We investigate here the impact of outer giants on terrestrial planet formation, under the assumption that the final assembly of the planetary system is set by a giant impact phase. Utilizing a state-of-the-art N-body simulation software, GENGA, we interpret how the late stage of terrestrial planet formation results in diversity within planetary systems. We design two global model setups: in the first we place a gas giant on the outer side of planetesimals and embryos disk, while the other only has planetesimals and embryos but no giant. For the model including the outer giant, we study the effect of different giant initial masses, in the range 1.0-3.0 Jupiter mass, and orbital radii, in the range 2.0-5.8 AU.We also study the influence of different initial positions of planetesimals and embryos on the results. Our N-body simulation time is approximately 50 Myr. The results show that the existence of outer giant will promote the interaction between planetesimals and embryos, making the orbits of the formed terrestrial planets more compact, but placing the giant planet too close to the planetesimals and embryos disk suppresses the formation of massive rocky planets. In addition, under the classical theory, where planetary embryos and planetesimals collide to form terrestrial planets, our results show that the presence of a giant planet actually decreases the gap complexity of the inner planetary system., Comment: 12 pages, 15 figures

The gaseous envelope of an accreting rocky planet becomes hot enough to sublimate silicates and other refractory minerals. For this work, we studied the effect of the resulting envelope enrichment with a heavy vapor species on the composition and temperature of the envelope. For simplification, we used the gas-phase molecule SiO to represent the sublimation of silicate material. We solved the equilibrium structure equations in 1D for planets in the mass range of $0.1$ to $3\,M_\oplus$. The convective stability criterion was extended to take the stabilizing effect of the condensation of SiO clouds into account. We assumed that the envelope is both in hydrostatic equilibrium and in vapor equilibrium with the underlying magma ocean. This means that pebbles do not undergo sublimation in the envelope and therefore survive until they plunge into the magma ocean. We find that the emergence of an inner radiative region, where SiO condensation suppresses convection, increases the pressure and temperature in the inner envelope compared to pure H$_2$/He envelopes once $M_\mathrm{pl} \gtrsim 0.3\,M_\oplus$. For $M_\mathrm{pl}>0.75\,M_\oplus$, the temperature and pressure close to the surface reach the supercritical point of SiO. The amount of SiO stored in the envelope is lower than the total planet mass for low mass planets. However, for $M_\mathrm{pl}>2.0\,M_\oplus$, all accreted pebble material must contribute to maintain the vapor equilibrium in the envelope. Therefore, the non-vapor mass of the planet ceases to increase beyond this threshold. Overall, our vapor equilibrium model of the planetary envelope allows for direct core growth by pebble accretion up to much higher masses than previously thought., Comment: 10 pages, 8 figures, accepted for publication in A&A

During protoplanetary disk formation, dust grains located in the outer disk retain their pristine icy composition, while solids in the inner stellar-heated disk undergo volatile loss. This process may have left a fossil record in Solar System material showing different nucelosynthetic imprints that have been attributed to different degrees of thermal processing. However, it remains unclear how a large mass fraction of thermally-processed inner-disk pebbles is produced and how these grains are subsequently transported throughout the disk. In this work we numerically investigate the evolution in time of a two-component pebble disk, consisting of pristine pebbles and those that underwent ice sublimation. We find that stellar outbursts exceeding 1000 times the solar luminosity are efficient in thermally altering, through ice sublimation, a large mass fraction of pebbles (around 80%). After the establishment of this initial radial dust composition gradient throughout the disk, the subsequent mixing and inward drift of pristine outer-disk pebbles alter the inner disk bulk composition from processed to more unprocessed in time. Therefore, if processed pebbles without ice mantles have an isotopic composition similar to ureilite meteorites from the inner Solar System, inner-disk minor bodies forming from the early pebble flux (<1Myr) will be isotopically ureilite-like, while later-formed bodies will be increasingly admixed with the signature of the late incoming CI chondrite-like unprocessed pebbles. This appears to be largely consistent with the trend seen between the accretion age of different meteoric classes and their different stable isotope composition anomalies (in $\mu$54Cr, $\mu$48Ca, $\mu$30Si, $\mu$58Ni), but further work may be needed to explain the role of isotopically anomalous refractory inclusions and anomaly trends in other elements., Comment: 15 pages, 10 figures, 2 tables, 5 appendices (4 additional pages, 6 additional figures). Accepted for publication in A&A

Planetesimals are believed to form by the gravitational collapse of aerodynamically concentrated clumps of pebbles. Many properties of the objects in the cold classical Kuiper belt -- such as binarity, rotation, and size distribution -- are in agreement with this gravitational collapse model. Further support comes from the pebble-pile structure inferred for comet nuclei. For this study, we simulated the final assembly of a planetesimal from the gravitational collapse of a rotating clump of pebbles. We implemented a numerical method from granular dynamics to follow the collapse that includes the transition from a pebble swarm to solid cells at a high density. We compared the shapes of the simulated planetesimals with the shapes of the lobes of contact binaries and bilobed Solar System objects. We find that the gravitational collapse of slowly rotating pebble clouds naturally explains the formation of flattened ellipsoidal bodies. This result agrees well with the flattened structure of the bilobed planetesimal Arrokoth and the shapes of the components of bilobed comets., Comment: accepted for publication in A&A

We propose a pebble-driven core accretion scenario to explain the formation of giant planets around the late-M dwarfs of $M_{\star}{=}0.1{-}0.2 \ M_{\odot}$. In order to explore the optimal disk conditions for giant planet, we perform N-body simulations to investigate the growth and dynamical evolution of both single and multiple protoplanets in the disks with both inner viscously heated and outer stellar irradiated regions. The initial masses of the protoplanets are either assumed to be equal to $0.01 \ M_{\oplus}$ or calculated based on the formula derived from streaming instability simulations. Our findings indicate that massive planets are more likely to form in disks with longer lifetimes, higher solid masses, moderate to high levels of disk turbulence, and larger initial masses of protoplanets. In the single protoplanet growth cases, the highest planet core mass that can be reached is generally lower than the threshold necessary to trigger rapid gas accretion, which impedes the formation of giant planets. Nonetheless, in multi-protoplanet cases, the cores can exceed the pebble isolation mass barrier aided by frequent planet-planet collisions. This consequently speeds up their gas accretion and promotes giant planet formation, making the optimal parameter space to grow giant planets substantially wider. Taken together, our results suggest that even around very low-mass stellar hosts, the giant planets with orbital periods of ${\lesssim}100$ days are still likely to form when lunar-mass protoplanets first emerge from planetesimal accretion and then grow rapidly by a combination of pebble accretion and planet-planet collisions in disks with a high supply of pebble reservoir ${>}50 \ M_{\oplus}$ and turbulent level of $\alpha_{\rm t} {\sim} 10^{-3}{-}10^{-2}$., Comment: 22 pages, 13 figures, 3 tables, Accepted for publication in A&A

The presence of distant protoplanets may explain the observed gaps in the dust emission of protoplanetary disks. Here, we derive a novel analytical model to describe the temporal decay of the pebble flux through a protoplanetary disk as the result of radial drift. This has allowed us to investigate the growth and migration of distant protoplanets throughout the lifespan of the disk. We find that Moon-mass protoplanets that formed early on can grow to their pebble isolation mass, between approximately $20$ and $80\,M_{\oplus}$, within less than $1\,\mathrm{Myr,}$ in the $20$ to $80\,\mathrm{AU}$ region around solar-like stars. The subsequent fast migration in the early stages of gas accretion, after pebble accretion ends, transports these giant planets into their final orbits at $<\,$$10\,\mathrm{AU}$. However, our pebble decay model allows us to include a new pathway that may trigger the transition from pebble accretion to gas accretion after the pebble flux has decayed substantially. With this pebble decay pathway, we show that it is also possible to form gas giants beyond $10\,\mathrm{AU}$. The occurrence of these wide-orbit gas giants should be relatively low, since their core must attain sufficient mass to accrete gas before the pebble flux decays, while avoiding excessive migration. Since these gas giants do not reach the pebble isolation mass, their heavy element content is typically less than $10\,M_{\oplus}$. Our results imply that the observed gaps in protoplanetary disks could be caused by distant protoplanets that reached the pebble isolation mass and then migrated, while gas giants in wide orbits, such as PDS 70 b and c, accreted their gas after the decay in the pebble flux., Comment: 19 pages, 14 figures, accepted in A&A

As stellar compositions evolve over time in the Milky Way, so will the resulting planet populations. In order to place planet formation in the context of Galactic chemical evolution, we make use of a large ($N = 5\,325$) stellar sample representing the thin and thick discs, defined chemically, and the halo, and we simulate planet formation by pebble accretion around these stars. We build a chemical model of their protoplanetary discs, taking into account the relevant chemical transitions between vapour and refractory minerals, in order to track the resulting compositions of formed planets. We find that the masses of our synthetic planets increase on average with increasing stellar metallicity [Fe/H] and that giant planets and super-Earths are most common around thin-disc ($\alpha$-poor) stars since these stars have an overall higher budget of solid particles. Giant planets are found to be very rare ($\lesssim$1\%) around thick-disc ($\alpha$-rich) stars and nearly non-existent around halo stars. This indicates that the planet population is more diverse for more metal-rich stars in the thin disc. Water-rich planets are less common around low-metallicity stars since their low metallicity prohibits efficient growth beyond the water ice line. If we allow water to oxidise iron in the protoplanetary disc, this results in decreasing core mass fractions with increasing [Fe/H]. Excluding iron oxidation from our condensation model instead results in higher core mass fractions, in better agreement with the core-mass fraction of Earth, that increase with increasing [Fe/H]. Our work demonstrates how the Galactic chemical evolution and stellar parameters, such as stellar mass and chemical composition, can shape the resulting planet population., Comment: 21 pages, 16 figures, accepted in A&A

Protoplanets growing within the protoplanetary disk by pebble accretion acquire hydrostatic gas envelopes. Due to accretion heating, the temperature in these envelopes can become high enough to sublimate refractory minerals which are the major components of the accreted pebbles. Here we study the sublimation of different mineral species and determine whether sublimation plays a role during the growth by pebble accretion. For each snapshot in the growth process, we calculate the envelope structure and sublimation temperature of a set of mineral species representing different levels of volatility. Sublimation lines are determined using and equilibrium scheme for the chemical reactions responsible for destruction and formation of the relevant minerals. We find that the envelope of the growing planet reaches temperatures high enough to sublimate all considered mineral species when the mass is larger than 0.4 Earth masses. The sublimation lines are located within the gravitionally bound envelope of the planet. We make a detailed analysis of the sublimation of FeS at around 720 K, beyond which the mineral is attacked by H2 to form gaseous H2S and solid Fe. We calculate the sulfur concentration in the planet under the assumption that all sulfur released as H2S is lost from the planet by diffusion back to the protoplanetary disk. Our calculated values are in good agreement with the slightly depleted sulfur abundance of Mars, while the model overpredicts the extensive sulfur depletion of Earth by a factor of approximately 2. We show that a collision with a sulfur-rich body akin to Mars in the moon-forming impact lifts the Earth's sulfur abundance to approximately 10% of the solar value for all impactor masses above 0.05 Earth masses., Comment: 17 pages, 10 figures, accepted in A&A

The radial transport, or drift, of dust has taken a critical role in giant planet formation theory. However, it has been challenging to identify dust drift pile ups in the hard-to-observe inner disk. We find that the IM Lup disk shows evidence that it has been shaped by an episode of dust drift. Using radiative transfer and dust dynamical modeling we study the radial and vertical dust distribution. We find that high dust drift rates exceeding 110 M_earth/Myr are necessary to explain both the dust and CO observations. Furthermore, the bulk of the large dust present in the inner 20 au needs to be vertically extended, implying high turbulence alpha_z > 10^{-3} and small grains (0.2-1 mm). We suggest that this increased level of particle stirring is consistent with the inner dust-rich disk undergoing turbulence triggered by the vertical shear instability. The conditions in the IM Lup disk imply that giant planet formation through pebble accretion is only effective outside 20 au. If such an early, high turbulence inner region is a natural consequence of high dust drift rates, then this has major implications for understanding the formation regions of giant planets including Jupiter and Saturn., Comment: 12 pages, 9 figures, Accepted by Astrophysical Journal Letters

Pebble accretion is recognized as a significant accelerator of planet formation. Yet, only formulae for single-sized (monodisperse) distribution have been derived in the literature. These can lead to significant underestimates for Bondi accretion, for which the best accreted pebble size may not be the one that dominates the mass distribution. We derive in this paper the polydisperse theory of pebble accretion. We consider a power-law distribution in pebble radius, and we find the resulting surface and volume number density distribution functions. We derive also the exact monodisperse analytical pebble accretion rate for which 3D and 2D accretion are limits. In addition, we find analytical solutions to the polydisperse 2D Hill and 3D Bondi limits. We integrate the polydisperse pebble accretion numerically for the MRN distribution, finding a slight decrease (by an exact factor 3/7) in the Hill regime compared to the monodisperse case. In contrast, in the Bondi regime, we find 1-2 orders of magnitude higher accretion rates compared to monodisperse, also extending the onset of pebble accretion to 1-2 order of magnitude lower in mass. We find Myr-timescales, within the disk lifetime, for Bondi accretion on top of planetary seeds of masses $10^{-6}-10^{-4} M_\oplus$, over a significant range of the parameter space. This mass range overlaps with the high mass end of the planetesimal initial mass function, and thus pebble accretion is possible directly following formation by streaming instability. This alleviates the need for mutual planetesimal collisions as a major contribution to planetary growth., Comment: 17 pages, 12 figures. Accepted to ApJ

Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$\mu m$ dust particles in self-gravitating young protoplanetary disks with various strengths of turbulent viscosity. We aim to understand the physical conditions that determine the formation and spatial distribution of pebbles when both disk self-gravity and turbulent viscosity can be concurrently at work. Methods. We perform the thin-disk hydrodynamics simulations of self-gravitating protoplanetary disks over an initial time period of 0.5 Myr using the FEOSAD code. Turbulent viscosity is parameterized in terms of the spatially and temporally constant $\alpha$-parameter, while the effects of gravitational instability on dust growth is accounted for by calculating the effective parameter $\alpha_{\rm GI}$. We consider the evolution of dust component including momentum exchange with gas, dust self-gravity, and also a simplified model of dust growth. Results. We find that the level of turbulent viscosity strongly affects the spatial distribution and total mass of pebbles in the disk. The $\alpha=10^{-2}$ model is viscosity-dominated, pebbles are completely absent, and dust-to-gas mass ratio deviates from the reference 1:100 value no more than by 30\% throughout the disk extent. On the contrary, the $\alpha=10^{-3}$ model and, especially, the $\alpha=10^{-4}$ model are dominated by gravitational instability. The effective parameter $\alpha+\alpha_{\rm GI}$ is now a strongly varying function of radial distance. As a consequence, a bottle neck effect develops in the innermost disk regions, which makes gas and dust accumulate in a ring-like structure. Abridged., Comment: Accepted for publication in Astronomy and Astrophysics

Capture into mean motion resonance (MMR) is an important dynamical mechanism as it shapes the final architecture of a planetary system. We simulate systems of two or three planets undergoing migration with varied initial parameters such as planetary mass and disk surface density and analyse the resulting resonant chains. In contrast to previous studies, our results show that the disk properties have the dominant impact on capture into mean motion resonance, while the total planetary mass barely affects the final system configuration as long as the planet does not open a gap in the disk. We confirm that the adiabatic resonant capture is the correct framework to understand the conditions leading to MMR formation, since its predictions are qualitatively similar to the numerical results. However, we find that the eccentricity damping can facilitate the capture in a given resonance. We find that under typical disk conditions, planets tend to be captured into 2:1 or 3:2 MMRs, which agrees well with the observed exoplanet MMRs. Our results predict two categories of systems: those that have uniform chains of wide resonances (2:1 or 3:2 MMRs) and those that have a more compact inner pair than the outer pair such as 4:3:2 chains. Both categories of resonant chains are present in observed exoplanet systems. On the other hand, chains with a wider inner pair than the outer one are very rare and emerge from stochastic capture. Our work here can be used to link current configuration of exoplanetary systems to the formation conditions within protoplanetary disks., Comment: Accepted in A&A

The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability., Comment: 17 pages, 12 figures, 2 tables; accepted for publication in A&A; comments welcome

We explore the growth of planetary embryos by planetesimal accretion up to and beyond the point where pebble accretion becomes efficient at the so-called Hill-transition mass. Both the transition mass and the characteristic mass of planetesimals formed by the streaming instability increase with increasing distance from the star. We developed a model for the growth of a large planetesimal (embryo) embedded in a population of smaller planetesimals formed in a filament by the streaming instability. The model includes in a self-consistent way the collisional mass growth of the embryo, the fragmentation of the planetesimals, the velocity evolution of all involved bodies, as well as the viscous spreading of the filament. We find that the embryo accretes all available material in the filament during the lifetime of the protoplanetary disc only in the inner regions of the disc. In contrast, we find little or no growth in the outer parts of the disc beyond 5--10 AU. Overall, our results demonstrate very long timescales for collisional growth of planetesimals in the regions of the protoplanetary disc where giant planets form. As such, in order to form giant planets in cold orbits, pebble accretion must act directly on the largest bodies present in the initial mass-function of planetesimals with little or no help from mutual collisions., Comment: 12 pages, 11 figures, accepted for publication in A&A

Volatile molecules containing hydrogen, carbon, and nitrogen are key components of planetary atmospheres. In the pebble accretion model for rocky planet formation, these volatile species are accreted during the main planetary formation phase. For this study, we modelled the partitioning of volatiles within a growing planet and the outgassing to the surface. The core stores more than 90\% of the hydrogen and carbon budgets of Earth for realistic values of the partition coefficients of H and C between metal and silicate melts. The magma oceans of Earth and Venus are sufficiently deep to undergo oxidation of ferrous Fe$^{2+}$ to ferric Fe$^{3+}$. This increased oxidation state leads to the outgassing of primarily CO$_2$ and H$_2$O from the magma ocean of Earth. In contrast, the oxidation state of Mars' mantle remains low and the main outgassed hydrogen carrier is H$_2$. This hydrogen easily escapes the atmosphere due to the irradiation from the young Sun in XUV wavelengths, dragging with it the majority of the CO, CO$_2$, H$_2$O, and N$_2$ contents of the atmosphere. A small amount of surface water is maintained on Mars, in agreement with proposed ancient ocean shorelines, for moderately low values of the mantle oxidation. Nitrogen partitions relatively evenly between the core and the atmosphere due to its extremely low solubility in magma; the burial of large reservoirs of nitrogen in the core is thus not possible. The overall low N contents of Earth disagree with the high abundance of N in all chondrite classes and favours a volatile delivery by pebble snow. Our model of rapid rocky planet formation by pebble accretion displays broad consistency with the volatile contents of the Sun's terrestrial planets. The diversity of the terrestrial planets can therefore be used as benchmark cases to calibrate models of extrasolar rocky planets and their atmospheres., Comment: Version accepted for Astronomy & Astrophysics

We explore the heating and differentiation of rocky planets that grow by rapid pebble accretion. Our terrestrial planets grow outside of the ice line and initially accrete 28\% water ice by mass. The accretion of water stops after the protoplanet reaches a mass of $0.01\,M_{\rm E}$ where the gas envelope becomes hot enough to sublimate the ice and transport the vapour back to the protoplanetary disc by recycling flows. The energy released by the decay of $^{26}$Al melts the accreted ice to form clay (phyllosilicates), oxidized iron (FeO), and a water surface layer with ten times the mass of Earth's modern oceans. The ocean--atmosphere system undergoes a run-away greenhouse effect after the effective accretion temperature crosses a threshold of around 300 K. The run-away greenhouse process vaporizes the water layer, thereby trapping the accretion heat and heating the surface to more than 6,000 K. This causes the upper part of the mantle to melt and form a global magma ocean. Metal melt separates from silicate melt and sediments towards the bottom of the magma ocean; the gravitational energy released by the sedimentation leads to positive feedback where the beginning differentiation of the planet causes the whole mantle to melt and differentiate. All rocky planets thus naturally experience a magma ocean stage. We demonstrate that Earth's small excess of $^{182}$W (the decay product of $^{182}$Hf) relative to the chondrites is consistent with such rapid core formation within 5 Myr followed by equilibration of the W reservoir in Earth's mantle with $^{182}$W-poor material from the core of a planetary-mass impactor, provided that the equilibration degree is at least 25%-50%, depending on the initial Hf/W ratio. The planetary collision must have occurred at least 35 Myr after the main accretion phase of the terrestrial planets., Comment: Version accepted for Astronomy & Astrophysics

We present a series of papers dedicated to modelling the accretion and differentiation of rocky planets that form by pebble accretion within the lifetime of the protoplanetary disc. In this first paper, we focus on how the accreted ice determines the distribution of iron between the mantle (oxidized FeO and FeO$_{1.5}$) and the core (metallic Fe and FeS). We find that an initial primitive composition of ice-rich material leads, upon heating by the decay of $^{26}$Al, to extensive water flow and the formation of clay minerals inside planetesimals. Metallic iron dissolves in liquid water and precipitates as oxidized magnetite Fe$_3$O$_4$. Further heating by $^{26}$Al destabilizes the clay at a temperature of around 900 K. The released supercritical water ejects the entire water content from the planetesimal. Upon reaching the silicate melting temperature of 1,700 K, planetesimals further differentiate into a core (made mainly of iron sulfide FeS) and a mantle with a high fraction of oxidized iron. We propose that the asteroid Vesta's significant FeO fraction in the mantle is a testimony of its original ice content. We consider Vesta to be a surviving member of the population of protoplanets from which Mars, Earth, and Venus grew by pebble accretion. We show that the increase in the core mass fraction and decrease in FeO contents with increasing planetary mass (in the sequence Vesta -- Mars -- Earth) is naturally explained by the growth of terrestrial planets outside of the water ice line through accretion of pebbles containing iron that was dominantly in metallic form with an intrinsically low oxidation degree., Comment: Version accepted for Astronomy & Astrophysics

Aims. In this work, we aim to make a differential comparison of the neutron-capture and p-process element molybdenum (Mo) in the stellar populations in the local disk(s) and the bulge, focusing on minimising possible systematic effects in the analysis. Methods. The stellar sample consists of 45 bulge and 291 local disk K-giants, observed with high-resolution optical spectra. The abundances are determined by fitting synthetic spectra using the SME-code. The disk sample is separated into thin- and thick-disk components using a combination of abundances and kinematics. The cosmic origin of Mo is investigated and discussed by comparing with previous published abundances of Mo and the neutron-capture elements cerium (Ce) and europium (Eu). Results. We determine reliable Mo abundances for 35 bulge and 282 disk giants with a typical uncertainty of [Mo/Fe]~0.2 and ~0.1 dex for the bulge and disk, respectively. Conclusions. We find that the bulge possibly is enhanced in [Mo/Fe] compared to the thick disk, which we do not observe in either [Ce/Fe] nor [Eu/Fe]. This might suggest a higher past star-formation rate in the bulge, however, since we do not observe the bulge to be enhanced in [Eu/Fe], the origin of the molybdenum enhancement is yet to be constrained. Although, the scatter is large, we may be observing evidence of the p-process contributing to the heavy element production in the chemical evolution of the bulge., Comment: Accepted for publication in A&A. 12 pages, 8 figures

Streaming instability is hypothesized to be triggered at particular protoplanetary disk locations where the volume density of the solid particles is enriched comparable to that of the gas. A ring of planetesimals thus forms when this condition is fulfilled locally. These planetesimals collide with each other and accrete inward drifting pebbles from the outer disk to further increase masses. We investigate the growth of the planetesimals that form in a ring-belt at various disk radii. Their initial mass distributions are calculated based on the formula summarized from the streaming instability simulations. We simulate the subsequent dynamical evolution of the planetesimals with a protoplanetary disk model based either on the minimum mass solar nebula (MMSN) or on the Toomre stability criterion. For the MMSN model, both pebble accretion and planetesimal accretion are efficient at a close-in orbit of $0.3$ AU, resulting in the emergence of several super-Earth mass planets after $1$ Myr. For comparison, only the most massive planetesimals undergo substantial mass growth when they are born at $r{=}3$ AU, while the planetesimals at $r{=}30$ AU experience little or no growth. On the other hand, in the denser Toomre disk, the most massive forming planets can reach Earth mass at $t{=}1$ Myr and reach a mass between that of Neptune and that of Saturn within $3$ Myr at $30$ AU and $100$ AU. Both the pebble and planetesimal accretion rate decrease with disk radial distance. Nevertheless, planetesimal accretion is less pronounced than pebble accretion at more distant disk regions. Taken together, the planets acquire higher masses when the disk has a higher gas density, a higher pebble flux, and/or a lower Stokes number of pebbles., Comment: 11 pages, 7 figures, Accepted for publication in A&A

Meteorites display an isotopic composition dichotomy between non-carbonaceous (NC) and carbonaceous (CC) groups, indicating that planetesimal formation in the solar protoplanetary disk occurred in two distinct reservoirs. The prevailing view is that a rapidly formed Jupiter acted as a barrier between these reservoirs. We show a fundamental inconsistency in this model: if Jupiter is an efficient blocker of drifting pebbles, then the interior NC reservoir is depleted by radial drift within a few hundred thousand years. If Jupiter lets material pass it, then the two reservoirs will be mixed. Instead, we demonstrate that the arrival of the CC pebbles in the inner disk is delayed for several million years by the viscous expansion of the protoplanetary disk. Our results support that Jupiter formed in the outer disk (>10 AU) and allowed a considerable amount of CC material to pass it and become accreted by the terrestrial planets., Comment: Authors' version, 13 pages, 4 figures. Published in Science Advances in the issue of April 22, 2022

The pathway to forming the iron-rich planet Mercury remains mysterious. Mercury's core makes up 70% of the planetary mass, which implies a significant enrichment of iron relative to silicates, while its mantle is strongly depleted in oxidized iron. The high core mass fraction is traditionally ascribed to evaporative loss of silicates, e.g. following a giant impact, but the high abundance of moderately volatile elements in the mantle of Mercury is inconsistent with reaching temperatures much above 1,000 K during its formation. Here we explore the nucleation of solid particles from a gas of solar composition that cools down in the hot inner regions of the protoplanetary disc. The high surface tension of iron causes iron particles to nucleate homogeneously (i.e., not on a more refractory substrate) under very high supersaturation. The low nucleation rates lead to depositional growth of large iron pebbles on a sparse population of nucleated iron nano-particles. Silicates in the form of iron-free MgSiO$_3$ nucleate at similar temperatures but obtain smaller sizes due to the much higher number of nucleated particles. This results in a chemical separation of large iron particles from silicate particles with ten times lower Stokes numbers. We propose that such conditions lead to the formation of iron-rich planetesimals by the streaming instability. In this view, Mercury formed by accretion of iron-rich planetesimals with a sub-solar abundance of highly reduced silicate material. Our results imply that the iron-rich planets known to orbit the Sun and other stars are not required to have experienced mantle-stripping impacts. Instead their formation could be a direct consequence of temperature fluctuations in protoplanetary discs and chemical separation of distinct crystal species through the ensuing nucleation process., Comment: Accepted for publication in Astronomy & Astrophysics

Observations and models of giant planets indicate that such objects are enriched in heavy elements compared to solar abundances. The prevailing view is that giant planets accreted multiple Earth masses of heavy elements after the end of core formation. Such late solid enrichment is commonly explained by the accretion of planetesimals. Planetesimals are expected to form at the edges of planetary gaps, and here we address the question of whether these planetesimals can be accreted in large enough amounts to explain the inferred high heavy element contents of giant planets. We perform a series of N-body simulations of the dynamics of planetesimals and planets during the planetary growth phase, taking into account gas drag as well as the enhanced collision cross-section caused by the extended envelopes. We consider the growth of Jupiter and Saturn via gas accretion after reaching the pebble isolation mass and we include their migration in an evolving disk. We find that the accretion efficiency of planetesimals formed at planetary gap edges is very low: less than 10% of the formed planetesimals are accreted even in the most favorable cases, which in our model corresponds to a few Earth-masses. When planetesimals are assumed to form beyond the feeding zone of the planets, extending to a few Hill radii from a planet, accretion becomes negligible. Furthermore, we find that the accretion efficiency increases when the planetary migration distance is increased and that the efficiency does not increase when the planetesimal radii are decreased. Based on these results we conclude that it is difficult to explain the large heavy element content of giant planets with planetesimal accretion during the gas accretion phase. Alternative processes most likely are required, e.g. accretion of vapor deposited by drifting pebbles., Comment: Accepted for publication in A&A

Transformer models have shown great success handling long-range interactions, making them a promising tool for modeling video. However, they lack inductive biases and scale quadratically with input length. These limitations are further exacerbated when dealing with the high dimensionality introduced by the temporal dimension. While there are surveys analyzing the advances of Transformers for vision, none focus on an in-depth analysis of video-specific designs. In this survey, we analyze the main contributions and trends of works leveraging Transformers to model video. Specifically, we delve into how videos are handled at the input level first. Then, we study the architectural changes made to deal with video more efficiently, reduce redundancy, re-introduce useful inductive biases, and capture long-term temporal dynamics. In addition, we provide an overview of different training regimes and explore effective self-supervised learning strategies for video. Finally, we conduct a performance comparison on the most common benchmark for Video Transformers (i.e., action classification), finding them to outperform 3D ConvNets even with less computational complexity.

The streaming instability provides an efficient way of overcoming the growth barriers in the initial stages of the planet formation process. Considering the realistic case of a particle size distribution, the dynamics of the system is altered compared to the outcome of single size models. In order to understand the outcome of the multi-species streaming instability in detail, we perform a large parameter study in terms of particle number, particle size distribution, particle size range, initial metallicity and initial particle scale height. We study vertically stratified systems and determine the metallicity threshold for filament formation. We compare these with a system where the initial particle distribution is unstratified and find that its evolution follows that of its stratified counterpart. We find that change in particle number does not result in significant variation in the efficiency and timing of filament formation. We also see that there is no clear trend for how varying the size distribution in combination with particle size range changes the outcome of the multi-species streaming instability. Finally, we find that an initial metallicity of $Z_{\rm{init}}=0.005$ and $Z_{\rm{init}}=0.01$ both result in similar critical metallicity values for the start of filament formation. Our results show that the inclusion of a particle size distribution into streaming instability simulations, while changing the dynamics as compared to mono-disperse systems, does not result in overall unfavorable conditions for solid growth. We attribute the sub-dominant role of multiple species to the high-density conditions in the midplane, conditions under which also linear stability analysis predict little difference between single and multiple species., Comment: 13 pages, 13 figures, 2 tables, accepted to A\&A

Inward migration of giant planets is predicted by hydrodynamical simulations during the gas phase of the protoplanetary disc. The phenomenon is also invoked to explain resonant and near-resonant exoplanetary system structures. The early inward migration may also have affected our Solar System and sculpted its different minor planet reservoirs. In this study we explore how the early inward migration of the giant planets shapes the Kuiper Belt. We test different scenarios with only Neptune and Uranus and with all the four giant planets, including also some models with the subsequent outward planetesimal-driven migration of Neptune after the gas dispersal. We find objects populating mean motion resonances even when Neptune and Uranus do not migrate at all or only migrate inwards. When the planets are fixed, planetesimals stick only temporarily to the mean motion resonances, while inwards migration yields a new channel to populate the resonances without invoking convergent migration. In these cases, however, it is hard to populate mean motion resonances that do not cross the planetesimal disc (such as 2:1 and 5:2) and there is a lack of resonant KBOs that cross Neptune's orbit. These Neptune crossers are an unambiguous signature of the outward migration of Neptune. The starting position and the growth rate of Neptune matters for the contamination of the classical Kuiper belt region from neighbouring regions. The eccentricity and inclination space of the hot classicals and the scattered disc region become much more populated when all the giant planets are included. The 5:2 resonance with Neptune becomes increasingly populated with deeper inward migrations of Neptune. The overall inclination distribution, however, is still narrower than from observations, as is generally the case for Kuiper belt population models., Comment: 13 pages, 11 figures, 1 table; accepted for publication in Astronomy & Astrophysics

Understanding the origin of chondritic components and their accretion pathways is critical to unravel the magnitude of mass transport in the protoplanetary disk, the accretionary history of the terrestrial planet region and, by extension, its prebiotic inventory. Here, we trace the heritage of pristine components from the relatively unaltered CV chondrite Leoville through their mass-independent Cr and mass-dependent Zn isotope compositions. Investigating these chondritic fractions in such detail reveals an onion-shell structure of chondrules, which is characterized by 54Cr- and 66Zn-poor cores surrounded by increasingly 54Cr- and 66Zn-rich igneous rims and an outer coating of fine-grained dust. This is interpreted as a progressive addition of 54Cr- and 66Zn-rich, CI-like material to the accretion region of these carbonaceous chondrites. Our findings show that the observed Cr isotopic range in chondrules from more altered CV chondrites is the result of chemical equilibration between chondrules and matrix during secondary alteration. The 54Cr-poor nature of the cores of Leoville chondrules implies formation in the inner Solar System and subsequent massive outward chondrule transport past the Jupiter barrier. At the same time, CI-like dust is transferred inwards. We propose that the accreting Earth acquired CI-like dust through this mechanism within the lifetime of the disk. This radial mixing of chondrules and matrix shows the limited capacity of Jupiter to act as an efficient barrier and maintain the proposed non-carbonaceous and carbonaceous chondrite dichotomy over time. Finally, also considering current astrophysical models, we explore both inner and outer Solar System origins for the CV chondrite parent body., Comment: accepted after peer-review

Pebbles of millimeter sizes are abundant in protoplanetary discs around young stars. Chondrules inside primitive meteorites - formed by melting of dust aggregate pebbles or in impacts between planetesimals - have similar sizes. The role of pebble accretion for terrestrial planet formation is nevertheless unclear. Here we present a model where inwards-drifting pebbles feed the growth of terrestrial planets. The masses and orbits of Venus, Earth, Theia (which later collided with the Earth to form the Moon) and Mars are all consistent with pebble accretion onto protoplanets that formed around Mars' orbit and migrated to their final positions while growing. The isotopic compositions of Earth and Mars are matched qualitatively by accretion of two generations of pebbles, carrying distinct isotopic signatures. Finally, we show that the water and carbon budget of Earth can be delivered by pebbles from the early generation before the gas envelope became hot enough to vaporise volatiles., Comment: Science Advances, in press

We investigate the possibility of erosion of planetesimals in a protoplanetary disk. We use theory and direct numerical simulations (Lattice Boltzmann Method) to calculate the erosion of large -- much larger than the mean-free-path of gas molecules -- bodies of different shapes in flows. We find that erosion follows a universal power-law in time, at intermediate times, independent of the Reynolds number of the flow and the initial shape of the body. Consequently, we estimate that planetesimals in eccentric orbits, of even very small eccentricity, rapidly (in about hundred years) erodes away if the semi-major axis of their orbit lies in the inner disk -- less than about $10$ au. Even planetesimals in circular orbits erode away in approximately ten thousand years if the semi-major axis of their orbits are $\lessapprox 0.6$au., Comment: Accepted for publication in Astrophysical J

The dynamics of planetesimals plays an important role in planet formation, because their velocity distribution sets the growth rate to larger bodies. When planetesimals form in protoplanetary discs, their orbits are nearly circular and planar due to the effect of gas drag. However, mutual close encounters of the planetesimals increase eccentricities and inclinations until an equilibrium between stirring and damping is reached. After disc dissipation, there is no more gas drag and mutual close encounters as well as encounters with planets stir the orbits again. The high number of planetesimals in protoplanetary discs renders it difficult to simulate their dynamics by means of direct N-body simulations of planet formation. Therefore, we developed a novel method for the dynamical evolution of planetesimals that is based on following close encounters between planetesimal-mass bodies and gravitational stirring by planet-mass bodies. To separate the orbital motion from the close encounters, we employ a Hamiltonian splitting scheme as used in symplectic N-body integrators. Close encounters are identified using a cell algorithm with linear scaling in the number of bodies. A grouping algorithm is used to create small groups of interacting bodies which are integrated separately. Our method allows simulating a high number of planetesimals interacting through gravity and collisions with low computational cost. The typical computational time is of the order of minutes or hours, up to a few days for more complex simulations, as compared to several hours or even weeks for the same setup with full N-body. The dynamical evolution of the bodies is sufficiently well reproduced. This will make it possible to study the growth of planetesimals through collisions and pebble accretion coupled to their dynamics for a much higher number of bodies than previously accessible with full N-body simulations., Comment: accepted for publication in A&A

We analyse the size evolution of pebbles accreted into the gaseous envelope of a protoplanet growing in a protoplanetary disc, taking into account collisions driven by the relative sedimentation speed as well as the convective gas motion. Using a simple estimate of the convective gas speed based on the pebble accretion luminosity, we find that the speed of the convective gas is higher than the sedimentation speed for all particles smaller than 1 mm. This implies that both pebbles and pebble fragments are strongly affected by the convective gas motion and will be transported by large-scale convection cells both towards and away from the protoplanet's surface. We present a simple scheme for evolving the characteristic size of the pebbles, taking into account the effects of erosion, mass transfer and fragmentation. Including the downwards motion of convective cells for the transport of pebbles with an initial radius of 1 millimeter, we find pebble sizes between 100 microns and 1 millimeter near the surface of the protoplanet. These sizes are generally amenable to accretion at the base of the convection flow. Small protoplanets far from the star (>30 AU) nevertheless erode their pebbles to sizes below 10 microns; future hydrodynamical simulations will be needed to determine whether such small fragments can detach from the convection flow and become accreted by the protoplanet., Comment: Accepted for publication in The Astrophysical Journal

The dynamical stability of tightly packed exoplanetary systems remains poorly understood. While for a two-planet system a sharp stability boundary exists, numerical simulations of three and more planet systems show that they can experience instability on timescales up to billions of years. Moreover, an exponential trend between the planet orbital separation measured in units of Hill radii and the survival time has been reported. While these findings have been observed in numerous numerical simulations, little is known of the actual mechanism leading to instability. Contrary to a constant diffusion process, planetary systems seem to remain dynamically quiescent for most of their lifetime before a very short unstable phase. In this work, we show how the slow chaotic diffusion due to the overlap of three-body resonances dominates the timescale leading to the instability for initially coplanar and circular orbits. While the last instability phase is related to scattering due to two-planet mean motion resonances (MMR), for circular orbits the two-planets MMR are too far separated to destabilize systems initially away from them. We develop an analytical model to generalize the empirical trend obtained for equal mass and equally-spaced planets to general systems. We obtain an analytical estimate of the survival time consistent with simulations over four orders of magnitude for the planet to star mass ratio $\epsilon$, and 6 to 8 orders of magnitude for the instability time. We also confirm that measuring the orbital spacing in terms of Hill radii is not adapted and that the right spacing unit scales as $\epsilon^{1/4}$. We predict that beyond a certain spacing, the three-planet resonances are not overlapped, which results in an increase of the survival time. We finally discuss the extension of our result to more general systems, containing more planets on initially non circular orbits., Comment: Accepted in A&A. Added numerical tests with test particles and non uniform systems. A notebook as well as python scripts to reproduce the article figure are available at https://github.com/acpetit/PlanetSysSurvivalTime

Single-binary scattering may lead to an exchange where the single object captures a component of the binary, forming a new binary. This has been well studied in encounters between a star--planet pair and a single star. Here we explore the application of the exchange mechanism to a planet--satellite pair and another planet in the gravitational potential of a central star. As a case study, we focus on encounters between a satellite-bearing object and Neptune. We investigate whether Neptune can capture satellites from that object and if the captured satellites have orbits analogous to the Neptunian moons Triton and Nereid. Using $N$-body simulations, we study the capture probability at different encounter distances. Post-capture, we use a simple analytical argument to estimate how the captured orbits evolve under collisional and tidal effects. We find that the average capture probability reaches $\sim$$10\%$ if Neptune penetrates the donor planet's satellite system. Most moons grabbed by Neptune acquire highly eccentric orbits. Post-capture, around half of those captured, especially those on tight orbits, can be circularised, either by tides only or by collisions+tides, turning into Triton-like objects. Captures further out, on the other hand, stay on wide and eccentric orbits like that of Nereid. Both moon types can be captured in the same encounter and they have wide distributions in orbital inclination. Therefore, Triton naturally has a $\sim$50\% chance of being retrograde. A similar process potentially applies to an exoplanetary system, and our model predicts that exomoons can jump from one planet to another during planetary scattering. Specifically, there should be two distinct populations of captured moons: one on close-in circular orbits and the other on far-out eccentric orbits. The two populations may have highly inclined prograde or retrograde orbits., Comment: A&A in press

The first stages of planet formation take place in protoplanetary disks that are largely made up of gas. Understanding how the gas affects planetesimals in the protoplanetary disk is therefore essential. In this paper, we discuss whether or not gas flow can erode planetesimals. We estimate how much shear stress is exerted onto the planetesimal surface by the gas as a function of disk and planetesimal properties. To determine whether erosion can take place, we compare this with previous measurements of the critical stress that a pebble-pile planetesimal can withstand before erosion begins. If erosion takes place, we estimate the erosion time of the affected planetesimals. We also illustrate our estimates with two-dimensional numerical simulations of flows around planetesimals using the lattice Boltzmann method. We find that the wall shear stress can overcome the critical stress of planetesimals in an eccentric orbit within the innermost regions of the disk. The high eccentricities needed to reach erosive stresses could be the result of shepherding by migrating planets. We also find that if a planetesimal erodes, it does so on short timescales. For planetesimals residing outside of $1 \ \rm{au}$, we find that they are mainly safe from erosion, even in the case of highly eccentric orbits.

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$)., Comment: 12 pages, 8 figures, accepted in A&A

Recent surveys have revealed that protoplanetary discs typically have dust masses that appear to be insufficient to account for the high occurrence rate of exoplanet systems. We demonstrate that this observed dust depletion is consistent with the radial drift of pebbles. Using a Monte Carlo method we simulate the evolution of a cluster of protoplanetary discs, using a 1D numerical method to viscously evolve each gas disc together with the radial drift of dust particles that have grown to 100 $\mu$m in size. For a 2 Myr old cluster of stars, we find a slightly sub-linear scaling between the gas disc mass and the gas accretion rate ($M_\mathrm{g}\propto\dot{M}^{0.9}$). However, for the dust mass we find that evolved dust discs have a much weaker scaling with the gas accretion rate, with the precise scaling depending on the age at which the cluster is sampled and the intrinsic age spread of the discs in the cluster. Ultimately, we find that the dust mass present in protoplanetary disc is on the order of 10-100 Earth masses in 1-3 Myr old star-forming regions, a factor of 10 to 100 depleted from the original dust budget. As the dust drains from the outer disc, pebbles pile up in the inner disc and locally increase the dust-to-gas ratio by a factor of up to 4 above the initial value. In these high dust-to-gas ratio regions we find conditions that are favourable for planetesimal formation via the streaming instability and subsequent growth by pebble accretion. We also find the following scaling relations with stellar mass within a 1-2 Myr old cluster: a slightly super-linear scaling between the gas accretion rate and stellar mass ($\dot{M}\propto M_\star^{1.4}$), a slightly super-linear scaling between the gas disc mass and the stellar mass ($M_\mathrm{g}\propto M_\star^{1.4}$) and a super-linear relation between the dust disc mass and stellar mass ($M_\mathrm{d}\propto M_\star^{1.4-4.1}$)., Comment: 18 pages, 18 figures, accepted for publication in A&A

We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is $30\rm~ms^{-1}$. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula (MMSN) regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset ($<10^5$ yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. (Abridged), Comment: 13 pages, 10 figures

Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting., Comment: 19 pages, 18 figures, Accepted to A&A