Showing 1–41 of 41 results for author: Lambrechts, M

Rocky planet formation in compact disks around M dwarfs

Authors: M. Sanchez, N. van der Marel, M. Lambrechts, G. D. Mulders, O. M. Guerra-Alvarado

Abstract: Rocky planets in compact configurations are the most common ones around M dwarfs. Many disks around very low mass stars (between 0.1 and 0.5 M$_\odot$) are rather compact and small (without observable substructures and radius less than 20 au), which favours the idea of an efficient radial drift that could enhance planet formation in compact orbits. We aim to investigate the potential formation pat… ▽ More Rocky planets in compact configurations are the most common ones around M dwarfs. Many disks around very low mass stars (between 0.1 and 0.5 M$_\odot$) are rather compact and small (without observable substructures and radius less than 20 au), which favours the idea of an efficient radial drift that could enhance planet formation in compact orbits. We aim to investigate the potential formation paths of the observed close-in rocky exoplanet population around M dwarfs, assuming that planet formation could take place in compact disks with an efficient dust radial drift. We developed N-body simulations that include a sample of embryos growing by pebble accretion exposed to planet-disk interactions, star-planet tidal interactions and general relativistic corrections. For a star of 0.1 M$_\odot$ we considered different gas disk viscosity and initial embryo distributions. We also explore planet formation by pebble accretion around stars of 0.3 and 0.5 M$_\odot$. Lastly, for each stellar mass, we run simulations that include a sample of embryos growing by planetesimal accretion. Our main result is that the sample of simulated planets that grow by pebble accretion in a gas disk with low viscosity ($α=10^{-4}$) can reproduce the close-in low-mass exoplanet population around M dwarfs in terms of multiplicity, masses and semi-major axis. Furthermore, we found that a gas disk with high viscosity ($α=10^{-3}$) can not reproduce the observed planet masses. Also, we show that planetesimal accretion favours the formation of smaller planets than the ones formed by pebble accretion. Rocky planet formation around M dwarfs can take place in compact and small dust disks driven by an efficient radial drift in a gas disk with low viscosity. This result points towards a new approach in the direction of the disk conditions needed for rocky planet formation around very low mass stars. △ Less

Submitted 9 July, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted for publication in A&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.

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)

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 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)

JWST reveals excess cool water near the snowline in compact disks, consistent with pebble drift

Authors: Andrea Banzatti, Klaus M. Pontoppidan, John Carr, Evan Jellison, Ilaria Pascucci, Joan Najita, Carlos E. Munoz-Romero, Karin I. Oberg, Anusha Kalyaan, Paola Pinilla, Sebastiaan Krijt, Feng Long, Michiel Lambrechts, Giovanni Rosotti, Gregory J. Herczeg, Colette Salyk, Ke Zhang, Edwin Bergin, Nick Ballering, Michael R. Meyer, Simon Bruderer, the JDISCS collaboration

Abstract: Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles dri… ▽ More Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snowline. Previous analyses of IRS water spectra, however, were uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snowline is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 K to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks, compared to the large disks, demonstrating an enhanced cool component with $T \approx$ 170-400 K and equivalent emitting radius $R_{\rm{eq}}\approx$ 1-10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snowline, suggesting that there is indeed a higher inwards mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST. △ Less

Submitted 3 September, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ Letters

The Effect of Dust Evolution and Traps on Inner Disk Water Enrichment

Authors: Anusha Kalyaan, Paola Pinilla, Sebastiaan Krijt, Andrea Banzatti, Giovanni Rosotti, Gijs D. Mulders, Michiel Lambrechts, Feng Long, Gregory J. Herczeg

Abstract: Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance wit… ▽ More Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance within the water snow line region (within a few au). In this study, we aim to understand what effect the presence of gaps in the outer gas disk may have on water vapor enrichment in the inner disk. Building on previous work, we employ a volatile-inclusive multi-Myr disk evolution model that considers an evolving ice-bearing drifting dust population, sensitive to dust-traps, which loses its icy content to sublimation upon reaching the snow line. We find that vapor abundance in the inner disk is strongly affected by fragmentation velocity (v$_{\rm f}$) and turbulence, which control how intense vapor enrichment from pebble delivery is, if present, and how long it may last. Generally, for disks with low to moderate turbulence ($α$ $\le$ 1 $\times$ 10$^{-3}$) and for a range of v$_{\rm f}$, radial location, and gap depth (especially that of the innermost gaps), can significantly alter enrichment. Shallow inner gaps may continuously leak material from beyond it, despite the presence of additional deep outer gaps. We finally find that the for realistic v$_{\rm f}$ ($\le$ 10 m s$^{-1}$), presence of gaps is more important than planetesimal formation beyond the snow line in regulating pebble and volatile delivery into the inner disk. △ Less

Submitted 4 July, 2023; originally announced July 2023.

Comments: 24 Pages, 15 Figures, Accepted to ApJ

