Showing 1–50 of 137 results for author: Johansen, A

The PLATO Mission

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

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

Submitted 8 June, 2024; originally announced June 2024.

Fast formation of large ice pebbles after FU Orionis outbursts

Authors: Katrin Ros, Anders Johansen

Abstract: 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 particl… ▽ More 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. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 8 pages, 6 figures

Journal ref: A&A 686, A237 (2024)

How the presence of a giant planet affects the outcome of terrestrial planet formation simulations

Authors: Zhihui Kong, Anders Johansen, Michiel Lambrechts, Jonathan H. Jiang, Zong-Hong Zhu

Abstract: 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 th… ▽ More 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. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 12 pages, 15 figures

Journal ref: A&A 687, A121 (2024)

Vapor equilibrium models of accreting rocky planets demonstrate direct core growth by pebble accretion

Authors: Marie-Luise Steinmeyer, Anders Johansen

Abstract: 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 equi… ▽ More 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. △ Less

Submitted 16 February, 2024; v1 submitted 14 February, 2024; originally announced February 2024.

Comments: 10 pages, 8 figures, accepted for publication in A&A

Thermal processing of primordial pebbles in evolving protoplanetary disks

Authors: Maria Jose Colmenares, Michiel Lambrechts, Elishevah van Kooten, Anders Johansen

Abstract: 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 l… ▽ More 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 $μ$54Cr, $μ$48Ca, $μ$30Si, $μ$58Ni), but further work may be needed to explain the role of isotopically anomalous refractory inclusions and anomaly trends in other elements. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 15 pages, 10 figures, 2 tables, 5 appendices (4 additional pages, 6 additional figures). Accepted for publication in A&A

Journal ref: A&A 685, A114 (2024)

Formation of flattened planetesimals by gravitational collapse of rotating pebble clouds

Authors: Sebastian Lorek, Anders Johansen

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: accepted for publication in A&A

Forming Giant Planets Around Late-M Dwarfs: Pebble Accretion and Planet-Planet Collision

Authors: Mengrui Pan, Beibei Liu, Anders Johansen, Masahiro Ogihara, Su Wang, Jianghui Ji, Sharon X. Wang, Fabo Feng, Ignasi Riba

Abstract: 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… ▽ More 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 $α_{\rm t} {\sim} 10^{-3}{-}10^{-2}$. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 22 pages, 13 figures, 3 tables, Accepted for publication in A&A

Formation of wide-orbit giant planets in protoplanetary disks with a decreasing pebble flux

Authors: Nerea Gurrutxaga, Anders Johansen, Michiel Lambrechts, Johan Appelgren

Abstract: 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 M… ▽ More 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. △ Less

Submitted 6 February, 2024; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: 19 pages, 14 figures, accepted in A&A

Journal ref: A&A, 682, A43 (2024)

Planet formation throughout the Milky Way: Planet populations in the context of Galactic chemical evolution

Authors: Jesper Nielsen, Matthew Raymond Gent, Maria Bergemann, Philipp Eitner, Anders Johansen

