-
Sequential giant planet formation initiated by disc substructure
Authors:
Tommy Chi Ho Lau,
Til Birnstiel,
Joanna Drążkowska,
Sebastian Markus Stammler
Abstract:
Planet formation models are necessary to understand the origins of diverse planetary systems. Circumstellar disc substructures have been proposed as preferred locations of planet formation but a complete formation scenario has not been covered by a single model so far. We aim to study the formation of giant planets facilitated by disc substructure and starting with sub-micron-sized dust. We connec…
▽ More
Planet formation models are necessary to understand the origins of diverse planetary systems. Circumstellar disc substructures have been proposed as preferred locations of planet formation but a complete formation scenario has not been covered by a single model so far. We aim to study the formation of giant planets facilitated by disc substructure and starting with sub-micron-sized dust. We connect dust coagulation and drift, planetesimal formation, $N$-body gravity, pebble accretion, planet migration, planetary gas accretion and gap opening in one consistent modelling framework. We find rapid formation of multiple gas giants from the initial disc substructure. The migration trap near the substructure allows the formation of cold gas giants. A new pressure maximum is created at the outer edge of the planetary gap, which triggers the next generation of planet formation resulting in a compact chain of giant planets. A high planet formation efficiency is achieved as the first gas giants are effective in preventing dust from drifting further inwards, which preserves materials for planet formation. Sequential planet formation is a promising framework to explain the formation of chains of gas and ice giants.
△ Less
Submitted 3 July, 2024; v1 submitted 18 June, 2024;
originally announced June 2024.
-
CAI formation in the early Solar System
Authors:
P. Woitke,
J. Drażkowska,
H. Lammer,
K. Kadam,
P. Marigo
Abstract:
Ca-Al-rich inclusions (CAIs) are the oldest dated solid materials in the solar system, found as light-coloured crystalline ingredients in meteorites. Their formation time is commonly associated with age zero of the Solar System. Yet, the physical and chemical processes that once led to the formation of these sub-millimetre to centimetre-sized mineral particles in the early solar nebula are still a…
▽ More
Ca-Al-rich inclusions (CAIs) are the oldest dated solid materials in the solar system, found as light-coloured crystalline ingredients in meteorites. Their formation time is commonly associated with age zero of the Solar System. Yet, the physical and chemical processes that once led to the formation of these sub-millimetre to centimetre-sized mineral particles in the early solar nebula are still a matter of debate. This paper proposes a pathway to form such inclusions during the earliest phases of disc evolution. We combine 1D viscous disc evolutionary models with 2D radiative transfer, equilibrium condensation, and new dust opacity calculations. We show that the viscous heating associated with the high accretion rates in the earliest evolutionary phases causes the midplane inside of about 0.5 au to heat up to limiting temperatures of about 1500-1700 K, but no further. These high temperatures force all refractory material components of the inherited interstellar dust grains to sublimate - except for a few Al-Ca-Ti oxides such as Al2O3, Ca2Al2SiO7, and CaTiO3. Once the Mg-Fe silicates are gone, the dust becomes more transparent and the heat is more efficiently transported to the disc surface, which prevents any further warm-up. This thermostat mechanism keeps these minerals above their annealing temperature for hundreds of thousands of years, which creates large, pure and crystalline particles. These particles are dragged out by the viscously spreading disc. Beyond about 0.5 au, the silicates re-condense on the Ca-Al-rich particles, adding an amorphous silicate matrix. We estimate that this mechanism to produce CAIs works during the first 50000 years of disc evolution. These particles then continue to move outward and populate the entire disc up to radii of about 50 au, before, eventually, the accretion rate subsides, the disc cools, and the particles start to drift inwards.
△ Less
Submitted 24 April, 2024;
originally announced April 2024.
-
Origin of isotopic diversity among carbonaceous chondrites
Authors:
Jan L. Hellmann,
Jonas M. Schneider,
Elias Wölfer,
Joanna Drążkowska,
Christian A. Jansen,
Timo Hopp,
Christoph Burkhardt,
Thorsten Kleine
Abstract:
Carbonaceous chondrites are some of the most primitive meteorites and derive from planetesimals that formed a few million years after the beginning of the solar system. Here, using new and previously published Cr, Ti, and Te isotopic data, we show that carbonaceous chondrites exhibit correlated isotopic variations that can be accounted for by mixing among three major constituents having distinct i…
▽ More
Carbonaceous chondrites are some of the most primitive meteorites and derive from planetesimals that formed a few million years after the beginning of the solar system. Here, using new and previously published Cr, Ti, and Te isotopic data, we show that carbonaceous chondrites exhibit correlated isotopic variations that can be accounted for by mixing among three major constituents having distinct isotopic compositions, namely refractory inclusions, chondrules, and CI chondrite-like matrix. The abundances of refractory inclusions and chondrules are coupled and systematically decrease with increasing amount of matrix. We propose that these correlated abundance variations reflect trapping of chondrule precursors, including refractory inclusions, in a pressure maximum in the disk, which is likely related to the water ice line and the ultimate formation location of Jupiter. The variable abundance of refractory inclusions/chondrules relative to matrix is the result of their distinct aerodynamical properties resulting in differential delivery rates and their preferential incorporation into chondrite parent bodies during the streaming instability, consistent with the early formation of matrix-poor and the later accretion of matrix-rich carbonaceous chondrites. Our results suggest that chondrules formed locally from isotopically heterogeneous dust aggregates which themselves derive from a wide area of the disk, implying that dust enrichment in a pressure trap was an important step to facilitate the accretion of carbonaceous chondrite parent bodies or, more generally, planetesimals in the outer solar system.
△ Less
Submitted 7 March, 2023;
originally announced March 2023.
-
Leaky dust traps: How fragmentation impacts dust filtering by planets
Authors:
Sebastian Markus Stammler,
Tim Lichtenberg,
Joanna Drążkowska,
Tilman Birnstiel
Abstract:
The nucleosynthetic isotope dichotomy between carbonaceous (CC) and non-carbonaceous (NC) meteorites has been interpreted as evidence for spatial separation and the coexistence of two distinct planet-forming reservoirs for several million years in the solar protoplanetary disk. The rapid formation of Jupiter's core within one million years after the formation of calcium-aluminium-rich inclusions (…
▽ More
The nucleosynthetic isotope dichotomy between carbonaceous (CC) and non-carbonaceous (NC) meteorites has been interpreted as evidence for spatial separation and the coexistence of two distinct planet-forming reservoirs for several million years in the solar protoplanetary disk. The rapid formation of Jupiter's core within one million years after the formation of calcium-aluminium-rich inclusions (CAIs) has been suggested as a potential mechanism for spatial and temporal separation. In this scenario, Jupiter's core would open a gap in the disk and trap inward-drifting dust grains in the pressure bump at the outer edge of the gap, separating the inner and outer disk materials from each other. We performed simulations of dust particles in a protoplanetary disk with a gap opened by an early-formed Jupiter core, including dust growth and fragmentation as well as dust transport, using the dust evolution software DustPy. Our numerical experiments indicate that particles trapped in the outer edge of the gap rapidly fragment and are transported through the gap, contaminating the inner disk with outer disk material on a timescale that is inconsistent with the meteoritic record. This suggests that other processes must have initiated or at least contributed to the isotopic separation between the inner and outer Solar System.
△ Less
Submitted 21 January, 2023; v1 submitted 13 January, 2023;
originally announced January 2023.
-
Rapid Formation of Exoplanetesimals Revealed by White Dwarfs
Authors:
A. Bonsor,
T. Lichtenberg,
J. Drazkowska,
A. M. Buchan
Abstract:
The timing of formation for the first planetesimals determines the mode of planetary accretion and their geophysical and compositional evolution. Astronomical observations of circumstellar discs and Solar System geochronology provide evidence for planetesimal formation during molecular cloud collapse, much earlier than previously estimated. Here, we present distinct observational evidence from whi…
▽ More
The timing of formation for the first planetesimals determines the mode of planetary accretion and their geophysical and compositional evolution. Astronomical observations of circumstellar discs and Solar System geochronology provide evidence for planetesimal formation during molecular cloud collapse, much earlier than previously estimated. Here, we present distinct observational evidence from white dwarf planetary systems for planetesimal formation occurring during the first few hundred thousand years after cloud collapse in exoplanetary systems. A significant fraction of white dwarfs have accreted planetary material rich in iron core or mantle material. In order for the exo-asteroids accreted by white dwarfs to form iron cores, substantial heating is required. By simulating planetesimal evolution and collisional evolution we show that the most likely heat source is short-lived radioactive nuclides such as Al-2 (half life of approximately 0.7 Myr). Core-rich materials in the atmospheres of white dwarfs, therefore, provide independent evidence for rapid planetesimal formation, concurrent with star formation.
△ Less
Submitted 14 November, 2022;
originally announced November 2022.
-
Rapid Formation of Massive Planetary Cores in a Pressure Bump
Authors:
Tommy Chi Ho Lau,
Joanna Drążkowska,
Sebastian M. Stammler,
Tilman Birnstiel,
Cornelis P. Dullemond
Abstract:
Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular s…
▽ More
Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular structures observed in protoplanetary disks, starting with sub-micron-sized dust grains. We connect the models of dust coagulation and drift, planetesimal formation in the streaming instability, gravitational interactions between planetesimals, pebble accretion, and planet migration, into one uniform framework. We find that planetesimals forming early at the massive end of the size distribution grow quickly dominantly by pebble accretion. These few massive bodies grow on the timescales of ~100 000 years and stir the planetesimals formed later preventing the emergence of further planetary cores. Additionally, a migration trap occurs allowing for retention of the growing cores. Pressure bumps are favourable locations for the emergence and rapid growth of planetary cores by pebble accretion as the dust density and grain size are increased and the pebble accretion onset mass is reduced compared to a smooth-disk model.
△ Less
Submitted 21 December, 2022; v1 submitted 8 November, 2022;
originally announced November 2022.
-
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.
-
Why do M dwarfs have more transiting planets?
Authors:
Gijs D. Mulders,
Joanna Drążkowska,
Nienke van der Marel,
Fred J. Ciesla,
Ilaria Pascucci
Abstract:
We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in t…
▽ More
We propose a planet formation scenario to explain the elevated occurrence rates of transiting planets around M dwarfs compared to sun-like stars discovered by Kepler. We use a pebble drift and accretion model to simulate the growth of planet cores inside and outside of the snow line. A smaller pebble size interior to the snow line delays the growth of super-Earths, allowing giant planet cores in the outer disk to form first. When those giant planets reach pebble isolation mass they cut off the flow of pebbles to the inner disk and prevent the formation of close-in super-Earths. We apply this model to stars with masses between 0.1 and 2 solar mass and for a range of initial disk masses. We find that the masses of hot super-Earths and of cold giant planets are anti-correlated. The fraction of our simulations that form hot super-Earths is higher around lower-mass stars and matches the exoplanet occurrence rates from Kepler. The fraction of simulations forming cold giant planets is consistent with the stellar mass dependence from radial velocity surveys. A key testable prediction of the pebble accretion hypothesis is that the occurrence rates of super-Earths should decrease again for M dwarfs near the sub-stellar boundary like Trappist-1.
△ Less
Submitted 6 October, 2021;
originally announced October 2021.
-
The radial structure of planetary bodies formed by the streaming instability
Authors:
Rico G. Visser,
Joanna Drążkowska,
Carsten Dominik
Abstract:
Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global an…
▽ More
Comets and small planetesimals are believed to contain primordial building blocks in the form of millimeter to centimeter sized pebbles. One of the viable growing mechanisms to form these small bodies is through the streaming instability (SI) in which pebbles cluster and gravitationally collapse towards a planetesimal or comet in the presence of gas drag. However, most SI simulations are global and lack the resolution to follow the final collapse stage of a pebble cloud within its Hill radius. We aim to track the collapse of a gravitationally bound pebble cloud subject to mutual collisions and gas drag with the representative particle approach. We determine the radial pebble size distribution of the collapsed core and the impact of mutual pebble collisions on the pebble size distribution. We find that virial equilibrium is never reached during the cloud evolution and that, in general, pebbles with given Stokes number (St) collapse towards an optically thick core in a sequence from aerodynamically largest to aerodynamically smallest. We show that at the location for which the core becomes optically thick, the terminal velocity is well below the fragmentation threshold velocity. While collisional processing is negligible during cloud evolution, the collisions that do occur are sticking. These results support the observations that comets and small planetary bodies are composed of primordial pebbles in the milimeter to centimeter size range
△ Less
Submitted 4 February, 2021; v1 submitted 22 January, 2021;
originally announced January 2021.
-
Bifurcation of planetary building blocks during Solar System formation
Authors:
Tim Lichtenberg,
Joanna Drazkowska,
Maria Schönbächler,
Gregor J. Golabek,
Thomas O. Hands
Abstract:
Geochemical and astronomical evidence demonstrate that planet formation occurred in two spatially and temporally separated reservoirs. The origin of this dichotomy is unknown. We use numerical models to investigate how the evolution of the solar protoplanetary disk influenced the timing of protoplanet formation and their internal evolution. Migration of the water snow line can generate two distinc…
▽ More
Geochemical and astronomical evidence demonstrate that planet formation occurred in two spatially and temporally separated reservoirs. The origin of this dichotomy is unknown. We use numerical models to investigate how the evolution of the solar protoplanetary disk influenced the timing of protoplanet formation and their internal evolution. Migration of the water snow line can generate two distinct bursts of planetesimal formation that sample different source regions. These reservoirs evolve in divergent geophysical modes and develop distinct volatile contents, consistent with constraints from accretion chronology, thermo-chemistry, and the mass divergence of inner and outer Solar System. Our simulations suggest that the compositional fractionation and isotopic dichotomy of the Solar System was initiated by the interplay between disk dynamics, heterogeneous accretion, and internal evolution of forming protoplanets.
△ Less
Submitted 21 January, 2021;
originally announced January 2021.
-
How dust fragmentation may be beneficial to planetary growth by pebble accretion
Authors:
Joanna Drazkowska,
Sebastian M. Stammler,
Til Birnstiel
Abstract:
Pebble accretion is an emerging paradigm for the fast growth of planetary cores. Pebble flux and pebble sizes are the key parameters used in the pebble accretion models. We aim to derive the pebble sizes and fluxes from state-of-the-art dust coagulation models, understand their dependence on disk parameters and the fragmentation threshold velocity, and the impact of those on the planetary growth b…
▽ More
Pebble accretion is an emerging paradigm for the fast growth of planetary cores. Pebble flux and pebble sizes are the key parameters used in the pebble accretion models. We aim to derive the pebble sizes and fluxes from state-of-the-art dust coagulation models, understand their dependence on disk parameters and the fragmentation threshold velocity, and the impact of those on the planetary growth by pebble accretion. We use a one-dimensional dust evolution model including dust growth and fragmentation to calculate realistic pebble sizes and mass flux. We use this information to integrate the growth of planetary embryos placed at various locations in the protoplanetary disk. Pebble flux strongly depends on disk properties, such as its size and turbulence level, as well as on the dust aggregates fragmentation threshold. We find that dust fragmentation may be beneficial to planetary growth in multiple ways. First of all, it prevents the solids from growing to very large sizes, for which the efficiency of pebble accretion drops. What is more, small pebbles are depleted at a slower rate, providing a long-lasting pebble flux. As the full coagulation models are computationally expensive, we provide a simple method of estimating pebble sizes and flux in any protoplanetary disk model without substructure and with any fragmentation threshold velocity.
△ Less
Submitted 5 January, 2021;
originally announced January 2021.
-
Planet-induced Vortices with Dust Coagulation in Protoplanetary Disks
Authors:
Ya-Ping Li,
Hui Li,
Shengtai Li,
Tilman Birnstiel,
Joanna Drazkowska,
Sebastian Stammler
Abstract:
In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the…
▽ More
In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation and drift (to the vortex center) are all quite efficient, producing dust particles ranging from micron to $\sim 1.0\ {\rm cm}$, as well as overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite non-uniform inside the vortex, with the mass weighted average dust size at the vortex center ($\sim 4.0$ mm) being a factor of $\sim10$ larger than other vortex regions. Both large ($\sim$ mm) and small (tens of micron) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. After the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size ($\sim$ mm). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at 1.33 mm band, while only distinct hot spots remain at 7.0 mm.
△ Less
Submitted 13 March, 2020;
originally announced March 2020.
-
Ring Morphology with Dust Coagulation in Protoplanetary Disks
Authors:
JT Laune,
Hui Li,
Shengtai Li,
Ya-Ping Li,
Levi G. Walls,
Tilman Birnstiel,
Joanna Drazkowska,
Sebastian Stammler
Abstract:
Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions,…
▽ More
Tidal interactions between the embedded planets and their surrounding protoplanetary disks are often postulated to produce the observed complex dust substructures, including rings, gaps, and asymmetries. In this Letter, we explore the consequences of dust coagulation on the dust dynamics and ring morphology. Coagulation of dust grains leads to dust size growth which, under typical disk conditions, produces faster radial drifts, potentially threatening the dust ring formation. Utilizing 2D hydrodynamical simulations of protoplanetary disks which include a full treatment of dust coagulation, we find that if the planet does not open a gap quickly enough, the formation of an inner ring is impeded due to dust coagulation and subsequent radial drift. Furthermore, we find that a "buildup" of sub-mm sized grains often appears in the dust emission at the outer edge of the dust disk.
△ Less
Submitted 23 December, 2019;
originally announced December 2019.
-
Including Dust Coagulation in Hydrodynamic Models of Protoplanetary Disks: Dust Evolution in the Vicinity of a Jupiter-mass Planet
Authors:
Joanna Drazkowska,
Shengtai Li,
Til Birnstiel,
Sebastian M. Stammler,
Hui Li
Abstract:
Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2-D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is o…
▽ More
Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2-D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is obtained by including multiple dust fluids in a single grid-based hydrodynamic simulation and solving the Smoluchowski equation for dust coagulation on top of solving for the hydrodynamic evolution. We find that fragmentation of dust aggregates trapped in a pressure bump outside of the planetary gap leads to an enhancement in density of small grains. We compare the results obtained from the full coagulation treatment to the commonly used, fixed dust size approach and to previously applied, less computationally intensive methods for including dust coagulation. We find that the full coagulation results cannot be reproduced using the fixed-size treatment, but some can be mimicked using a relatively simple method for estimating the characteristic dust size in every grid cell.
△ Less
Submitted 23 September, 2019;
originally announced September 2019.
-
The DSHARP Rings: Evidence of Ongoing Planetesimal Formation?
Authors:
Sebastian M. Stammler,
Joanna Drazkowska,
Til Birnstiel,
Hubert Klahr,
Cornelis P. Dullemond,
Sean M. Andrews
Abstract:
Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of…
▽ More
Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process.
We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations.
We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Gas accretion damped by dust back-reaction at the snow line
Authors:
Matías Gárate,
Til Birnstiel,
Joanna Drazkowska,
Sebastian Markus Stammler
Abstract:
Context. The water snowline divides dry and icy solid material in protoplanetary disks, and has been thought to significantly affect planet formation at all stages. If dry particles break up more easily than icy ones, then the snowline causes a traffic jam, because small grains drift inward at lower speeds than larger pebbles. Aims. We aim to evaluate the effect of high dust concentrations around…
▽ More
Context. The water snowline divides dry and icy solid material in protoplanetary disks, and has been thought to significantly affect planet formation at all stages. If dry particles break up more easily than icy ones, then the snowline causes a traffic jam, because small grains drift inward at lower speeds than larger pebbles. Aims. We aim to evaluate the effect of high dust concentrations around the snowline onto the gas dynamics. Methods. Using numerical simulations, we model the global radial evolution of an axisymmetric protoplanetary disk. Our model includes particle growth, evaporation and recondensation of water, and the back-reaction of dust onto the gas, taking into account the vertical distribution of dust particles. Results. We find that the dust back-reaction can stop and even reverse the net flux of gas outside the snowline, decreasing the gas accretion rate onto the star to under $50\%$ of its initial value. At the same time the dust accumulates at the snowline, reaching dust-to-gas ratios of $ε\gtrsim 0.8$, and delivers large amounts of water vapor towards the inner disk, as the icy particles cross the snowline. However, the accumulation of dust at the snowline and the decrease in the gas accretion rate only take place if the global dust-to-gas ratio is high ($\varepsilon_0 \gtrsim 0.03$), if the viscous turbulence is low ($α_ν\lesssim 10^{-3} $), if the disk is large enough ($r_c \gtrsim 100\, \textrm{au}$), and only during the early phases of the disk evolution ($t \lesssim 1\, \textrm{Myr}$). Otherwise the dust back-reaction fails to perturb the gas motion.
△ Less
Submitted 27 February, 2020; v1 submitted 18 June, 2019;
originally announced June 2019.
-
Effects of Ringed Structures and Dust Size Growth on Millimeter Observations of Protoplanetary Disks
Authors:
Ya-Ping Li,
Hui Li,
Luca Ricci,
Shengtai Li,
Tilman Birnstiel,
Andrea Isella,
Megan Ansdell,
Feng Yuan,
Joanna Drazkowska,
Sebastian Stammler
Abstract:
The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We ha…
▽ More
The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.
△ Less
Submitted 3 May, 2019;
originally announced May 2019.
-
Dust evolution and satellitesimal formation in circumplanetary disks
Authors:
Joanna Drazkowska,
Judit Szulágyi
Abstract:
It is believed that satellites of giant planets form in circumplanetary disks. Many of the previous contributions assumed that their formation process proceeds similarly to rocky planet formation, via accretion of the satellite seeds, called satellitesimals. However, the satellitesimal formation itself poses a nontrivial problem as the dust evolution in the circumplanetary disk is heavily impacted…
▽ More
It is believed that satellites of giant planets form in circumplanetary disks. Many of the previous contributions assumed that their formation process proceeds similarly to rocky planet formation, via accretion of the satellite seeds, called satellitesimals. However, the satellitesimal formation itself poses a nontrivial problem as the dust evolution in the circumplanetary disk is heavily impacted by fast radial drift and thus dust growth to satellitesimals is hindered. To address this problem, we connected state-of-the-art hydrodynamical simulations of a circumplanetary disk around a Jupiter-mass planet with dust growth and drift model in a post-processing step. We found that there is an efficient pathway to satellitesimal formation if there is a dust trap forming within the disk. Thanks to the natural existence of an outward gas flow region in the hydrodynamical simulation, a significant dust trap arises at the radial distance of 85~R$_{\rm J}$ from the planet, where the dust-to-gas ratio becomes high enough to trigger streaming instability. The streaming instability leads to the efficient formation of the satellite seeds. Because of the constant infall of material from the circumstellar disk and the very short timescale of dust evolution, the circumplanetary disk acts as a satellitesimal factory, constantly processing the infalling dust to pebbles that gather in the dust trap and undergo the streaming instability.
△ Less
Submitted 11 September, 2018; v1 submitted 7 July, 2018;
originally announced July 2018.
-
A sub-grid model for the growth of dust particles in hydrodynamical simulations of protoplanetary disks
Authors:
Tomas Tamfal,
Joanna Drazkowska,
Lucio Mayer,
Clément Surville
Abstract:
We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust…
▽ More
We present the first 2D hydrodynamical finite volume simulations in which dust is fully coupled with the gas, including its back-reaction onto it, and at the same time the dust size is evolving according to coagulation and fragmentation based on a sub-grid model. The aim of this analysis is to present the differences occurring when dust evolution is included relative to simulations with fixed dust size, with and without an embedded Jupiter-mass planet that triggers gap formation. We use the two-fluid polar Godunov-type code RoSSBi developed by Surville et al. 2016 combined with a new local sub-grid method for dust evolution based on the model by Birnstiel et al. 2012. We find striking differences between simulations with variable and fixed dust sizes. The timescales for dust depletion differ significantly and yield a completely different evolution of the dust surface density. In general sharp features such as pile-ups of dust in the inner disk and near gap edges, when a massive planet is present, become much weaker. This has important implications on the interpretation of observed substructure in disks, suggesting that the presence of a massive planet does not necessarily cause sharp gaps and rings in the dust component. Also, particles with different dust sizes show a different distribution, pointing to the importance of multi-wavelength synthetic observations in order to compare with observations by ALMA and other instruments. We also find that simulations adopting fixed intermediate particle sizes, in the range $10^{-2} - 10^{-1}$ cm, best approximate the surface density evolution seen in simulations with dust evolution.
△ Less
Submitted 14 September, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
Planetesimal formation during protoplanetary disk buildup
Authors:
Joanna Drazkowska,
Cornelis P. Dullemond
Abstract:
Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during…
▽ More
Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during the disk buildup stage, in which gas and dust fall onto the disk from its parent molecular cloud. We couple our earlier planetesimal formation model at the water snow line to a simple model of disk formation and evolution. We find that under most conditions planetesimals only form after the buildup stage when the disk becomes less massive and less hot. However, there are parameters for which planetesimals already form during the disk buildup. This occurs when the viscosity driving the disk evolution is intermediate ($α_v \sim 10^{-3}-10^{-2}$) while the turbulent mixing of the dust is reduced compared to that ($α_t \le 0.03 \cdot α_v$), and with the assumption that water vapor is vertically well-mixed with the gas. Such $α_t \ll α_v$ scenario could be expected for layered accretion, where the gas flow is mostly driven by the active surface layers, while the midplane layers, where most of the dust resides, are quiescent.
△ Less
Submitted 23 February, 2023; v1 submitted 1 March, 2018;
originally announced March 2018.
-
Satellites Form Fast & Late: a Population Synthesis for the Galilean Moons
Authors:
M. Cilibrasi,
J. Szulágyi,
L. Mayer,
J. Drażkowska,
Y. Miguel,
P. Inderbitzi
Abstract:
The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical s…
▽ More
The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical simulations. The disc evolution is taken into account within the population synthesis. The satellitesimals were assumed to grow via streaming instability. We find that the moons form fast, often within $10^4$ years, due to the short orbital timescales in the circumplanetary disc. They form in sequence, and many are lost into the planet due to fast type I migration, polluting Jupiter's envelope with typically 15 Earth-masses of metals. The last generation of moons can form very late in the evolution of the giant planet, when the disc has already lost more than the 99% of its mass. The late circumplanetary disc is cold enough to sustain water ice, hence not surprisingly the 85% of the moon population has icy composition. The distribution of the satellite-masses is peaking slightly above Galilean masses, up until a few Earth-masses, in a regime which is observable with the current instrumentation around Jupiter-analog exoplanets orbiting sufficiently close to their host stars. We also find that systems with Galilean-like masses occur in 20% of the cases and they are more likely when discs have long dispersion timescales and high dust-to-gas ratios.
△ Less
Submitted 23 August, 2018; v1 submitted 18 January, 2018;
originally announced January 2018.
-
Planetesimal formation starts at the snow line
Authors:
Joanna Drazkowska,
Yann Alibert
Abstract:
Planetesimal formation stage represents a major gap in our understanding of the planet formation process. The late-stage planet accretion models typically make arbitrary assumptions about planetesimals and pebbles distribution while the dust evolution models predict that planetesimal formation is only possible at some orbital distances. We want to test the importance of water snow line for trigger…
▽ More
Planetesimal formation stage represents a major gap in our understanding of the planet formation process. The late-stage planet accretion models typically make arbitrary assumptions about planetesimals and pebbles distribution while the dust evolution models predict that planetesimal formation is only possible at some orbital distances. We want to test the importance of water snow line for triggering formation of the first planetesimals during the gas-rich phase of protoplanetary disk, when cores of giant planets have to form. We connect prescriptions for gas disk evolution, dust growth and fragmentation, water ice evaporation and recondensation, as well as transport of both solids and water vapor, and planetesimal formation via streaming instability into a single, one-dimensional model for protoplanetary disk evolution. We find that processes taking place around the snow line facilitate planetesimal formation in two ways. First, due to the change of sticking properties between wet and dry aggregates, there is a "traffic jam" inside of the snow line that slows down the fall of solids onto the star. Second, ice evaporation and outward diffusion of water followed by its recondensation increases the abundance of icy pebbles that trigger planetesimal formation via streaming instability just outside of the snow line. Planetesimal formation is hindered by growth barriers and radial drift and thus requires particular conditions to take place. Snow line is a favorable location where planetesimal formation is possible for a wide range of conditions, but still not in every protoplanetary disk model. This process is particularly promoted in large, cool disks with low intrinsic turbulence and increased initial dust-to-gas ratio.
△ Less
Submitted 29 September, 2017;
originally announced October 2017.
-
Close-in planetesimal formation by pile-up of drifting pebbles
Authors:
Joanna Drazkowska,
Yann Alibert,
Ben Moore
Abstract:
The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and…
▽ More
The consistency of planet formation models suffers from the disconnection between the regime of small and large bodies. This is primarily caused by so-called growth barriers: the direct growth of larger bodies is halted at centimetre-sized objects and particular conditions are required for the formation of larger, gravitationally bound planetesimals. We aim to connect models of dust evolution and planetesimal formation to identify regions of protoplanetary discs that are favourable for the formation of kilometre-sized bodies and the first planetary embryos. We combine semi-analytical models of viscous protoplanetary disc evolution, dust growth and drift including backreaction of the dust particles on the gas, and planetesimal formation via the streaming instability into one numerical code. We investigate how planetesimal formation is affected by the mass of the protoplanetary disc, its initial dust content, and the stickiness of dust aggregates. We find that the dust growth and drift leads to a global redistribution of solids. The pile-up of pebbles in the inner disc provides local conditions where the streaming instability is effective. Planetesimals form in an annulus with its inner edge lying between 0.3 AU and 1 AU and its width ranging from 0.3 AU to 3 AU. The resulting surface density of planetesimals follows a radial profile that is much steeper than the initial disc profile. These results support formation of terrestrial planets in the solar system from a narrow annulus of planetesimals, which reproduces their peculiar mass ratios.
△ Less
Submitted 5 September, 2016; v1 submitted 19 July, 2016;
originally announced July 2016.
-
Can dust coagulation trigger streaming instability?
Authors:
Joanna Drazkowska,
Cornelis P. Dullemond
Abstract:
Streaming instability can be a very efficient way of overcoming growth and drift barriers to planetesimal formation. However, it was shown that strong clumping, which leads to planetesimal formation, requires a considerable number of large grains. State-of-the-art streaming instability models do not take into account realistic size distributions resulting from the collisional evolution of dust. We…
▽ More
Streaming instability can be a very efficient way of overcoming growth and drift barriers to planetesimal formation. However, it was shown that strong clumping, which leads to planetesimal formation, requires a considerable number of large grains. State-of-the-art streaming instability models do not take into account realistic size distributions resulting from the collisional evolution of dust. We investigate whether a sufficient quantity of large aggregates can be produced by sticking and what the interplay of dust coagulation and planetesimal formation is. We develop a semi-analytical prescription of planetesimal formation by streaming instability and implement it in our dust coagulation code based on the Monte Carlo algorithm with the representative particles approach. We find that planetesimal formation by streaming instability may preferentially work outside the snow line, where sticky icy aggregates are present. The efficiency of the process depends strongly on local dust abundance and radial pressure gradient, and requires a super-solar metallicity. If planetesimal formation is possible, the dust coagulation and settling typically need ~100 orbits to produce sufficiently large and settled grains and planetesimal formation lasts another ~1000 orbits. We present a simple analytical model that computes the amount of dust that can be turned into planetesimals given the parameters of the disk model.
△ Less
Submitted 16 October, 2014; v1 submitted 14 October, 2014;
originally announced October 2014.
-
Modeling dust growth in protoplanetary disks: The breakthrough case
Authors:
Joanna Drazkowska,
Fredrik Windmark,
Cornelis P. Dullemond
Abstract:
Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of…
▽ More
Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of the coagulation problems. However, for the sweep-up growth driven by the "lucky" breakthrough mechanism, the methods exhibit very different resolution dependencies. With too few mass bins, the Smoluchowski algorithm tends to overestimate the growth rate and the probability of breakthrough. The Monte Carlo method is less resolution dependent in the growth timescale aspect but it tends to underestimate the breakthrough chance due to its limited dynamic mass range. We discuss the features and drawbacks of both the approaches, which may limit their astrophysical applications.
△ Less
Submitted 3 June, 2014;
originally announced June 2014.
-
Planetesimal formation via sweep-up growth at the inner edge of dead zones
Authors:
J. Drazkowska,
F. Windmark,
C. P. Dullemond
Abstract:
The early stages of planet formation are still not well understood. Coagulation models have revealed numerous obstacles to the dust growth, such as the bouncing, fragmentation and radial drift barriers. We study the interplay between dust coagulation and drift in order to determine the conditions in protoplanetary disk that support the formation of planetesimals. We focus on planetesimal formation…
▽ More
The early stages of planet formation are still not well understood. Coagulation models have revealed numerous obstacles to the dust growth, such as the bouncing, fragmentation and radial drift barriers. We study the interplay between dust coagulation and drift in order to determine the conditions in protoplanetary disk that support the formation of planetesimals. We focus on planetesimal formation via sweep-up and investigate whether it can take place in a realistic protoplanetary disk. We have developed a new numerical model that resolves spatial distribution of dust in the radial and vertical dimension. The model uses representative particles approach to follow the dust evolution in protoplanetary disk. The coagulation and fragmentation of solids is taken into account using Monte Carlo method. A collision model adopting the mass transfer effect, that can occur for different-sized dust aggregate collisions, is implemented. We focus on a protoplanetary disk including a pressure bump caused by a steep decline of turbulent viscosity around the snow line. Our results show that sufficient resolution of the vertical disk structure in dust coagulation codes is necessary to obtain adequately short growth timescales, especially in the case of a low turbulence region. We find that a sharp radial variation of the turbulence strength at the inner edge of dead zone promotes planetesimal formation in several ways. It provides a pressure bump that efficiently prevents the dust from drifting inwards. It also causes a radial variation in the size of aggregates at which growth barriers occur, favoring the growth of large aggregates via sweeping up of small particles. In our model, by employing an ad hoc alpha viscosity change near the snow line, it is possible to grow planetesimals by incremental growth on timescales of approximately 10^5 years.
△ Less
Submitted 12 July, 2013; v1 submitted 14 June, 2013;
originally announced June 2013.
-
Comparison between particle and fluid approximations to dust dynamics
Authors:
Joanna Drążkowska,
Michał Hanasz,
Kacper Kowalik
Abstract:
We present a new particle module of the magnetohydrodynamic (MHD) Piernik code. The original multi-fluid grid code based on the Relaxing Total Variation Diminishing (RTVD) scheme has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. In thi…
▽ More
We present a new particle module of the magnetohydrodynamic (MHD) Piernik code. The original multi-fluid grid code based on the Relaxing Total Variation Diminishing (RTVD) scheme has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. In this poster we introduce the scheme used to solve equations of motion for the particles and present the first results coming from the module. The results of test problems are also compared with the results coming from fluid simulations made with Piernik-MHD code. The comparison shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.
△ Less
Submitted 3 September, 2010;
originally announced September 2010.
-
Particle module of Piernik MHD code
Authors:
Joanna Drążkowska,
Michał Hanasz,
Kacper Kowalik
Abstract:
Piernik is a multi-fluid grid magnetohydrodynamic (MHD) code based on the Relaxing Total Variation Diminishing (RTVD) conservative scheme. The original code has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. The comparison between the t…
▽ More
Piernik is a multi-fluid grid magnetohydrodynamic (MHD) code based on the Relaxing Total Variation Diminishing (RTVD) conservative scheme. The original code has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. The comparison between the test problem results and the results coming from fluid simulations made with Piernik code shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.
△ Less
Submitted 30 August, 2010;
originally announced August 2010.