Disc population synthesis: Decrease in the solid mass reservoir through pebble drift

Authors: Johan Appelgren, Michiel Lambrechts, Nienke van der Marel

Abstract: Surveys of star-forming regions reveal that the dust mass of protoplanetary discs decreases by several orders of magnitude on a timescale of a few million years. This decrease in the mass budget of solids is likely due to the gas-drag-induced radial drift of mm-sized solids, called pebbles. However, quantifying the evolution of this dust component in young stellar clusters is difficult due to the… ▽ More Surveys of star-forming regions reveal that the dust mass of protoplanetary discs decreases by several orders of magnitude on a timescale of a few million years. This decrease in the mass budget of solids is likely due to the gas-drag-induced radial drift of mm-sized solids, called pebbles. However, quantifying the evolution of this dust component in young stellar clusters is difficult due to the inherent large spread in stellar masses and formation times. Therefore, we aim to model the collective evolution of a cluster to investigate the effectiveness of radial drift in clearing the discs of mm-sized particles. We use a protoplanetary disc model that numerically solves for disc formation, and the viscous evolution and photoevaporative clearing of the gas component, while also including the drift of particles limited in size by fragmentation. We find that discs are born with dust masses between 50 Earth masses and 1000 Earth masses, for stars with, respectively, masses between 0.1 solar masses and 1 solar masses. The majority of this initial dust reservoir is typically lost through drift before photoevaporation opens a gap in the gas disc for models both with and without strong X-ray-driven mass loss rates. We conclude that the decrease in time of the mass locked in fragmentation-limited pebbles is consistent with the evolution of dust masses and ages inferred from nearby star-forming regions when assuming viscous evolution rates corresponding to mean gas disc lifetimes between 3 Myr and 8 Myr. △ Less

Submitted 1 August, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Comments: 16 pages, 11 figures, accepted for publication in A&A. Addressed additional language comments

Journal ref: A&A 673, A139 (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

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

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

A highly settled disk around Oph 163131

Authors: M. Villenave, K. R. Stapelfeldt, G. Duchene, F. Menard, M. Lambrechts, A. Sierra, C. Flores, W. R. F. Dent, S. Wolff, A. Ribas, M. Benisty, N. Cuello, C. Pinte

Abstract: High dust density in the midplane of protoplanetary disks is favorable for efficient grain growth and can allow fast formation of planetesimals and planets, before disks dissipate. Vertical settling and dust trapping in pressure maxima are two mechanisms allowing dust to concentrate in geometrically thin and high density regions. In this work, we aim to study these mechanisms in the highly incline… ▽ More High dust density in the midplane of protoplanetary disks is favorable for efficient grain growth and can allow fast formation of planetesimals and planets, before disks dissipate. Vertical settling and dust trapping in pressure maxima are two mechanisms allowing dust to concentrate in geometrically thin and high density regions. In this work, we aim to study these mechanisms in the highly inclined protoplanetary disk SSTC2D J163131.2-242627 (Oph163131, i~84deg). We present new high angular resolution continuum and 12CO ALMA observations of Oph163131. The gas emission appears significantly more extended in the vertical and radial direction compared to the dust emission, consistent with vertical settling and possibly radial drift. In addition, the new continuum observations reveal two clear rings. The outer ring, located at ~100 au, is well resolved in the observations, which allows us to put stringent constraints on the vertical extent of millimeter dust particles. We model the disk using radiative transfer and find that the scale height of millimeter sized grains is 0.5au or less at 100au from the central star. This value is about one order of magnitude smaller than the scale height of smaller micron-sized dust grains constrained by previous modeling, which implies that efficient settling of the large grains is occurring in the disk. When adopting a parametric dust settling prescription, we find that the observations are consistent with a turbulent viscosity coefficient of about alpha<=10^-5 at 100au. Finally, we find that the thin dust scale height measured in Oph163131 is favorable for planetary growth by pebble accretion: a 10 M_Earth planet may grow within less than 10 Myr, even in orbits exceeding 50au. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ

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

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

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

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

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

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)

Hot Jupiters, cold kinematics: High phase space densities of host stars reflect an age bias

Authors: Alexander J. Mustill, Michiel Lambrechts, Melvyn B. Davies