Abstract: 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 chemica… ▽ More 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 ($α$-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 ($α$-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. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 21 pages, 16 figures, accepted in A&A

Journal ref: A&A 678, A74 (2023)

Sublimation of refractory minerals in the gas envelopes of accreting rocky planets

Authors: Marie-Luise Steinmeyer, Peter Woitke, Anders Johansen

Abstract: 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 gro… ▽ More 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. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 17 pages, 10 figures, accepted in A&A

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

A potential site for wide-orbit giant planet formation in the IM Lup disk

Authors: Arthur Bosman, Johan Appelgren, Edwin A. Bergin, Michiel Lambrechts, Anders Johansen

Abstract: 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… ▽ More 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. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 12 pages, 9 figures, Accepted by Astrophysical Journal Letters

An Analytical Theory for the Growth from Planetesimals to Planets by Polydisperse Pebble Accretion

Authors: Wladimir Lyra, Anders Johansen, Manuel H. Cañas, Chao-Chin Yang

Abstract: 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 pebb… ▽ More 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. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: 17 pages, 12 figures. Accepted to ApJ

Formation of pebbles in (gravito-)viscous protoplanetary disks with various turbulent strengths

Authors: Eduard I. Vorobyov, Vardan G. Elbakyan, Anders Johansen, Michiel Lambrechts, Aleksandr M. Skliarevskii, O. P. Stoyanovskaya

Abstract: Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μ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… ▽ More Aims. Dust plays a crucial role in the evolution of protoplanetary disks. We study the dynamics and growth of initially sub-$μ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 $α$-parameter, while the effects of gravitational instability on dust growth is accounted for by calculating the effective parameter $α_{\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 $α=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 $α=10^{-3}$ model and, especially, the $α=10^{-4}$ model are dominated by gravitational instability. The effective parameter $α+α_{\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. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Mean motion resonance capture in the context of type-I migration

Authors: Kaltrina Kajtazi, Antoine C. Petit, Anders Johansen

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 11 November, 2022; originally announced November 2022.

Comments: Accepted in A&A

Journal ref: A&A 669, A44 (2023)

The coexistence of the streaming instability and the vertical shear instability in protoplanetary disks: Planetesimal formation thresholds explored in two-dimensional global models

Authors: Urs Schäfer, Anders Johansen

Abstract: 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… ▽ More 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. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: 17 pages, 12 figures, 2 tables; accepted for publication in A&A; comments welcome

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

Growing the seeds of pebble accretion through planetesimal accretion

Authors: Sebastian Lorek, Anders Johansen

Abstract: 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)… ▽ More 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. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 12 pages, 11 figures, accepted for publication in A&A

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

Anatomy of rocky planets formed by rapid pebble accretion III. Partitioning of volatiles between planetary core, mantle, and atmosphere

Authors: Anders Johansen, Thomas Ronnet, Martin Schiller, Zhengbin Deng, Martin Bizzarro

Abstract: 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 hy… ▽ More 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. △ Less

Submitted 26 January, 2023; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Version accepted for Astronomy & Astrophysics

Journal ref: A&A 671, A76 (2023)

Anatomy of rocky planets formed by rapid pebble accretion II. Differentiation by accretion energy and thermal blanketing

Authors: Anders Johansen, Thomas Ronnet, Martin Schiller, Zhengbin Deng, Martin Bizzarro

Abstract: 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 di… ▽ More 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. △ Less

Submitted 26 January, 2023; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Version accepted for Astronomy & Astrophysics

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

Anatomy of rocky planets formed by rapid pebble accretion I. How icy pebbles determine the core fraction and FeO contents

Authors: Anders Johansen, Thomas Ronnet, Martin Schiller, Zhengbin Deng, Martin Bizzarro

Abstract: 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 composit… ▽ More 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. △ Less

Submitted 26 January, 2023; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: Version accepted for Astronomy & Astrophysics

Journal ref: A&A 671, A74 (2023)

Abundances of disk and bulge giants from high-resolution optical spectra V. Molybdenum -- the p-process element

Authors: Rebecca Forsberg, Nils Ryde, Henrik Jönsson, R Michael Rich, Anders Johansen

Abstract: 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 det… ▽ More 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. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 12 pages, 8 figures

Growth after the streaming instability: The radial distance dependence of the planetary growth

Authors: Hyerin Jang, Beibei Liu, Anders Johansen

Abstract: 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 inv… ▽ More 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. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 11 pages, 7 figures, Accepted for publication in A&A

Journal ref: A&A 664, A86 (2022)

Natural separation of two primordial planetary reservoirs in an expanding solar protoplanetary disk

Authors: Beibei Liu, Anders Johansen, Michiel Lambrechts, Martin Bizzarro, Troels Haugbølle

Abstract: 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 bl… ▽ More 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. △ Less

Submitted 22 April, 2022; originally announced April 2022.

Comments: Authors' version, 13 pages, 4 figures. Published in Science Advances in the issue of April 22, 2022

Nucleation and growth of iron pebbles explains the formation of iron-rich planets akin to Mercury

Authors: Anders Johansen, Caroline Dorn

Abstract: 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 modera… ▽ More 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. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

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

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

Authors: G. Lesur, B. Ercolano, M. Flock, M. -K. Lin, C. -C. Yang, J. A. Barranco, P. Benitez-Llambay, J. Goodman, A. Johansen, H. Klahr, G. Laibe, W. Lyra, P. Marcus, R. P. Nelson, J. Squire, J. B. Simon, N. Turner, O. M. Umurhan, A. N. Youdin

Abstract: The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t… ▽ More The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

A low accretion efficiency of planetesimals formed at planetary gap edges

Authors: Linn E. J. Eriksson, Thomas Ronnet, Anders Johansen, Ravit Helled, Claudio Valletta, Antoine C. Petit

Abstract: 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 ga… ▽ More 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. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 661, A73 (2022)

Streaming instability of multiple particle species II -- Numerical convergence with increasing particle number

Authors: Noemi Schaffer, Anders Johansen, Michiel Lambrechts

Abstract: 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 pa… ▽ More 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. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 13 pages, 13 figures, 2 tables, accepted to A\&A

Journal ref: A&A 653, A14 (2021)

How the formation of Neptune shapes the Kuiper belt

Authors: Simona Pirani, Anders Johansen, Alexander J. Mustill

Abstract: 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 migratio… ▽ More 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. △ Less

Submitted 25 April, 2021; originally announced April 2021.

Comments: 13 pages, 11 figures, 1 table; accepted for publication in Astronomy & Astrophysics

Hybrid accretion of carbonaceous chondrites by radial transport across the Jupiter barrier

Authors: Elishevah van Kooten, Martin Schiller, Frederic Moynier, Anders Johansen, Troels Haugboelle, Martin Bizzarro

Abstract: 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… ▽ More 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. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: accepted after peer-review

Journal ref: The Astrophysical Journal, 2021

A pebble accretion model for the formation of the terrestrial planets in the Solar System

Authors: Anders Johansen, Thomas Ronnet, Martin Bizzarro, Martin Schiller, Michiel Lambrechts, Åke Nordlund, Helmut Lammer

Abstract: 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… ▽ More 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. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: Science Advances, in press

Formation of Venus, Earth and Mars: Constrained by isotopes

Authors: H. Lammer, R. Brasser, A. Johansen, M. Scherf, M. Leitzinger

Abstract: 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 an… ▽ More 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. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 49 pages, 7 figures. This is a preprint of an article published in Space Science Reviews. The final authenticated version can be found online at: https://doi.org/10.1007/s11214-020-00778-4

Journal ref: Space Sci Rev 217, 7 (2021)

Planetesimals on eccentric orbits erode rapidly

Authors: Lukas Cedenblad, Noemi Schaffer, Anders Johansen, B. Mehlig, Dhrubaditya Mitra

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2021; v1 submitted 29 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astrophysical J

Report number: NORDITA 2021-084

Journal ref: The Astrophysical Journal 921 (2021) 123

The fate of planetesimals formed at planetary gap edges

Authors: L. E. J. Eriksson, T. Ronnet, A. Johansen

Abstract: The presence of rings and gaps in protoplanetary discs are often ascribed to planet-disc interactions, where dust and pebbles are trapped at the edges of planetary induced gas gaps. Recent work has shown that these are likely sites for planetesimal formation via the streaming instability. Given the large amount of planetesimals that potentially form at gap edges, we address the question of their f… ▽ More The presence of rings and gaps in protoplanetary discs are often ascribed to planet-disc interactions, where dust and pebbles are trapped at the edges of planetary induced gas gaps. Recent work has shown that these are likely sites for planetesimal formation via the streaming instability. Given the large amount of planetesimals that potentially form at gap edges, we address the question of their fate and ability to radially transport solids in protoplanetary discs. We perform a series of N-body simulations of planetesimal orbits, taking into account the effect of gas drag and mass loss via ablation. We consider two planetary systems: one akin to the young Solar System, and another one inspired by HL Tau. In both systems, the close proximity to the gap-opening planets results in large orbital excitations, causing the planetesimals to leave their birth locations and spread out across the disc soon after formation. Planetesimals that end up on eccentric orbits interior of 10au experience efficient ablation, and lose all mass before they reach the innermost disc region. In our nominal Solar System simulation with $\dot{M}_0=10^{-7}\, M_{\odot}\, \textrm{yr}^{-1}$ and $α=10^{-2}$, we find that 70% of the initial planetesimal mass has been ablated after 500kyr. The ablation rate in HL Tau is lower, and only 11% of the initial planetesimal mass has been ablated after 1Myr. The ablated material consist of a mixture of solid grains and vaporized ices, where a large fraction of the vaporized ices re-condense to form solid ice. Assuming that the solids grow to pebbles in the disc midplane, this results in a pebble flux of $\sim 10-100\,M_{\oplus}\textrm{Myr}^{-1}$ through the inner disc. Our results demonstrate that scattered planetesimals can carry a significant flux of solids past planetary-induced gaps in young and massive protoplanetary discs. △ Less

Submitted 8 March, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 648, A112 (2021)

A close-encounter method for simulating the dynamics of planetesimals

Authors: Sebastian Lorek, Anders Johansen

Abstract: 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 st… ▽ More 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. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: accepted for publication in A&A

Journal ref: A&A 644, A14 (2020)

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

Authors: A. Brandenburg, A. Johansen, P. A. Bourdin, W. Dobler, W. Lyra, M. Rheinhardt, S. Bingert, N. E. L. Haugen, A. Mee, F. Gent, N. Babkovskaia, C. -C. Yang, T. Heinemann, B. Dintrans, D. Mitra, S. Candelaresi, J. Warnecke, P. J. Käpylä, A. Schreiber, P. Chatterjee, M. J. Käpylä, X. -Y. Li, J. Krüger, J. R. Aarnes, G. R. Sarson , et al. (12 additional authors not shown)

Abstract: The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the… ▽ More The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 7 pages, submitted to the Journal for Open Source Software (JOSS)

Report number: NORDITA-2020-087

Journal ref: Journal of Open Source Software 6, 2807 (2021)

Transport, destruction and growth of pebbles in the gas envelope of a protoplanet

Authors: Anders Johansen, Åke Nordlund

Abstract: 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 sed… ▽ More 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. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: Accepted for publication in The Astrophysical Journal

The path to instability in compact multi-planetary systems

Authors: Antoine C. Petit, Gabriele Pichierri, Melvyn B. Davies, Anders Johansen

Abstract: 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… ▽ More 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 $ε$, 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 $ε^{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. △ Less

Submitted 18 August, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: 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

Journal ref: A&A 641, A176 (2020)

Capture of satellites during planetary encounters A case study of the Neptunian moons Triton and Nereid

Authors: Daohai Li, Anders Johansen, Alexander J. Mustill, Melvyn B. Davies, Apostolos A. Christou

Abstract: 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… ▽ More 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. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: A&A in press

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

Erosion of planetesimals by gas flow

Authors: Noemi Schaffer, Anders Johansen, Lukas Cedenblad, Bernhard Mehling, Dhrubaditya Mitra

Abstract: 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 planetesi… ▽ More 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. △ Less

Submitted 16 May, 2020; originally announced May 2020.

Journal ref: A&A 639, A39 (2020)

Pebble-driven Planet Formation around Very Low-mass Stars and Brown Dwarfs

Authors: Beibei Liu, Michiel Lambrechts, Anders Johansen, Ilaria Pascucci, Thomas Henning

Abstract: 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 in… ▽ More 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\%$). △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 12 pages, 8 figures, accepted in A&A

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

Dust clearing by radial drift in evolving protoplanetary discs

Authors: Johan Appelgren, Michiel Lambrechts, Anders Johansen

Abstract: 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 viscou… ▽ More 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 $μ$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}$). △ Less

Submitted 29 September, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 18 pages, 18 figures, accepted for publication in A&A

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

Gravitoviscous protoplanetary disks with a dust component. III. Evolution of gas, dust, and pebbles

Authors: Vardan G. Elbakyan, Anders Johansen, Michiel Lambrechts, Vitaly Akimkin, Eduard I. Vorobyov

Abstract: 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… ▽ More 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) △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 13 pages, 10 figures

Journal ref: A&A 637, A5 (2020)

On the structure and mass delivery towards circumplanetary discs

Authors: Matthäus Schulik, Anders Johansen, Bertram Bitsch, Elena Lega, Michiel Lambrechts

Abstract: 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 FA… ▽ More 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. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 19 pages, 18 figures, Accepted to A&A

Journal ref: A&A 642, A187 (2020)

Resonance in the K2-19 system is at odds with its high reported eccentricities

Authors: Antoine C. Petit, Erik A. Petigura, Melvyn B. Davies, Anders Johansen

Abstract: K2-19 hosts a planetary system composed of two outer planets, b and c, with size of $7.0\pm 0.2~R_\oplus$ and $4.1\pm0.2~R_\oplus$ , and an inner planet, d, with a radius of $1.11\pm 0.05 R_\oplus$. A recent analysis of Transit-Timing Variations (TTVs) suggested b and c are close to but not in 3:2 mean motion resonance (MMR) because the classical resonant angles circulate. Such an architecture cha… ▽ More K2-19 hosts a planetary system composed of two outer planets, b and c, with size of $7.0\pm 0.2~R_\oplus$ and $4.1\pm0.2~R_\oplus$ , and an inner planet, d, with a radius of $1.11\pm 0.05 R_\oplus$. A recent analysis of Transit-Timing Variations (TTVs) suggested b and c are close to but not in 3:2 mean motion resonance (MMR) because the classical resonant angles circulate. Such an architecture challenges our understanding of planet formation. Indeed, planet migration through the protoplanetary disk should lead to a capture into the MMR. Here, we show that the planets are in fact, locked into the 3:2 resonance despite circulation of the conventional resonant angles and aligned periapses. However, we show that such an orbital configuration cannot be maintained for more than a few hundred million years due to the tidal dissipation experienced by planet d. The tidal dissipation remains efficient because of a secular forcing of the innermost planet eccentricity by planets b and c. While the observations strongly rule out an orbital solution where the three planets are on close to circular orbits, it remains possible that a fourth planet is affecting the TTVs such that the four planet system is consistent with the tidal constraints. △ Less