Abstract: Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution, through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, Hot Jupiters are found… ▽ More Context. The birth environments of planetary systems are thought to influence planet formation and orbital evolution, through external photoevaporation and stellar flybys. Recent work has claimed observational support for this, in the form of a correlation between the properties of planetary systems and the local Galactic phase space density of the host star. In particular, Hot Jupiters are found overwhelmingly around stars in regions of high phase space density, which may reflect a formation environment with high stellar density. Aims. We instead investigate whether the high phase space density may have a galactic kinematic origin: Hot Jupiter hosts may be biased towards being young and therefore kinematically cold, because tidal inspiral leads to the destruction of the planets on Gyr timescales, and the velocity dispersion of stars in the Galaxy increases on similar timescales. Methods. We use 6D positions and kinematics from Gaia for the Hot Jupiter hosts and their neighbours, and construct distributions of the phase space density. We investigate correlations between the stars' local phase space density and peculiar velocity. Results. We find a strong anticorrelation between the phase space density and the host star's peculiar velocity with respect to the Local Standard of Rest. Therefore, most stars in "high-density" regions are kinematically cold, which may be caused by the aforementioned bias towards detecting Hot Jupiters around young stars before the planets' tidal destruction. Conclusions. We do not find evidence in the data for Hot Jupiter hosts preferentially being in phase space overdensities compared to other stars of similar kinematics, nor therefore for their originating in birth environments of high stellar density. △ Less

Submitted 29 November, 2021; v1 submitted 29 March, 2021; originally announced March 2021.

Comments: Accepted to A&A. Code and ancillary files to reproduce these results can be found at https://github.com/AJMustill/HJGalaxy

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

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

Influence of grain growth on the thermal structure of protoplanetary discs

Authors: Sofia Savvidou, Bertram Bitsch, Michiel Lambrechts

Abstract: The thermal structure of a protoplanetary disc is regulated by the opacity that dust grains provide. However, previous works have often considered simplified prescriptions for the dust opacity in hydrodynamical disc simulations, e.g. by considering only a single particle size. In the present work we perform 2D hydrodynamical simulations of protoplanetary discs where the opacity is self-consistentl… ▽ More The thermal structure of a protoplanetary disc is regulated by the opacity that dust grains provide. However, previous works have often considered simplified prescriptions for the dust opacity in hydrodynamical disc simulations, e.g. by considering only a single particle size. In the present work we perform 2D hydrodynamical simulations of protoplanetary discs where the opacity is self-consistently calculated for the dust population, taking into account the particle size, composition and abundance. We first compare simulations using single grain sizes to two different multi-grain size distributions at different levels of turbulence strengths, parameterized through the $α$-viscosity, and different gas surface densities. Assuming a single dust size leads to inaccurate calculations of the thermal structure of discs, because the grain size dominating the opacity increases with orbital radius. Overall the two grain size distributions, one limited by fragmentation only and the other determined from a more complete fragmentation-coagulation equilibrium, give similar results for the thermal structure. We find that both grain size distributions give less steep opacity gradients that result in less steep aspect ratio gradients, in comparison to discs with only micrometer sized dust. Moreover, in the discs with a grain size distribution, the innermost outward migration region is removed and planets embedded is such discs experience lower migration rates. We also investigate the dependency of the water iceline position on the alpha-viscosity, the initial gas surface density at 1 AU and the dust-to-gas ratio and find $r_{ice} \propto α^{0.61} Σ_{g,0}^{0.8} f_{DG}^{0.37}$ independently of the distribution used. The inclusion of the feedback loop between grain growth, opacities and disc thermodynamics allows for more self-consistent simulations of accretion discs and planet formation. △ Less

Submitted 28 May, 2020; originally announced May 2020.

Comments: Accepted by A&A, 27 pages, 19 figures

Journal ref: A&A 640, A63 (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)

Super-Earth masses sculpted by pebble isolation around stars of different masses

Authors: Beibei Liu, Michiel Lambrechts, Anders Johansen, Fan Liu

Abstract: We develop a pebble-driven model to study the formation and evolution of planets around stars in the mass range of 0.08 and 1 solar mass. The growth and migration of a large number of individual protoplanetary embryos are simulated in a population synthesis manner. We test two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplaneta… ▽ More We develop a pebble-driven model to study the formation and evolution of planets around stars in the mass range of 0.08 and 1 solar mass. The growth and migration of a large number of individual protoplanetary embryos are simulated in a population synthesis manner. We test two hypotheses for the birth locations of embryos: at the water ice line or log-uniformly distributed over entire protoplanetary disks. Two types of disks with different turbulent viscous parameters alpha of 1e-3 and 1e-4 are investigated, to shed light on the role of outward migration of protoplanets. The forming planets are compared with the observed exoplanets in terms of masses, semimajor axes, metallicities, and water contents. We find that gas giant planets are likely to form when the characteristic disk sizes are larger, the disk accretion rates are higher, the disks are more metal-rich and/or their stellar hosts are more massive. Our model shows that 1) the characteristic mass of super-Earth is set by the pebble isolation mass. Super-Earth masses increase linearly with the mass of its stellar host, corresponding to one Earth mass around a late M-dwarf star and 20 Earth masses around a solar-mass star. 2) The low-mass planets up to 20 Earth masses can form around stars with a wide range of metallicities, while massive gas giant planets are preferred to grow around metal-rich stars. 3) Super-Earth planets that are mainly composed of silicates, with relatively low water fractions can form from protoplanetary embryos at the water ice line in weakly turbulent disks where outward migration is suppressed. However, if the embryos are formed over a wide range of radial distances, the super-Earths would end up having a distinctive, bimodal composition in water mass. Our model succeeds in quantitatively reproducing several important observed properties of exoplanets and correlations with their stellar hosts. △ Less