Submitted 1 July, 2020; v1 submitted 10 March, 2020; originally announced March 2020.

Comments: Accepted in MNRAS, replaced with published version

Evolution of MU69 from a binary planetesimal into contact by Kozai-Lidov oscillations and nebular drag

Authors: Wladimir Lyra, Andrew N. Youdin, Anders Johansen

Abstract: The New Horizons flyby of the cold classical Kuiper Belt object MU69 showed it to be a contact binary. The existence of other contact binaries in the 1-10km range raises the question of how common these bodies are and how they evolved into contact. Here we consider that the lobes of MU69 formed as a binary in the Solar nebula, and calculate its orbital evolution in the presence of gas. We find tha… ▽ More The New Horizons flyby of the cold classical Kuiper Belt object MU69 showed it to be a contact binary. The existence of other contact binaries in the 1-10km range raises the question of how common these bodies are and how they evolved into contact. Here we consider that the lobes of MU69 formed as a binary in the Solar nebula, and calculate its orbital evolution in the presence of gas. We find that the sub-Keplerian wind of the disk brings the drag timescales for 10km bodies under 1 Myr for quadratic-velocity drag. In the Kuiper belt, however, the drag is linear with velocity and the effect of the wind cancels out as the angular momentum gained in half an orbit is lost in the other half: the drag timescales for 10 km bodies remain over 10 Myr. In this situation we find that a combination of drag and Kozai-Lidov oscillations is a promising channel for collapse. We analytically solve the hierarchical three-body problem with drag and implement it into a Kozai plus tidal friction model. The permanent quadrupoles of the lobes make the Kozai oscillations stochastic, and as drag shrinks the semimajor axis it more easily allows the fluctuations to bring the system into contact. Evolution to contact happens very rapidly (within $10^4$ yr) in the pure, double-average quadrupole, Kozai region between $\approx 85-95^\circ$, and within 3 Myr in the drag-assisted region beyond it. The synergy between $J_2$ and gas drag widens the window of contact to 80-100$^\circ$ initial inclination, over a larger range of semimajor axes. As such, the model predicts an initial contact binary fraction of about 10\% for the cold classicals in the Kuiper belt. The speed at contact deviates from the escape velocity only because of the oblateness. For MU69, the contact velocity should be 3.3-4.2 m/s, in line with the observational evidence from the lack of deformation features and estimate of the tensile strength. △ Less