Submitted 2 September, 2019; originally announced September 2019.

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

Journal ref: A&A 632, A7 (2019)

Quasi-static contraction during runaway gas accretion onto giant planets

Authors: Michiel Lambrechts, Elena Lega, Richard P. Nelson, Aurélien Crida, Alessandro Morbidelli

Abstract: Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 1… ▽ More Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: Accepted by A&A

Journal ref: A&A 630, A82 (2019)

Oscillatory migration of accreting protoplanets driven by a 3D distortion of the gas flow

Authors: Ondřej Chrenko, Michiel Lambrechts

Abstract: Context. The dynamics of a low-mass protoplanet accreting solids is influenced by the heating torque, which was found to suppress inward migration in protoplanetary disks with constant opacities. Aims. We investigate the differences of the heating torque between disks with constant and temperature-dependent opacities. Methods. Interactions of a super-Earth-sized protoplanet with the gas disk a… ▽ More Context. The dynamics of a low-mass protoplanet accreting solids is influenced by the heating torque, which was found to suppress inward migration in protoplanetary disks with constant opacities. Aims. We investigate the differences of the heating torque between disks with constant and temperature-dependent opacities. Methods. Interactions of a super-Earth-sized protoplanet with the gas disk are explored using 3D radiation hydrodynamic simulations. Results. Accretion heating of the protoplanet creates a hot underdense region in the surrounding gas, leading to misalignment of the local density and pressure gradients. As a result, the 3D gas flow is perturbed and some of the streamlines form a retrograde spiral rising above the protoplanet. In the constant-opacity disk, the perturbed flow reaches a steady state and the underdense gas responsible for the heating torque remains distributed in accordance with previous studies. If the opacity is non-uniform, however, the differences in the disk structure can lead to more vigorous streamline distortion and eventually to a flow instability. The underdense gas develops a one-sided asymmetry which circulates around the protoplanet in a retrograde fashion. The heating torque thus strongly oscillates in time and does not on average counteract inward migration. Conclusions. The torque variations make the radial drift of the protoplanet oscillatory, consisting of short intervals of alternating rapid inward and outward migration. We speculate that transitions between the positive and oscillatory heating torque may occur in specific disk regions susceptible to vertical convection, resulting in the convergent migration of multiple planetary embryos. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: Accepted for publication in A&A, 19 pages, 18 figures

Journal ref: A&A 626, A109 (2019)

Formation of planetary systems by pebble accretion and migration: Hot super-Earth systems from breaking compact resonant chains

Authors: André Izidoro, Bertram Bitsch, Sean N. Raymond, Anders Johansen, Alessandro Morbidelli, Michiel Lambrechts, Seth A. Jacobson

Abstract: At least 30\% of main sequence stars host planets with sizes of between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk--planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-N… ▽ More At least 30\% of main sequence stars host planets with sizes of between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk--planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes ($\lesssim15M_{\oplus}$) and super-massive planetary cores potentially able to become gas giant planets ($\gtrsim15M_{\oplus}$). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the inner edge of the disk forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to $\sim$95\% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the original resonant configurations. Our simulations naturally match observations when they produce a dominant fraction ($\gtrsim95\%$) of unstable systems with a sprinkling ($\lesssim5\%$) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (${\rm N\geq2}$) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet--planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths (abridged). △ Less

Submitted 19 April, 2021; v1 submitted 23 February, 2019; originally announced February 2019.

Comments: Accepted in A&A, version including language editing

Journal ref: A&A 650, A152 (2021)

Formation of planetary systems by pebble accretion and migration: Growth of gas giants

Authors: Bertram Bitsch, André Izidoro, Anders Johansen, Sean N. Raymond, Alessandro Morbidelli, Michiel Lambrechts, Seth A. Jacobson

Abstract: Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nomi… ▽ More Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r<1AU) and form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explore slower migration rates which allow the formation of gas giants from embryos originating from the 5-10AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. Giant planet growth requires a sufficiently-high pebble flux to enable growth to out-compete migration. Even higher pebble fluxes produce systems with multiple gas giants. We show that planetary embryos starting interior to 5AU do not grow into gas giants, even if migration is slow and the pebble flux is large. Instead they grow to the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is independent of the pebble flux. Additionally we show that the long term evolution of our formed planetary systems can produce systems with hot super-Earths and outer gas giants as well as systems of giants on eccentric orbits (abridged). △ Less

Submitted 1 March, 2019; v1 submitted 23 February, 2019; originally announced February 2019.

Comments: 25 pages, accepted by A&A, see companion papers by Izidoro et al. (2019) and Lambrechts et al. (2019)

Journal ref: A&A 623, A88 (2019)

Formation of planetary systems by pebble accretion and migration: How the radial pebble flux determines a terrestrial-planet or super-Earth growth mode

Authors: Michiel Lambrechts, Alessandro Morbidelli, Seth A. Jacobson, Anders Johansen, Bertram Bitsch, Andre Izidoro, Sean N. Raymond

Abstract: Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation… ▽ More Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles. △ Less

Submitted 22 February, 2019; originally announced February 2019.

Comments: Accepted in A&A

Journal ref: A&A 627, A83 (2019)

Observation of aerodynamic instability in the flow of a particle stream in a dilute gas

Authors: Holly L. Capelo, Jan Molacek, Michiel Lambrechts, John Lawson, Anders Johansen, Jurgen Blum, Eberhard Bodenschatz, Haitao Xu

Abstract: Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability… ▽ More Forming macroscopic solid bodies in circumstellar discs requires local dust concentration levels significantly higher than the mean. Interactions of the dust particles with the gas must serve to augment local particle densities, and facilitate growth past barriers in the metre size range. Amongst a number of mechanisms that can amplify the local density of solids, aerodynamic streaming instability (SI) is one of the most promising. This work tests the physical assumptions of models that lead to SI in protoplanetary disks (PPDs). We conduct laboratory experiments in which we track the three-dimensional motion of spherical solid particles fluidized in a low-pressure, laminar, incompressible, gas stream. The particle sizes span the Stokes-Epstein drag regime transition and the overall dust-to-gas mass density ratio is close to unity. Lambrechts et al. (2016) established the similarity of the laboratory flow to a simplified PPD model flow. We obtain experimental results suggesting an instability due to particle-gas interaction: i) there exist variations in particle concentration in the direction of the mean drag forces; ii) the particles have a tendency to 'catch up' to one another when they are in proximity; iii) particle clumping occurs on very small scales, which implies local enhancements above the background dust-to-gas mass density ratio by factors of several tens; v) the presence of these density enhancements occurs for a mean dust-to-gas mass density ratio approaching or greater than 1; v) we find evidence for collective particle drag reduction when the local particle number density becomes high and when the background gas pressure is high so that the drag is in the continuum regime. The experiments presented here are precedent-setting for observing SI under controlled conditions and may lead to a deeper understanding of how it operates in nature. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: Accepted for publication in A&A; abstract abridged for arXiv

Journal ref: A&A 622, A151 (2019)

Dynamics of multiple protoplanets embedded in gas/pebble disks and its dependence on $Σ$ and $ν$ parameters

Authors: Miroslav Brož, Ondřej Chrenko, David Nesvorný, Michiel Lambrechts