Submitted 26 November, 2020; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 13 pages (main article) + 15 pg appendix; 14 figures. Accepted by Icarus. To appear in Icarus' special issue on Arrokoth. Proof published online on May 22nd, 2020, at https://www.sciencedirect.com/science/article/abs/pii/S001910352030213X?via%3Dihub

On the coexistence of the streaming instability and the vertical shear instability in protoplanetary disks

Authors: Urs Schäfer, Anders Johansen, Robi Banerjee

Abstract: The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction w… ▽ More The streaming instability is a leading candidate mechanism to explain the formation of planetesimals. Yet, the role of this instability in the driving of turbulence in protoplanetary disks, given its fundamental nature as a linear hydrodynamical instability, has so far not been investigated in detail. We study the turbulence that is induced by the streaming instability as well as its interaction with the vertical shear instability. For this purpose, we employ the FLASH Code to conduct two-dimensional axisymmetric global disk simulations spanning radii from $1$ au to $100$ au, including the mutual drag between gas and dust as well as the radial and vertical stellar gravity. If the streaming instability and the vertical shear instability start their growth at the same time, we find the turbulence in the dust mid-plane layer to be primarily driven by the streaming instability. It gives rise to vertical gas motions with a Mach number of up to ${\sim}10^{-2}$. The dust scale height is set in a self-regulatory manner to about $1\%$ of the gas scale height. In contrast, if the vertical shear instability is allowed to saturate before the dust is introduced into our simulations, then it continues to be the main source of the turbulence in the dust layer. The vertical shear instability induces turbulence with a Mach number of ${\sim}10^{-1}$ and thus impedes dust sedimentation. Nonetheless, we find the vertical shear instability and the streaming instability in combination to lead to radial dust concentration in long-lived accumulations which are significantly denser than those formed by the streaming instability alone. Thus, the vertical shear instability may promote planetesimal formation by creating weak overdensities that act as seeds for the streaming instability. △ Less