Abstract: Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydroco… ▽ More Protoplanets of Super-Earth sizes may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase eccentricities and prevent resonant captures (Chrenko et al. 2017). We study populations of embryos accreting pebbles using Fargo-Thorin 2D hydrocode. We find that embryos in a disk with high surface density ($Σ_0 = 990\,{\rm g}\,{\rm cm}^{-2}$) undergo `unsuccessful' two-body encounters which do not lead to a merger. Only when a 3rd protoplanet arrives to the convergence zone, three-body encounters lead to mergers. For a low-viscosity disk ($ν= 5\times10^{13}\,{\rm cm}^2\,{\rm s}^{-1}$) a massive coorbital is a possible outcome, for which a pebble isolation develops and the coorbital is stabilised. For more massive protoplanets ($5\,M_\oplus$), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times, until the coorbital configuration is disrupted and planets merge. This may be a pathway how to solve the problem that coorbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets ($0.1\,M_\oplus$) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disk due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surroundings. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory. △ Less

Submitted 5 December, 2018; v1 submitted 8 October, 2018; originally announced October 2018.

Comments: 15 pages, 25 figures, Astronomy & Astrophysics, accepted

Journal ref: A&A 620, A157 (2018)

The WFIRST Exoplanet Microlensing Survey

Authors: David P. Bennett, Rachel Akeson, Jay Anderson, Lee Armus, Etienne Bachelet, Vanessa Bailey, Thomas Barclay, Richard Barry, Jean-Phillipe Beaulieu, Andrea Belini, Dominic J. Benford, Aparna Bhattacharya, Padi Boyd, Valerio Bozza, Sebastiano Calchi Novati, Kenneth Carpenter, Arnaud Cassan, David Ciardi, Andrew Cole, Knicole Colon, Christian Coutures, Martin Dominik, Pascal Fouque, Kevin Grady, Tyler Groff , et al. (49 additional authors not shown)

Abstract: The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefr… ▽ More The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap." △ Less

Submitted 26 March, 2018; v1 submitted 22 March, 2018; originally announced March 2018.

Comments: White paper submitted to the National Academy Committee on an Exoplanet Science Strategy; 6 pages (typo fixed)

Pebble isolation mass --- scaling law and implications for the formation of super-Earths and gas giants

Authors: Bertram Bitsch, Alessandro Morbidelli, Anders Johansen, Elena Lega, Michiel Lambrechts, Aurélien Crida

Abstract: The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be… ▽ More The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $α$. We find that small pebbles, coupled to the gas, with Stokes number $τ_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $α$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System. △ Less

Submitted 12 January, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: 16 pages, accepted by A&A, now with language correction

Journal ref: A&A 612, A30 (2018)

Reduced gas accretion on super-Earths and ice giants

Authors: Michiel Lambrechts, Elena Lega

Abstract: A large fraction of giant planets have gaseous envelopes that are limited to about 10 % of their total mass budget. Such planets are present in the Solar System (Uranus, Neptune) and are frequently observed in short periods around other stars (the so-called Super-Earths). In contrast to these observations, theoretical calculations based on the evolution of hydrostatic envelopes argue that such low… ▽ More A large fraction of giant planets have gaseous envelopes that are limited to about 10 % of their total mass budget. Such planets are present in the Solar System (Uranus, Neptune) and are frequently observed in short periods around other stars (the so-called Super-Earths). In contrast to these observations, theoretical calculations based on the evolution of hydrostatic envelopes argue that such low mass envelopes cannot be maintained around cores exceeding five Earth masses. Instead, under nominal disc conditions, these planets would acquire massive envelopes through runaway gas accretion within the lifetime of the protoplanetary disc. In this work, we show that planetary envelopes are not in hydrostatic balance, which slows down envelope growth. A series of 3-dimensional, global, radiative hydrodynamical simulations reveal a steady state gas flow, which enters through the poles and exits in the disc midplane. Gas is pushed through the outer envelope in about 10 orbital timescales. In regions of the disc that are not significantly dust-depleted, envelope accretion onto cores of about five Earth masses can get stalled as the gas flow enters the deep interior. Accreted solids sublimate deep in the convective interior, but small opacity-providing grains are trapped in the flow and do not settle, which further prevents rapid envelope accretion. The transition to runaway gas accretion can however be reached when cores grow larger than typical Super-Earths, beyond 15 Earth masses, and preferably when disc opacities are below kappa=1 cm^2/g. These findings offer an explanation for the typical low-mass envelopes around the cores of Super-Earths. △ Less

Submitted 2 August, 2017; originally announced August 2017.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 606, A146 (2017)

Eccentricity excitation and merging of planetary embryos heated by pebble accretion

Authors: Ondřej Chrenko, Miroslav Brož, Michiel Lambrechts

Abstract: Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish r… ▽ More Context: Planetary embryos can continue to grow by pebble accretion until they become giant planet cores. Simultaneously, these embryos mutually interact and also migrate due to torques arising from the protoplanetary disk. Aims: Our aim is to investigate how pebble accretion alters the orbital evolution of embryos undergoing the Type-I migration. In particular, we study whether they establish resonant chains, whether these chains are prone to instabilities and if giant planet cores form through embryo merging, thus occurring more rapidly than by pebble accretion alone. Methods: For the first time, we perform self-consistent global-scale radiative hydrodynamic simulations of a two-fluid protoplanetary disk consisting of gas and pebbles, the latter being accreted by embedded embryos. Accretion heating, along with other radiative processes, is accounted for to correctly model the Type-I migration. Results: We track the evolution of four super-Earth-like embryos, initially located in a region where the disk structure allows for a convergent migration. Generally, embryo merging is facilitated by rapidly increasing embryo masses and breaks the otherwise oligarchic growth. Moreover, we find that the orbital eccentricity of each embryo is considerably excited ($\simeq0.03$) due to the presence of an asymmetric underdense lobe of gas, a so-called `hot trail', produced by accretion heating of the embryo's vicinity. Eccentric orbits lead the embryos to frequent close encounters and make resonant locking more difficult. Conclusions: Embryo merging typically produces one massive core ($\gtrsim 10\,\mathrm{M_{E}}$) in our simulations, orbiting near $10\,\mathrm{AU}$. Pebble accretion is naturally accompanied by occurrence of eccentric orbits which should be considered in future efforts to explain the structure of exoplanetary systems. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: Accepted for publication in A&A, 25 pages, 17 figures

Journal ref: A&A 606, A114 (2017)

Spontaneous concentrations of solids through two-way drag forces between gas and sedimenting particles

Authors: Michiel Lambrechts, Anders Johansen, Holly L. Capelo, Jürgen Blum, Eberhard Bodenschatz

Abstract: The behaviour of sedimenting particles depends on the dust-to-gas ratio of the fluid. Linear stability analysis shows that solids settling in the Epstein drag regime would remain homogeneously distributed in non-rotating incompressible fluids, even when dust-to-gas ratios reach unity. However, the non-linear evolution has not been probed before. Here, we present numerical calculations indicating t… ▽ More The behaviour of sedimenting particles depends on the dust-to-gas ratio of the fluid. Linear stability analysis shows that solids settling in the Epstein drag regime would remain homogeneously distributed in non-rotating incompressible fluids, even when dust-to-gas ratios reach unity. However, the non-linear evolution has not been probed before. Here, we present numerical calculations indicating that in a particle-dense mixture solids spontaneously mix out of the fluid and form swarms overdense in particles by at least a factor 10. The instability is caused by mass-loaded regions locally breaking the equilibrium background stratification. The driving mechanism depends on non-linear perturbations of the background flow and shares some similarity to the streaming instability in accretion discs. The resulting particle-rich swarms may stimulate particle growth by coagulation. In the context of protoplanetary discs, the instability could be relevant for aiding small particles to settle to the midplane in the outer disc. Inside the gas envelopes of protoplanets, enhanced settling may lead to a reduced dust opacity, which facilitates the contraction of the envelope. We show that the relevant physical set up can be recreated in a laboratory setting. This will allow our numerical calculations to be investigated experimentally in the future. △ Less

Submitted 4 April, 2016; originally announced April 2016.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 591, A133 (2016)

Fossilized condensation lines in the Solar System protoplanetary disk

Authors: A. Morbidelli, B. Bitsch, A. Crida, M. Gounelle, T. Guillot, S. Jacobson, A. Johansen, M. Lambrechts, E. Lega

Abstract: The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation o… ▽ More The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: In press in Icarus

The growth of planets by pebble accretion in evolving protoplanetary discs

Authors: Bertram Bitsch, Michiel Lambrechts, Anders Johansen

Abstract: The formation of planets depends on the underlying protoplanetary disc structure, which influences both the accretion and migration rates of embedded planets. The disc itself evolves on time-scales of several Myr during which both temperature and density profiles change as matter accretes onto the central star. Here we use a detailed model of an evolving disc to determine the growth of planets by… ▽ More The formation of planets depends on the underlying protoplanetary disc structure, which influences both the accretion and migration rates of embedded planets. The disc itself evolves on time-scales of several Myr during which both temperature and density profiles change as matter accretes onto the central star. Here we use a detailed model of an evolving disc to determine the growth of planets by pebble accretion and their migration through the disc. Cores that reach their pebble isolation mass accrete gas to finally form giant planets with extensive gas envelopes, while planets that do not reach pebble isolation mass are stranded as ice giants and ice planets containing only minor amounts of gas in their envelopes. Unlike earlier population synthesis models, our model works without any artificial reductions in migration speed and for protoplanetary discs with gas and dust column densities similar to those inferred from observations. We find that in our nominal disc model the emergence of planetary embryos preferably occurs after approximately 2 Myr in order to not exclusively form gas giants, but also ice giants and smaller planets. The high pebble accretion rates ensure that critical core masses for gas accretion can be reached at all orbital distances. Gas giant planets nevertheless experience significant reduction in semi-major axes by migration. Considering instead planetesimal accretion for planetary growth, we show that formation time-scales are too long to compete with the migration time-scales and the dissipation time of the protoplanetary disc. Altogether, we find that pebble accretion overcomes many of the challenges in the formation of ice and gas giants in evolving protoplanetary discs. △ Less

Submitted 7 August, 2015; v1 submitted 18 July, 2015; originally announced July 2015.

Comments: Accepted by A&A, now with language editing

Journal ref: A&A 609, C2 (2018)

The great dichotomy of the Solar System: small terrestrial embryos and massive giant planet cores

Authors: A. Morbidelli, M. Lambrechts, S. Jacobson, B. Bitsch

Abstract: The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk.… ▽ More The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer. △ Less

Submitted 4 June, 2015; originally announced June 2015.

Comments: In press in Icarus

Journal ref: Icarus 258 (2015) 418-429

The structure of protoplanetary discs around evolving young stars

Authors: Bertram Bitsch, Anders Johansen, Michiel Lambrechts, Alessandro Morbidelli

Abstract: The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying… ▽ More The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying structure of the protoplanetary disc, specifically the profiles of temperature, gas scale height and density. The commonly used disc structure of the Minimum Mass Solar Nebular (MMSN) is a simple power law in all these quantities. However, protoplanetary disc models with both viscous and stellar heating show several bumps and dips in temperature, scale height and density caused by transitions in opacity, which are missing in the MMSN model. These play an important role in the formation of planets, as they can act as sweet spots for the formation of planetesimals via the streaming instability and affect the direction and magnitude of type-I-migration. We present 2D simulations of accretion discs that feature radiative cooling, viscous and stellar heating, and are linked to the observed evolutionary stages of protoplanetary discs and their host stars. These models allow us to identify preferred planetesimal and planet formation regions in the protoplanetary disc as a function of the disc's metallicity, accretion rate and lifetime. We derive simple fitting formulae that feature all structural characteristics of protoplanetary discs during the evolution of several Myr. These fits are straightforward to apply for modelling any growth stage of planets where detailed knowledge of the underlying disc structure is required. △ Less