Submitted 17 February, 2020; originally announced February 2020.

Comments: 20 pages, 14 figures; accepted for publication in A&A; comments welcome

Journal ref: A&A 635, A190 (2020)

Pebble drift and planetesimal formation in protoplanetary discs with embedded planets

Authors: Linn E. J. Eriksson, Anders Johansen, Beibei Liu

Abstract: Nearly-axisymmetric gaps and rings are commonly observed in protoplanetary discs. The leading theory regarding the origin of these patterns is that they are due to dust trapping at the edges of gas gaps induced by the gravitational torques from embedded planets. If the concentration of solids at the gap edges becomes high enough, it could potentially result in planetesimal formation by the streami… ▽ More Nearly-axisymmetric gaps and rings are commonly observed in protoplanetary discs. The leading theory regarding the origin of these patterns is that they are due to dust trapping at the edges of gas gaps induced by the gravitational torques from embedded planets. If the concentration of solids at the gap edges becomes high enough, it could potentially result in planetesimal formation by the streaming instability. We test this hypothesis by performing global 1-D simulations of dust evolution and planetesimal formation in a protoplanetary disc that is perturbed by multiple planets. We explore different combinations of particle sizes, disc parameters, and planetary masses, and find that planetesimals form in all these cases. We also compare the spatial distribution of pebbles from our simulations with protoplanetary disc observations. Planets larger than one pebble isolation mass catch drifting pebbles efficiently at the edge of their gas gaps, and depending on the efficiency of planetesimal formation at the gap edges, the protoplanetary disc transforms within a few 100,000 years to either a transition disc with a large inner hole devoid of dust or to a disc with narrow bright rings. For simulations with planetary masses lower than the pebble isolation mass, the outcome is a disc with a series of weak ring patterns but no strong depletion between the rings. Lowering the pebble size artificially to 100 micrometer-sized "silt", we find that regions between planets get depleted of their pebble mass on a longer time-scale of up to 0.5 million years. These simulations also produce fewer planetesimals than in the nominal model with millimeter-sized particles and always have at least two rings of pebbles still visible after 1 Myr. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: Accepted in A&A

Journal ref: A&A 635, A110 (2020)

Formation of moon systems around giant planets: Capture and ablation of planetesimals as foundation for a pebble accretion scenario

Authors: Thomas Ronnet, Anders Johansen

Abstract: The four major satellites of Jupiter, known as the Galilean moons, and Saturn's most massive satellite, Titan, are believed to have formed in a predominantly gaseous circum-planetary disk, during the last stages of formation of their parent planet. Pebbles from the protoplanetary disk are blocked from flowing into the circumplanetary disk by the positive pressure gradient at the outer edge of the… ▽ More The four major satellites of Jupiter, known as the Galilean moons, and Saturn's most massive satellite, Titan, are believed to have formed in a predominantly gaseous circum-planetary disk, during the last stages of formation of their parent planet. Pebbles from the protoplanetary disk are blocked from flowing into the circumplanetary disk by the positive pressure gradient at the outer edge of the planetary gap, so the gas drag assisted capture of planetesimals should be the main contributor to the delivery of solids onto circum-planetary disks. However, a consistent framework for the subsequent accretion of the moons remains to be built. Here we use numerical integrations to show that most planetesimals being captured within a circum-planetary disk are strongly ablated due to the frictional heating they experience, thus supplying the disk with small dust grains, whereas only a small fraction 'survives' their capture. We then construct a simple model of a circum-planetary disk supplied by ablation, where the flux of solids through the disk is at equilibrium with the ablation supply rate, and investigate the formation of moons in such disks. We show that the growth of satellites is driven mainly by accretion of the pebbles that coagulate from the ablated material. The pebble-accreting protosatellites rapidly migrate inward and pile up in resonant chains at the inner edge of the circum-planetary disk. We propose that dynamical instabilities in these resonant chains are at the origin of the different architectures of Jupiter's and Saturn's moon systems. The assembly of moon systems through pebble accretion can therefore be seen as a down-scaled manifestation of the same process that forms systems of super-Earths and terrestrial-mass planets around solar-type stars and M-dwarfs. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 20 pages, 13 figures; accepted for publication in A&A

Journal ref: A&A 633, A93 (2020)

On the inclinations of the Jupiter Trojans

Authors: Simona Pirani, Anders Johansen, Alexander J. Mustill