Submitted 25 December, 2014; v1 submitted 12 November, 2014; originally announced November 2014.

Comments: Accepted by A&A, v3 corrected small typo in the fitting formulae

Journal ref: A&A 575, A28 (2015)

Forming the cores of giant planets from the radial pebble flux in protoplanetary discs

Authors: Michiel Lambrechts, Anders Johansen

Abstract: The formation of planetary cores must proceed rapidly in order for the giant planets to accrete their gaseous envelopes before the dissipation of the protoplanetary gas disc (<3 Myr). In orbits beyond 10 AU, direct accumulation of planetesimals by the cores is too slow. Fragments of planetesimals could be accreted faster, but planetesimals are likely too large for fragmentation to be efficient, an… ▽ More The formation of planetary cores must proceed rapidly in order for the giant planets to accrete their gaseous envelopes before the dissipation of the protoplanetary gas disc (<3 Myr). In orbits beyond 10 AU, direct accumulation of planetesimals by the cores is too slow. Fragments of planetesimals could be accreted faster, but planetesimals are likely too large for fragmentation to be efficient, and resonant trapping poses a further hurdle. Here we instead investigate the accretion of small pebbles (mm-cm sizes) that are the natural outcome of an equilibrium between the growth and radial drift of particles. We construct a simplified analytical model of dust coagulation and pebble drift in the outer disc, between 5 AU and 100 AU, which gives the temporal evolution of the solid surface density and the dominant particle size. These two key quantities determine how core growth proceeds at various orbital distances. We find that pebble surface densities are sufficiently high to achieve the inside-out formation of planetary cores within the disc lifetime. The overall efficiency by which dust gets converted to planets can be high, close to 50 % for planetary architectures similar to the Solar System. Growth by pebble accretion in the outer disc is sufficiently fast to overcome catastrophic Type I migration of the cores. These results require protoplanetary discs with large radial extent (~100 AU) and assume a low number of initial seed embryos. Our findings imply that protoplanetary discs with low disc masses, as expected around low-mass stars (<1 M_sun), or with sub-solar dust-to-gas ratios, do not easily form gas-giant planets (M > 100 M_E), but preferentially form Neptune-mass planets or smaller (< 10 M_E). This is consistent with exoplanet surveys which show that gas giants are relatively uncommon around stars of low mass or low metallicity. △ Less

Submitted 4 April, 2016; v1 submitted 26 August, 2014; originally announced August 2014.

Comments: Version accepted for publication in A&A

Journal ref: A&A 572, A107 (2014)

Separating gas-giant and ice-giant planets by halting pebble accretion

Authors: Michiel Lambrechts, Anders Johansen, Alessandro Morbidelli

Abstract: In the Solar System giant planets come in two flavours: 'gas giants' (Jupiter and Saturn) with massive gas envelopes and 'ice giants' (Uranus and Neptune) with much thinner envelopes around their cores. It is poorly understood how these two classes of planets formed. High solid accretion rates, necessary to form the cores of giant planets within the life-time of protoplanetary discs, heat the enve… ▽ More In the Solar System giant planets come in two flavours: 'gas giants' (Jupiter and Saturn) with massive gas envelopes and 'ice giants' (Uranus and Neptune) with much thinner envelopes around their cores. It is poorly understood how these two classes of planets formed. High solid accretion rates, necessary to form the cores of giant planets within the life-time of protoplanetary discs, heat the envelope and prevent rapid gas contraction onto the core, unless accretion is halted. We find that, in fact, accretion of pebbles (~ cm-sized particles) is self-limiting: when a core becomes massive enough it carves a gap in the pebble disc. This halt in pebble accretion subsequently triggers the rapid collapse of the super-critical gas envelope. As opposed to gas giants, ice giants do not reach this threshold mass and can only bind low-mass envelopes that are highly enriched by water vapour from sublimated icy pebbles. This offers an explanation for the compositional difference between gas giants and ice giants in the Solar System. Furthermore, as opposed to planetesimal-driven accretion scenarios, our model allows core formation and envelope attraction within disc life-times, provided that solids in protoplanetary discs are predominantly in pebbles. Our results imply that the outer regions of planetary systems, where the mass required to halt pebble accretion is large, are dominated by ice giants and that gas-giant exoplanets in wide orbits are enriched by more than 50 Earth masses of solids. △ Less

Submitted 26 August, 2014; originally announced August 2014.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 572, A35 (2014)

Rapid growth of gas-giant cores by pebble accretion

Authors: Michiel Lambrechts, Anders Johansen

Abstract: The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized plane… ▽ More The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time-scale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals. △ Less

Submitted 13 June, 2012; v1 submitted 14 May, 2012; originally announced May 2012.

Comments: Version accepted for publication in A&A