Abstract: Jupiter Trojans are are characterized by dark photometric colors, high inclinations and an asymmetry in number of bodies between the two swarms. Different models have been proposed to explain the high inclination of the Trojans and to interpret their relation with the TNOs, but none of them can also satisfactorily explain the asymmetry ratio. Recently it has been found that the asymmetry can arise… ▽ More Jupiter Trojans are are characterized by dark photometric colors, high inclinations and an asymmetry in number of bodies between the two swarms. Different models have been proposed to explain the high inclination of the Trojans and to interpret their relation with the TNOs, but none of them can also satisfactorily explain the asymmetry ratio. Recently it has been found that the asymmetry can arise if Jupiter has migrated inwards by at least a few au during its growth. The asymmetry and the dark colors of the Trojans are natural outcomes of this model, but simulations with massless unperturbed disc particles led to a flat distribution of the Trojan inclinations and a final total mass that was 3-4 orders of magnitude larger than the current one. In our work, we investigate the possible origin of the peculiar inclination distribution of the Trojans in the scenario where Jupiter migrates inwards. We analyze different possibilities: (a) the secular evolution of an initially flat Trojan population, (b) the presence of planetary embryos among the Trojans and (c) capture of the Trojans from a pre-stirred planetesimal population. We find that the secular evolution of the Trojans and Saturn do not affect the Trojan inclinations appreciably, nor is there any significant mass depletion. Embryos embedded in the swarms, in contrast, can stir the Trojan inclinations and can also deplete the swarms efficiently, but it turns out that it is very difficult to get rid of all of the massive bodies. We propose that the disc where Jupiter's core was forming was already excited by the presence of other embryos competing in Jupiter's core's feeding zone. We show that the trapped Trojans preserve their high inclination through the gas phase of the disc and that Saturn's perturbations are more effective on highly inclined Trojans, leading to a lower capture efficiency and to a substantial depletion of the swarms. △ Less

Submitted 3 October, 2019; originally announced October 2019.

Comments: 10 pages, 9 figures. Accepted for publication in A&A

A giant exoplanet orbiting a very low-mass star challenges planet formation models

Authors: J. C. Morales, A. J. Mustill, I. Ribas, M. B. Davies, A. Reiners, F. F. Bauer, D. Kossakowski, E. Herrero, E. Rodríguez, M. J. López-González, C. Rodríguez-López, V. J. S. Béjar, L. González-Cuesta, R. Luque, E. Pallé, M. Perger, D. Baroch, A. Johansen, H. Klahr, C. Mordasini, G. Anglada-Escudé, J. A. Caballero, M. Cortés-Contreras, S. Dreizler, M. Lafarga , et al. (157 additional authors not shown)

Abstract: Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an… ▽ More Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought. △ Less

Submitted 26 September, 2019; originally announced September 2019.

Comments: Manuscript author version. 41 pages, 11 figures

Exploring the conditions for forming cold gas giants through planetesimal accretion

Authors: Anders Johansen, Bertram Bitsch

Abstract: The formation of cold gas giants similar to Jupiter and Saturn in orbit and mass is a great challenge for planetesimal-driven core accretion models because the core growth rates far from the star are low. Here we model the growth and migration of single protoplanets that accrete planetesimals and gas. We integrated the core growth rate using fits in the literature to $N$-body simulations, which pr… ▽ More The formation of cold gas giants similar to Jupiter and Saturn in orbit and mass is a great challenge for planetesimal-driven core accretion models because the core growth rates far from the star are low. Here we model the growth and migration of single protoplanets that accrete planetesimals and gas. We integrated the core growth rate using fits in the literature to $N$-body simulations, which provide the efficiency of accreting the planetesimals that a protoplanet migrates through. We take into account three constraints from the solar system and from protoplanetary discs: (1) the masses of the terrestrial planets and the comet reservoirs in Neptune's scattered disc and the Oort cloud are consistent with a primordial planetesimal population of a few Earth masses per AU, (2) evidence from the asteroid belt and the Kuiper belt indicates that the characteristic planetesimal diameter is 100 km, and (3) observations of protoplanetary discs indicate that the dust is stirred by weak turbulence; this gas turbulence also excites the inclinations of planetesimals. Our nominal model built on these constraints results in maximum protoplanet masses of $0.1$ Earth masses. Ignoring constraint (1) above, we show that even a planetesimal population of 1,000 Earth masses, corresponding to 50 Earth masses per AU, fails to produce cold gas giants (although it successfully forms hot and warm gas giants). We conclude that a massive planetesimal reservoir is in itself insufficient to produce cold gas giants. The formation of cold gas giants by planetesimal accretion additionally requires that planetesimals are small and that the turbulent stirring is very weak, thereby violating all three above constraints. △ Less

Submitted 11 October, 2019; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: Version after language editing by the journal