-
Surviving the Heat: multi-wavelength analysis of V883 Ori reveals that dust aggregates survive the sublimation of their ice mantles
Authors:
Adrien Houge,
Enrique Macías,
Sebastiaan Krijt
Abstract:
Investigating the response of icy dust aggregates to water ice sublimation is essential for understanding the formation and properties of planetesimals in protoplanetary discs. However, their fate remains unclear, as previous studies suggest aggregates could either survive or completely fall apart to (sub)μm-sized grains. Protoplanetary discs around stars undergoing accretion outbursts represent a…
▽ More
Investigating the response of icy dust aggregates to water ice sublimation is essential for understanding the formation and properties of planetesimals in protoplanetary discs. However, their fate remains unclear, as previous studies suggest aggregates could either survive or completely fall apart to (sub)μm-sized grains. Protoplanetary discs around stars undergoing accretion outbursts represent a unique laboratory to study the ice sublimation process, as the water snowline is pushed outward to regions accessible to current observatories. In this work, we aim to understand the aggregates' response to ice sublimation by focusing on V883 Ori, a system currently undergoing a powerful accretion outburst. We present new analysis of archival high resolution ALMA observations of the protoplanetary disc of V883 Ori at 0.88, 1.3, 2.0, and 3.1 mm, and derive new radial spectral index profiles, which we compare with predictions from one-dimensional dust evolution simulations. In the region of V883 Ori where water ice has sublimated, we find lower spectral indices than previously obtained, indicating the presence of cm-sized particles. Coupled with our dust evolution models, we find that the only way to explain their presence is to assume they formed before the outburst, and survived the sublimation process. The resilience of dust aggregates to such intense events leads us to speculate that it may extend to other environments with more gentle heating, such as pebbles drifting through the water snowline in quiescent protoplanetary discs. In that case, it may alter the formation pathway of dry planetesimals interior to the snowline.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
High-contrast JWST-MIRI spectroscopy of planet-forming disks for the JDISC Survey
Authors:
Klaus M. Pontoppidan,
Colette Salyk,
Andrea Banzatti,
Ke Zhang,
Ilaria Pascucci,
Karin I. Oberg,
Feng Long,
Carlos Munoz-Romero,
John Carr,
Joan Najita,
Geoffrey A. Blake,
Nicole Arulanantham,
Sean Andrews,
Nicholas P. Ballering,
Edwin Bergin,
Jenny Calahan,
Douglas Cobb,
Maria Jose Colmenares,
Annie Dickson-Vandervelde,
Anna Dignan,
Joel Green,
Phoebe Heretz,
Greg Herczeg,
Anusha Kalyaan,
Sebastian Krijt
, et al. (4 additional authors not shown)
Abstract:
The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4…
▽ More
The JWST Disk Infrared Spectral Chemistry Survey (JDISCS) aims to understand the evolution of the chemistry of inner protoplanetary disks using the Mid-InfraRed Instrument (MIRI) on the James Webb Space Telescope (JWST). With a growing sample of >30 disks, the survey implements a custom method to calibrate the MIRI Medium Resolution Spectrometer (MRS) to contrasts of better than 1:300 across its 4.9-28 micron spectral range. This is achieved using observations of Themis-family asteroids as precise empirical reference sources. High spectral contrast enables precise retrievals of physical parameters, searches for rare molecular species and isotopologues, and constraints on the inventories of carbon- and nitrogen-bearing species. JDISCS also offers significant improvements to the MRS wavelength and resolving power calibration. We describe the JDISCS calibrated data and demonstrate its quality using observations of the disk around the solar-mass young star FZ Tau. The FZ Tau MIRI spectrum is dominated by strong emission from warm water vapor. We show that the water and CO line emission originates from the disk surface and traces a range of gas temperatures of ~500-1500 K. We retrieve parameters for the observed CO and H2O lines, and show that they are consistent with a radial distribution represented by two temperature components. A high water abundance of n(H2O)~10^-4 fills the disk surface at least out to the 350 K isotherm at 1.5 au. We search the FZ Tau environs for extended emission detecting a large (radius of ~300 au) ring of emission from H2 gas surrounding FZ Tau, and discuss its origin.
△ Less
Submitted 16 January, 2024; v1 submitted 28 November, 2023;
originally announced November 2023.
-
Water-Rich Disks around Late M-stars Unveiled: Exploring the Remarkable Case of Sz114
Authors:
Chengyan Xie,
Ilaria Pascucci,
Feng Long,
Klaus M. Pontoppidan,
Andrea Banzatti,
Anusha Kalyaan,
Colette Salyk,
Yao Liu,
Joan R. Najita,
Paola Pinilla,
Nicole Arulanantham,
Gregory J. Herczeg,
John Carr,
Edwin A. Bergin,
Nicholas P. Ballering,
Sebastiaan Krijt,
Geoffrey A. Blake,
Ke Zhang,
Karin I. Oberg,
Joel D. Green,
the JDISC collaboration
Abstract:
We present an analysis of the JDISC JWST/MIRI-MRS spectrum of Sz~114, an accreting M5 star surrounded by a large dust disk with a shallow gap at $\sim 39$ au. The spectrum is molecular-rich: we report the detection of water, CO, CO$_2$, HCN, C$_2$H$_2$, and H$_2$. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water…
▽ More
We present an analysis of the JDISC JWST/MIRI-MRS spectrum of Sz~114, an accreting M5 star surrounded by a large dust disk with a shallow gap at $\sim 39$ au. The spectrum is molecular-rich: we report the detection of water, CO, CO$_2$, HCN, C$_2$H$_2$, and H$_2$. The only identified atomic/ionic transition is from [NeII] at 12.81 micron. A distinct feature of this spectrum is the forest of water lines with the 17.22 micron emission surpassing that of most mid-to-late M-star disks by an order of magnitude in flux and aligning instead with disks of earlier-type stars. Moreover, flux ratios of C$_2$H$_2$/H$_2$O and HCN/H$_2$O in Sz~114 also resemble those of earlier-type disks, with a slightly elevated CO$_2$/H$_2$O ratio. While accretional heating can boost all infrared lines, the unusual properties of Sz~114 could be explained by the young age of the source, its formation under unusual initial conditions (a large massive disk), and the presence of dust substructures. The latter delays the inward drift of icy pebbles and help preserve a lower C/O ratio over an extended period. In contrast, mid-to-late M-star disks--which are typically faint, small in size, and likely lack significant substructures--may have more quickly depleted the outer icy reservoir and already evolved out of a water-rich inner disk phase. Our findings underscore the unexpected diversity within mid-infrared spectra of mid-to-late M-star disks, highlighting the need to expand the observational sample for a comprehensive understanding of their variations and thoroughly test pebble drift and planet formation models.
△ Less
Submitted 28 November, 2023; v1 submitted 19 October, 2023;
originally announced October 2023.
-
Chemical footprints of giant planet formation. Role of planet accretion in shaping the C/O ratio of protoplanetary disks
Authors:
Haochang Jiang,
Yu Wang,
Chris W. Ormel,
Sebastiaan Krijt,
Ruobing Dong
Abstract:
Protoplanetary disks, the birthplaces of planets, commonly feature bright rings and dark gaps in both continuum and line emission maps. Accreting planets are interacting with the disk, not only through gravity, but also by changing the local irradiation and elemental abundances, which are essential ingredients for disk chemistry. We propose that giant planet accretion can leave chemical footprints…
▽ More
Protoplanetary disks, the birthplaces of planets, commonly feature bright rings and dark gaps in both continuum and line emission maps. Accreting planets are interacting with the disk, not only through gravity, but also by changing the local irradiation and elemental abundances, which are essential ingredients for disk chemistry. We propose that giant planet accretion can leave chemical footprints in the gas local to the planet, which potentially leads to the spatial coincidence of molecular emissions with the planet in ALMA observation. Through 2D multi-fluid hydrodynamical simulations in Athena++ with built-in sublimation, we simulate the process of an accreting planet locally heating up its vicinity, opening a gas gap in the disk, and creating the conditions for C-photochemistry. An accreting planet located outside the methane snowline can render the surrounding gas hot enough to sublimate the C-rich organics off pebbles before they are accreted by the planet. This locally elevates the disk gas-phase C/O ratio, providing a potential explanation for the C$_2$H line-emission rings observed with ALMA. In particular, our findings provide an explanation for the MWC480 disk, where previous work has identified a statistically significant spatial coincidence of line-emission rings inside a continuum gap. Our findings present a novel view of linking the gas accretion of giant planets and their natal disks through the chemistry signals. This model demonstrates that giant planets can actively shape their forming chemical environment, moving beyond the traditional understanding of the direct mapping of primordial disk chemistry onto planets.
△ Less
Submitted 17 July, 2023;
originally announced July 2023.
-
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.
-
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.
-
Collisional evolution of dust and water ice in protoplanetary discs during and after an accretion outburst
Authors:
Adrien Houge,
Sebastiaan Krijt
Abstract:
Most protoplanetary discs are thought to undergo violent and frequent accretion outbursts, during which the accretion rate and central luminosity are elevated for several decades. This temporarily increases the disc temperature, leading to the sublimation of ice species as snowlines move outwards. In this paper, we investigate how an FUor-type accretion outburst alters the growth and appearance of…
▽ More
Most protoplanetary discs are thought to undergo violent and frequent accretion outbursts, during which the accretion rate and central luminosity are elevated for several decades. This temporarily increases the disc temperature, leading to the sublimation of ice species as snowlines move outwards. In this paper, we investigate how an FUor-type accretion outburst alters the growth and appearance of dust aggregates at different locations in protoplanetary discs. We develop a model based on the Monte Carlo approach to simulate locally the coagulation and fragmentation of icy dust particles and investigate different designs for their structure and response to sublimation. Our main finding is that the evolution of dust grains located between the quiescent and outburst water snowlines is driven by significant changes in composition and porosity. The time required for the dust population to recover from the outburst and return to a coagulation/fragmentation equilibrium depends on the complex interplay of coagulation physics and outburst properties, and can take up to 4500 yr at 5 au. Pebble-sized particles, the building blocks of planetesimals, are either deprecated in water ice or completely destroyed, respectively resulting in drier planetesimals or halting their formation altogether. When accretion outbursts are frequent events, the dust can be far from collisional equilibrium for a significant fraction of time, offering opportunities to track past outbursts in discs at millimetre wavelengths. Our results highlight the importance of including accretion outbursts in models of dust coagulation and planet formation.
△ Less
Submitted 20 March, 2023;
originally announced March 2023.
-
Chemical evolution in planet-forming regions with growing grains
Authors:
Christian Eistrup,
L. Ilsedore Cleeves,
Sebastiaan Krijt
Abstract:
[Abridged] Planets and their atmospheres are built from gas and solid material in protoplanetary disks. This solid material grows from smaller, micron-sized grains to larger sizes in the disks, during the process of planet formation. Our goal is to model the compositional evolution of volatile ices on grains of different sizes, assuming both time-dependent grain growth and several constant grain s…
▽ More
[Abridged] Planets and their atmospheres are built from gas and solid material in protoplanetary disks. This solid material grows from smaller, micron-sized grains to larger sizes in the disks, during the process of planet formation. Our goal is to model the compositional evolution of volatile ices on grains of different sizes, assuming both time-dependent grain growth and several constant grain sizes. The state-of-the-art Walsh chemical kinetics code is utilised for modeling chemical evolution. This code has been upgraded to account for the time-evolving sizes of solids. Chemical evolution is modelled locally at four different radii in a protoplanetary disk midplane for up to 10Myr. The evolution is modelled for five different constant grain sizes, and one model where the grain size changes with time according to a grain growth model appropriate for the disk midplane. Local grain growth, with conservation of total grain mass and the assumption of spherical grains, acts to reduced the total grain-surface area that is available for ice-phase reactions. This reduces these reactions efficiency compared to a chemical scenario with a conventional grain-size choice of 0.1$μ$m. The modelled chemical evolution with grain growth leads to increased abundances of H$_{2}$O ice. For carbon in the inner disk, grain growth leads CO gas to overtake CO$_{2}$ ice as dominant carrier, and in the outer disk, CH$_{4}$ ice to become the dominant carrier. Overall, a constant grain size adopted from a grain evolution model leads to almost identical chemical evolution, when compared with chemical evolution with evolving grain sizes. A constant grain size choice, albeit larger than 0.1$μ$m, may therefore be an appropriate simplification when approximating the impact of grain growth on chemical evolution.
△ Less
Submitted 26 July, 2022;
originally announced July 2022.
-
Chemical Habitability: Supply and Retention of Life's Essential Elements During Planet Formation
Authors:
Sebastiaan Krijt,
Mihkel Kama,
Melissa McClure,
Johanna Teske,
Edwin A. Bergin,
Oliver Shorttle,
Kevin J. Walsh,
Sean N. Raymond
Abstract:
Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimatel…
▽ More
Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds.
△ Less
Submitted 18 March, 2022;
originally announced March 2022.
-
Linking Outer Disk Pebble Dynamics and Gaps to Inner Disk Water Enrichment
Authors:
A. Kalyaan,
P. Pinilla,
S. Krijt,
G. D. Mulders,
A. Banzatti
Abstract:
Millimeter continuum imaging of protoplanetary disks reveals the distribution of solid particles and the presence of substructures (gaps and rings) beyond 5-10 au, while infrared (IR) spectra provide access to abundances of gaseous species at smaller disk radii. Building on recent observational findings of an anti-correlation between the inner disk water luminosity and outer dust disk radius, we a…
▽ More
Millimeter continuum imaging of protoplanetary disks reveals the distribution of solid particles and the presence of substructures (gaps and rings) beyond 5-10 au, while infrared (IR) spectra provide access to abundances of gaseous species at smaller disk radii. Building on recent observational findings of an anti-correlation between the inner disk water luminosity and outer dust disk radius, we aim here at investigating the dynamics of icy solids that drift from the outer disk and sublimate their ice inside the snow line, enriching the water vapor that is observed in the IR. We use a volatile-inclusive disk evolution model to explore a range of conditions (gap location, particle size, disk mass, and alpha-viscosity) under which gaps in the outer disk efficiently block the inward drift of icy solids. We find that inner-disk vapor enrichment is highly sensitive to the location of a disk gap, yielding for each particle size a radial "sweet spot" that reduces the inner-disk vapor enrichment to a minimum. For pebbles of 1-10 mm in size, which carry the most mass, this sweet spot is at 7-15 au, suggesting that inner gaps may have a key role in reducing ice delivery to the inner disk and may not allow the formation of Earths and super-Earths. This highlights the importance of observationally determining the presence and properties of inner gaps in disks. Finally, we argue that the inner water vapor abundance can be used as a proxy for estimating the pebble drift efficiency and mass-flux entering the inner disk.
△ Less
Submitted 6 September, 2021;
originally announced September 2021.
-
If you like C/O variations, you should have put a ring on it
Authors:
Nienke van der Marel,
Arthur Bosman,
Sebastiaan Krijt,
Gijs D. Mulders,
Jennifer B. Bergner
Abstract:
The C/O-ratio as traced with C$_2$H emission in protoplanetary disks is fundamental for constraining the formation mechanisms of exoplanets and our understanding of volatile depletion in disks, but current C$_2$H observations show an apparent bimodal distribution which is not well understood, indicating that the C/O distribution is not described by a simple radial dependence. The transport of icy…
▽ More
The C/O-ratio as traced with C$_2$H emission in protoplanetary disks is fundamental for constraining the formation mechanisms of exoplanets and our understanding of volatile depletion in disks, but current C$_2$H observations show an apparent bimodal distribution which is not well understood, indicating that the C/O distribution is not described by a simple radial dependence. The transport of icy pebbles has been suggested to alter the local elemental abundances in protoplanetary disks, through settling, drift and trapping in pressure bumps resulting in a depletion of volatiles in the surface and an increase of the elemental C/O. We combine all disks with spatially resolved ALMA C$_2$H observations with high-resolution continuum images and constraints on the CO snowline to determine if the C$_2$H emission is indeed related to the location of the icy pebbles. We report a possible correlation between the presence of a significant CO-icy dust reservoir and high C$_2$H emission, which is only found in disks with dust rings outside the CO snowline. In contrast, compact dust disks (without pressure bumps) and warm transition disks (with their dust ring inside the CO snowline) are not detected in C$_2$H, suggesting that such disks may never have contained a significant CO ice reservoir. This correlation provides evidence for the regulation of the C/O profile by the complex interplay of CO snowline and pressure bump locations in the disk. These results demonstrate the importance of including dust transport in chemical disk models, for a proper interpretation of exoplanet atmospheric compositions, and a better understanding of volatile depletion in disks, in particular the use of CO isotopologues to determine gas surface densities.
△ Less
Submitted 17 August, 2021;
originally announced August 2021.
-
System-level fractionation of carbon from disk and planetesimal processing
Authors:
Tim Lichtenberg,
Sebastiaan Krijt
Abstract:
Finding and characterizing extrasolar Earth analogs will rely on interpretation of the planetary system's environmental context. The total budget and fractionation between C-H-O species sensitively affect the climatic and geodynamic state of terrestrial worlds, but their main delivery channels are poorly constrained. We connect numerical models of volatile chemistry and pebble coagulation in the c…
▽ More
Finding and characterizing extrasolar Earth analogs will rely on interpretation of the planetary system's environmental context. The total budget and fractionation between C-H-O species sensitively affect the climatic and geodynamic state of terrestrial worlds, but their main delivery channels are poorly constrained. We connect numerical models of volatile chemistry and pebble coagulation in the circumstellar disk with the internal compositional evolution of planetesimals during the primary accretion phase. Our simulations demonstrate that disk chemistry and degassing from planetesimals operate on comparable timescales and can fractionate the relative abundances of major water and carbon carriers by orders of magnitude. As a result, individual planetary systems with significant planetesimal processing display increased correlation in the volatile budget of planetary building blocks relative to no internal heating. Planetesimal processing in a subset of systems increases the variance of volatile contents across planetary systems. Our simulations thus suggest that exoplanetary atmospheric compositions may provide constraints on $when$ a specific planet formed.
△ Less
Submitted 17 May, 2021; v1 submitted 13 May, 2021;
originally announced May 2021.
-
On the stickiness of CO$_{2}$ and H$_{2}$O ice particles
Authors:
Sota Arakawa,
Sebastiaan Krijt
Abstract:
Laboratory experiments revealed that CO$_{2}$ ice particles stick less efficiently than H$_{2}$O ice particles, and there is an order of magnitude difference in the threshold velocity for sticking. However, the surface energies and elastic moduli of CO$_{2}$ and H$_{2}$O ices are comparable, and the reason why CO$_{2}$ ice particles were poorly sticky compared to H$_{2}$O ice particles was unclear…
▽ More
Laboratory experiments revealed that CO$_{2}$ ice particles stick less efficiently than H$_{2}$O ice particles, and there is an order of magnitude difference in the threshold velocity for sticking. However, the surface energies and elastic moduli of CO$_{2}$ and H$_{2}$O ices are comparable, and the reason why CO$_{2}$ ice particles were poorly sticky compared to H$_{2}$O ice particles was unclear. Here we investigate the effects of viscoelastic dissipation on the threshold velocity for sticking of ice particles using the viscoelastic contact model derived by Krijt et al. We find that the threshold velocity for sticking of CO$_{2}$ ice particles reported in experimental studies is comparable to that predicted for perfectly elastic spheres. In contrast, the threshold velocity for sticking of H$_{2}$O ice particles is an order of magnitude higher than that predicted for perfectly elastic spheres. Therefore, we conclude that the large difference in stickiness between CO$_{2}$ and H$_{2}$O ice particles would mainly originate from the difference in the strength of viscoelastic dissipation, which is controlled by the viscoelastic relaxation time.
△ Less
Submitted 12 February, 2021;
originally announced February 2021.
-
Hints for icy pebble migration feeding an oxygen-rich chemistry in the inner planet-forming region of disks
Authors:
Andrea Banzatti,
Ilaria Pascucci,
Arthur D. Bosman,
Paola Pinilla,
Colette Salyk,
Greg J. Herczeg,
Klaus M. Pontoppidan,
Ivan Vazquez,
Andrew Watkins,
Sebastiaan Krijt,
Nathan Hendler,
Feng Long
Abstract:
We present a synergic study of protoplanetary disks to investigate links between inner disk gas molecules and the large-scale migration of solid pebbles. The sample includes 63 disks where two types of measurements are available: i) spatially-resolved disk images revealing the radial distribution of disk pebbles (mm-cm dust grains), from millimeter observations with ALMA or the SMA, and ii) infrar…
▽ More
We present a synergic study of protoplanetary disks to investigate links between inner disk gas molecules and the large-scale migration of solid pebbles. The sample includes 63 disks where two types of measurements are available: i) spatially-resolved disk images revealing the radial distribution of disk pebbles (mm-cm dust grains), from millimeter observations with ALMA or the SMA, and ii) infrared molecular emission spectra as observed with Spitzer. The line flux ratios of H2O with HCN, C2H2, and CO2 all anti-correlate with the dust disk radius R$_{dust}$, expanding previous results found by Najita et al. (2013) for HCN/H2O and the dust disk mass. By normalization with the dependence on accretion luminosity common to all molecules, only the H2O luminosity maintains a detectable anti-correlation with disk radius, suggesting that the strongest underlying relation is between H2O and R$_{dust}$. If R$_{dust}$ is set by large-scale pebble drift, and if molecular luminosities trace the elemental budgets of inner disk warm gas, these results can be naturally explained with scenarios where the inner disk chemistry is fed by sublimation of oxygen-rich icy pebbles migrating inward from the outer disk. Anti-correlations are also detected between all molecular luminosities and the infrared index n$_{13-30}$, which is sensitive to the presence and size of an inner disk dust cavity. Overall, these relations suggest a physical interconnection between dust and gas evolution both locally and across disk scales. We discuss fundamental predictions to test this interpretation and study the interplay between pebble drift, inner disk depletion, and the chemistry of planet-forming material.
△ Less
Submitted 28 September, 2020;
originally announced September 2020.
-
CO Depletion in Protoplanetary Disks: A Unified Picture Combining Physical Sequestration and Chemical Processing
Authors:
Sebastiaan Krijt,
Arthur D. Bosman,
Ke Zhang,
Kamber R. Schwarz,
Fred J. Ciesla,
Edwin A. Bergin
Abstract:
The gas-phase CO abundance (relative to hydrogen) in protoplanetary disks decreases by up to 2 orders of magnitude from its ISM value ${\sim}10^{-4}$, even after accounting for freeze-out and photo-dissociation. Previous studies have shown that while local chemical processing of CO and the sequestration of CO ice on solids in the midplane can both contribute, neither of these processes appears cap…
▽ More
The gas-phase CO abundance (relative to hydrogen) in protoplanetary disks decreases by up to 2 orders of magnitude from its ISM value ${\sim}10^{-4}$, even after accounting for freeze-out and photo-dissociation. Previous studies have shown that while local chemical processing of CO and the sequestration of CO ice on solids in the midplane can both contribute, neither of these processes appears capable of consistently reaching the observed depletion factors on the relevant timescale of $1{-}3\mathrm{~Myr}$. In this study, we model these processes simultaneously by including a compact chemical network (centered on carbon and oxygen) to 2D ($r+z$) simulations of the outer ($r>20\mathrm{~au}$) disk regions that include turbulent diffusion, pebble formation, and pebble dynamics. In general, we find that the CO/H$_2$ abundance is a complex function of time and location. Focusing on CO in the warm molecular layer, we find that only the most complete model (with chemistry and pebble evolution included) can reach depletion factors consistent with observations. In the absence of pressure traps, highly-efficient planetesimal formation, or high cosmic ray ionization rates, this model also predicts a resurgence of CO vapor interior to the CO snowline. We show the impact of physical and chemical processes on the elemental (C/O) and (C/H) ratios (in the gas and ice phases), discuss the use of CO as a disk mass tracer, and, finally, connect our predicted pebble ice compositions to those of pristine planetesimals as found in the Cold Classical Kuiper Belt and debris disks.
△ Less
Submitted 18 July, 2020;
originally announced July 2020.
-
Color, Composition, and Thermal Environment of Kuiper Belt Object (486958) Arrokoth
Authors:
W. M. Grundy,
M. K. Bird,
D. T. Britt,
J. C. Cook,
D. P. Cruikshank,
C. J. A. Howett,
S. Krijt,
I. R. Linscott,
C. B. Olkin,
A. H. Parker,
S. Protopapa,
M. Ruaud,
O. M. Umurhan,
L. A. Young,
C. M. Dalle Ore,
J. J. Kavelaars,
J. T. Keane,
Y. J. Pendleton,
S. B. Porter,
F. Scipioni,
J. R. Spencer,
S. A. Stern,
A. J. Verbiscer,
H. A. Weaver,
R. P. Binzel
, et al. (24 additional authors not shown)
Abstract:
The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrog…
▽ More
The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.
△ Less
Submitted 16 February, 2020;
originally announced February 2020.
-
Tracking Dust Grains During Transport and Growth in Protoplanetary Disks
Authors:
William Misener,
Sebastiaan Krijt,
Fred J. Ciesla
Abstract:
Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribu…
▽ More
Protoplanetary disks are dynamic objects, within which dust grains and gas are expected to be redistributed over large distances. Evidence for this redistribution is seen both in other protoplanetary disks and in our own Solar System, with high-temperature materials thought to originate close to the central star found in the cold, outer regions of the disks. While models have shown this redistribution is possible through a variety of mechanisms, these models have generally ignored the possible growth of solids via grain-grain collisions that would occur during transit. Here we investigate the interplay of coagulation and radial and vertical transport of solids in protoplanetary disks, considering cases where growth is limited by bouncing or by fragmentation. We find that in all cases, growth effectively limits the ability for materials to be carried outward or preserved at large distances from the star. This is due to solids being incorporated into large aggregates which drift inwards rapidly under the effects of gas drag. We discuss the implications for mixing in protoplanetary disks, and how the preservation of high temperature materials in outer disks may require structures or outward flow patterns to avoid them being lost via radial drift.
△ Less
Submitted 1 October, 2019;
originally announced October 2019.
-
Systematic Variations of CO Gas Abundance with Radius in Gas-rich Protoplanetary Disks
Authors:
Ke Zhang,
Edwin A. Bergin,
Kamber R. Schwarz,
Sebastiaan Krijt,
Fred Ciesla
Abstract:
CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks a…
▽ More
CO is the most widely used gas tracer of protoplanetary disks. Its abundance is usually assumed to be an interstellar ratio throughout the warm molecular layer of the disk. But recent observations of low CO gas abundance in many protoplanetary disks challenge our understanding of physical and chemical evolutions in disks. Here we investigate the CO abundance structures in four well-studied disks and compare their structures with predictions of chemical processing of CO and transport of CO ice-coated dust grains in disks. We use spatially resolved CO isotopologue line observations and detailed thermo-chemical models to derive CO abundance structures. We find that the CO abundance varies with radius by an order of magnitude in these disks. We show that although chemical processes can efficiently reduce the total column of CO gas within 1 Myr under an ISM level of cosmic-ray ionization rate, the depletion mostly occurs at the deep region of a disk. Without sufficient vertical mixing, the surface layer is not depleted enough to reproduce weak CO emissions observed. The radial profiles of CO depletion in three disks are qualitatively consistent with predictions of pebble formation, settling, and drifting in disks. But the dust evolution alone cannot fully explain the high depletion observed in some disks. These results suggest that dust evolution may play a significant role in transporting volatile materials and a coupled chemical-dynamical study is necessary to understand what raw materials are available for planet formation at different distances from the central star.
△ Less
Submitted 8 August, 2019;
originally announced August 2019.
-
Dust Transport and Processing in Centrifugally Driven Protoplanetary Disk Winds
Authors:
Steven Giacalone,
Seth Teitler,
Arieh Königl,
Sebastiaan Krijt,
Fred J. Ciesla
Abstract:
There is evidence that protoplanetary disks--including the protosolar one--contain crystalline dust grains on spatial scales where the dust temperature is lower than the threshold value for their formation through thermal annealing of amorphous interstellar silicates. We interpret these observations in terms of an extended, magnetocentrifugally driven disk wind that transports grains from the inne…
▽ More
There is evidence that protoplanetary disks--including the protosolar one--contain crystalline dust grains on spatial scales where the dust temperature is lower than the threshold value for their formation through thermal annealing of amorphous interstellar silicates. We interpret these observations in terms of an extended, magnetocentrifugally driven disk wind that transports grains from the inner disk--where they are thermally processed by the stellar radiation after being uplifted from the disk surfaces--to the outer disk regions. For any disk radius $r$ there is a maximum grain size $a_\mathrm{max}(r)$ that can be uplifted from that location: grains of size $a$$\ll$$a_\mathrm{max}$ are carried away by the wind, whereas those with $a$$\lesssim$$a_\mathrm{max}$ reenter the disk at larger radii. A significant portion of the reentering grains converge to--and subsequently accumulate in--a narrow region just beyond $r_\mathrm{max}(a)$, the maximum radius from which grains of size $a$ can be uplifted. We show that this model can account for the inferred crystallinity fractions in classical T Tauri and Herbig Ae disks and for their indicated near constancy after being established early in the disk evolution. It is also consistent with the reported radial gradients in the mean grain size, crystallinity, and crystal composition. In addition, this model yields the properties of the grains that remain embedded in the outflows from protoplanetary disks and naturally explains the inferred persistence of small grains in the surface layers of these disks.
△ Less
Submitted 10 July, 2019;
originally announced July 2019.
-
Planet formation: The case for large efforts on the computational side
Authors:
Wladimir Lyra,
Thomas Haworth,
Bertram Bitsch,
Simon Casassus,
Nicolás Cuello,
Thayne Currie,
Andras Gáspár,
Hannah Jang-Condell,
Hubert Klahr,
Nathan Leigh,
Giuseppe Lodato,
Mordecai-Mark Mac Low,
Sarah Maddison,
George Mamatsashvili,
Colin McNally,
Andrea Isella,
Sebastián Pérez,
Luca Ricci,
Debanjan Sengupta,
Dimitris Stamatellos,
Judit Szulágyi,
Richard Teague,
Neal Turner,
Orkan Umurhan,
Jacob White
, et al. (32 additional authors not shown)
Abstract:
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been a…
▽ More
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
The Efficiency of Noble Gas Trapping in Astrophysical Environments
Authors:
Fred Ciesla,
Sebastiaan Krijt,
Reika Yokochi,
Scott Sandford
Abstract:
Amorphous ice has long been invoked as a means for trapping extreme volatiles into solids, explaining the abundances of these species in comets and planetary atmospheres. Experiments have shown that such trapping is possible and have been used to estimate the abundances of each species in primitive ices after they formed. However, these experiments have been carried out at deposition rates which e…
▽ More
Amorphous ice has long been invoked as a means for trapping extreme volatiles into solids, explaining the abundances of these species in comets and planetary atmospheres. Experiments have shown that such trapping is possible and have been used to estimate the abundances of each species in primitive ices after they formed. However, these experiments have been carried out at deposition rates which exceed those expected in a molecular cloud or solar nebula by many orders of magnitude. Here we develop a numerical model which reproduces the experimental results and apply it to those conditions expected in molecular clouds and protoplanetary disks. We find that two regimes of ice trapping exist: `burial trapping' where the ratio of trapped species to water in the ice reflects that same ratio in the gas and `equilibrium trapping' where the ratio in the ice depends only on the partial pressure of the trapped species in the gas. The boundary between these two regimes is set by both the temperature and rate of ice deposition. Such effects must be accounted for when determining the source of trapped volatiles during planet formation.
△ Less
Submitted 14 November, 2018;
originally announced November 2018.
-
A Lagrangian Model for Dust Evolution in Protoplanetary Disks: Formation of Wet and Dry Planetesimals at Different Stellar Masses
Authors:
Djoeke Schoonenberg,
Chris W. Ormel,
Sebastiaan Krijt
Abstract:
We introduce a new Lagrangian smooth-particle method to model the growth and drift of pebbles in protoplanetary disks. The Lagrangian nature of the model makes it especially suited to follow characteristics of individual (groups of) particles, such as their composition. In this work we focus on the water content of solid particles. Planetesimal formation via streaming instability is taken into acc…
▽ More
We introduce a new Lagrangian smooth-particle method to model the growth and drift of pebbles in protoplanetary disks. The Lagrangian nature of the model makes it especially suited to follow characteristics of individual (groups of) particles, such as their composition. In this work we focus on the water content of solid particles. Planetesimal formation via streaming instability is taken into account, partly based on previous results on streaming instability outside the water snowline that were presented in Schoonenberg & Ormel (2017). We validate our model by reproducing earlier results from the literature and apply our model to steady-state viscous gas disks (with constant gas accretion rate) around stars with different masses. We also present various other models where we explore the effects of pebble accretion, the fragmentation velocity threshold, the global metallicity of the disk, and a time-dependent gas accretion rate. We find that planetesimals preferentially form in a local annulus outside the water snowline, at early times in the lifetime of the disk ($\lesssim$$10^{5} \: \rm{yr}$), when the pebble mass fluxes are high enough to trigger the streaming instability. During this first phase in the planet formation process, the snowline location hardly changes due to slow viscous evolution, and we conclude that assuming a constant gas accretion rate is justified in this first stage. The efficiency of converting the solids reservoir of the disk to planetesimals depends on the location of the water snowline. Cooler disks with a closer-in water snowline are more efficient at producing planetesimals than hotter disks where the water snowline is located further away from the star. Therefore, low-mass stars tend to form planetesimals more efficiently, but any correlation may be overshadowed by the spread in disk properties.
△ Less
Submitted 4 October, 2018;
originally announced October 2018.
-
Transport of CO in Protoplanetary Disks: Consequences of Pebble Formation, Settling, and Radial Drift
Authors:
Sebastiaan Krijt,
Kamber R. Schwarz,
Edwin A. Bergin,
Fred J. Ciesla
Abstract:
Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline.…
▽ More
Current models of (exo)planet formation often rely on a large influx of so-called `pebbles' from the outer disk into the planet formation region. In this paper, we investigate how the formation of pebbles in the cold outer regions of protoplanetary disks and their subsequent migration to the inner disk can alter the gas-phase CO distribution both interior and exterior to the midplane CO snowline. By simulating the resulting CO abundances in the midplane as well as the warm surface layer, we identify observable signatures of large-scale pebble formation and migration that can be used as `smoking guns' for these important processes. Specifically, we find that after $1\mathrm{~Myr}$, the formation and settling of icy pebbles results in the removal of up to $80\%$ of the CO vapor in the warm ($T>22\mathrm{~K}$) disk layers outside the CO snowline, while the radial migration of pebbles results in the generation of a plume of CO vapor interior the snowline, increasing the CO abundance by a factor ${\sim}2{-}6$ depending on the strength of the turbulence and the sizes of the individual pebbles. The absence of this plume of CO vapor in young nearby disks could indicate efficient conversion of CO into a more refractory species, or a reduction in the radial mass flux of pebbles by, for example, disk inhomogeneities or early planetesimal formation.
△ Less
Submitted 6 August, 2018;
originally announced August 2018.
-
The physics of protoplanetesimal dust agglomerates. X. High-velocity collisions between small and large dust agglomerates as growth barrier
Authors:
Rainer Schräpler,
Jürgen Blum,
Sebastiaan Krijt,
Jan-Hendrik Raabe
Abstract:
In a protoplanetary disk, dust aggregates in the $μ$m to mm size range possess mean collision velocities of 10 to 60 ms$^{-1}$ with respect to dm- to m-size bodies. We performed laboratory collision experiments to explore this parameter regime and found a size- and velocity-dependent threshold between erosion and growth. By using a local Monte Carlo coagulation calculation and complementary a simp…
▽ More
In a protoplanetary disk, dust aggregates in the $μ$m to mm size range possess mean collision velocities of 10 to 60 ms$^{-1}$ with respect to dm- to m-size bodies. We performed laboratory collision experiments to explore this parameter regime and found a size- and velocity-dependent threshold between erosion and growth. By using a local Monte Carlo coagulation calculation and complementary a simple semi-analytical timescale approach, we show that erosion considerably limits particle growth in protoplanetary disks and leads to a steady-state dust-size distribution from $μ$m to dm sized particles.
△ Less
Submitted 16 January, 2018; v1 submitted 12 December, 2017;
originally announced December 2017.
-
Fast litho-panspermia in the habitable zone of the TRAPPIST-1 system
Authors:
Sebastiaan Krijt,
Timothy J. Bowling,
Richard J. Lyons,
Fred J. Ciesla
Abstract:
With several short-period, Earth-mass planets in the habitable zone, the TRAPPIST-1 system potentially allows litho-panspermia to take place on very short timescales. We investigate the efficiency and speed of inter-planetary material transfer resulting from impacts onto the habitable zone planets. By simulating trajectories of impact ejecta from their moment of ejection until (re-)accretion, we f…
▽ More
With several short-period, Earth-mass planets in the habitable zone, the TRAPPIST-1 system potentially allows litho-panspermia to take place on very short timescales. We investigate the efficiency and speed of inter-planetary material transfer resulting from impacts onto the habitable zone planets. By simulating trajectories of impact ejecta from their moment of ejection until (re-)accretion, we find that transport between the habitable zone planets is fastest for ejection velocities around and just above planetary escape velocity. At these ejection velocities, ${\sim}10\%$ of the ejected material reaches another habitable zone planet within $10^2\mathrm{~yr}$, indicating litho-panspermia can be 4 to 5 orders of magnitude faster in TRAPPIST-1 than in the Solar System.
△ Less
Submitted 5 April, 2017;
originally announced April 2017.
-
Tracing water vapor and ice during dust growth
Authors:
Sebastiaan Krijt,
Fred J. Ciesla,
Edwin A. Bergin
Abstract:
The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust gro…
▽ More
The processes that govern the evolution of dust and water (in the form of vapor or ice) in protoplanetary disks are intimately connected. We have developed a model that simulates dust coagulation, dust dynamics (settling, turbulent mixing), vapor diffusion, and condensation/sublimation of volatiles onto grains in a vertical column of a protoplanetary disk. We employ the model to study how dust growth and dynamics influence the vertical distribution of water vapor and water ice in the region just outside the radial snowline. Our main finding is that coagulation (boosted by the enhanced stickiness of icy grains) and the ensuing vertical settling of solids results in water vapor being depleted, but not totally removed, from the region above the snowline on a timescale commensurate with the vertical turbulent mixing timescale. Depending on the strength of the turbulence and the temperature, the depletion can reach factors of up to ${\sim}50$ in the disk atmosphere. In our isothermal column, this vapor depletion results in the vertical snowline moving closer to the midplane (by up to 2 gas scale heights) and the gas-phase C/O ratio above the vertical snowline increasing. Our findings illustrate the importance of dynamical effects and the need for understanding coevolutionary dynamics of gas and solids in planet-forming environments.
△ Less
Submitted 20 October, 2016;
originally announced October 2016.
-
A Steeper than Linear Disk Mass-Stellar Mass Scaling Relation
Authors:
I. Pascucci,
L. Testi,
G. J. Herczeg,
F. Long,
C. F. Manara,
N. Hendler,
G. D. Mulders,
S. Krijt,
F. Ciesla,
Th. Henning,
S. Mohanty,
E. Drabek-Maunder,
D. Apai,
L. Szucs,
G. Sacco,
J. Olofsson
Abstract:
The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and ide…
▽ More
The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.
△ Less
Submitted 11 August, 2016;
originally announced August 2016.
-
Dust diffusion and settling in the presence of collisions: Trapping (sub)micron grains in the midplane
Authors:
Sebastiaan Krijt,
Fred J. Ciesla
Abstract:
In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with th…
▽ More
In protoplanetary disks, the distribution and abundance of small (sub)micron grains are important for a range of physical and chemical processes. For example, they dominate the optical depth at short wavelengths and their surfaces are the sites of many important chemical reactions such as the formation of water. Based on their aerodynamical properties (i.e., their strong dynamical coupling with the surrounding gas) it is often assumed that these small grains are well-mixed with the gas. Our goal is to study the vertical (re)distribution of grains taking into account settling, turbulent diffusion, as well as collisions with other dust grains. Assuming a fragmentation-limited background dust population, we developed a Monte Carlo approach that follows single monomers as they move through a vertical column of gas and become incorporated in different aggregates as they undergo sticking and fragmenting collisions. We find that (sub)micron grains are not necessarily well-mixed vertically, but can become trapped in a thin layer with a scale-height that is significantly smaller than that of the gas. This collisional trapping occurs when the timescale for diffusion is comparable to or longer than the collision timescale in the midplane and its effect is strongest when the most massive particles in the size-distribution show significant settling. Based on simulations and analytical considerations we conclude that for typical dust-to-gas ratios and turbulence levels, the collisional trapping of small grains should be a relatively common phenomenon. The absence of trapping could then indicate a low dust-to-gas ratio, possibly because a large portion of the dust mass has been removed through radial drift or is locked up in planetesimals.
△ Less
Submitted 8 March, 2016;
originally announced March 2016.
-
A panoptic model for planetesimal formation and pebble delivery
Authors:
Sebastiaan Krijt,
Chris W. Ormel,
Carsten Dominik,
Alexander G. G. M. Tielens
Abstract:
The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the…
▽ More
The journey from dust particle to planetesimal involves physical processes acting on scales ranging from micrometers (the sticking and restructuring of aggregates) to hundreds of astronomical units (the size of the turbulent protoplanetary nebula). Considering these processes simultaneously is essential when studying planetesimal formation. We develop a novel, global, semi-analytical model for the evolution of the mass-dominating dust particles in a turbulent protoplanetary disk that takes into account the evolution of the dust surface density while preserving the essential characteristics of the porous coagulation process. This panoptic model is used to study the growth from submicron to planetesimal sizes in disks around Sun-like stars. For highly porous ices, unaffected by collisional fragmentation and erosion, rapid growth to planetesimal sizes is possible in a zone stretching out to ${\sim}10\mathrm{~AU}$ for massive disks. When porous coagulation is limited by erosive collisions, the formation of planetesimals through direct coagulation is not possible, but the creation of a large population of aggregates with Stokes numbers close to unity might trigger the streaming instability (SI). However, we find that reaching conditions necessary for SI is difficult and limited to dust-rich disks, (very) cold disks, or disks with weak turbulence. Behind the snow-line, porosity-driven aggregation of icy grains results in rapid (${\sim}10^{4}\mathrm{~yr}$) formation of planetesimals. If erosive collisions prevent this, SI might be triggered for specific disk conditions. The numerical approach introduced in this work is ideally suited for studying planetesimal formation and pebble delivery simultaneously and will help build a coherent picture of the start of the planet formation process.
△ Less
Submitted 24 November, 2015;
originally announced November 2015.
-
Erosion and the limits to planetesimal growth
Authors:
Sebastiaan Krijt,
Chris W. Ormel,
Carsten Dominik,
Alexander G. G. M. Tielens
Abstract:
The coagulation of microscopic dust into planetesimals is the first step towards planet formation. The size and shape of the growing aggregates determine the efficiency of this early growth. It has been proposed that fluffy ice aggregates can grow very efficiently, suffering less from the bouncing and radial drift barriers. While the collision velocity between icy aggregates of similar size is tho…
▽ More
The coagulation of microscopic dust into planetesimals is the first step towards planet formation. The size and shape of the growing aggregates determine the efficiency of this early growth. It has been proposed that fluffy ice aggregates can grow very efficiently, suffering less from the bouncing and radial drift barriers. While the collision velocity between icy aggregates of similar size is thought to stay below the fragmentation threshold, they may nonetheless lose mass from collisions with much smaller projectiles. We investigate the effect of these erosive collisions on the ability of porous ice aggregates to cross the radial drift barrier. We develop a Monte Carlo code that calculates the evolution of the growing aggregates, while resolving the entire mass distribution at all times. The aggregate's porosity is treated independently of its mass, and is determined by collisions, gas compaction, and eventually self-gravity compaction. For erosion threshold velocities of 20-40 m/s, high-velocity collisions with small projectiles prevent the largest aggregates from growing when they start to drift. In these cases, our local simulations result in a steady-state distribution, with the majority of the dust mass in particles with Stokes numbers close to unity. Only for the highest erosion threshold considered (60 m/s), do porous aggregates manage to cross the radial drift barrier in the inner 10 AU of MMSN-like disks. Erosive collisions are more effective in limiting the growth than fragmentary collisions between similar-size particles. Conceivably, erosion limits the growth before the radial drift barrier, although the robustness of this statement depends on (uncertain) material properties of icy aggregates. If erosion inhibits planetesimal formation through direct sticking, the sea of ${\sim}10^9$ g, highly porous particles appears well-suited for triggering streaming instability.
△ Less
Submitted 11 December, 2014;
originally announced December 2014.
-
A dearth of small particles in debris disks: An energy-constrained smallest fragment size
Authors:
Sebastiaan Krijt,
Mihkel Kama
Abstract:
A prescription for the fragment size distribution resulting from dust grain collisions is essential when modelling a range of astrophysical systems, such as debris disks and planetary rings. While the slope of the fragment size distribution and the size of the largest fragment are well known, the behaviour of the distribution at the small size end is theoretically and experimentally poorly underst…
▽ More
A prescription for the fragment size distribution resulting from dust grain collisions is essential when modelling a range of astrophysical systems, such as debris disks and planetary rings. While the slope of the fragment size distribution and the size of the largest fragment are well known, the behaviour of the distribution at the small size end is theoretically and experimentally poorly understood. This leads debris disk codes to generally assume a limit equal to, or below, the radiation blow-out size. We use energy conservation to analytically derive a lower boundary of the fragment size distribution for a range of collider mass ratios. Focussing on collisions between equal-sized bodies, we apply the method to debris disks. For a given collider mass, the size of the smallest fragments is found to depend on collision velocity, material parameters, and the size of the largest fragment. We provide a physically motivated recipe for the calculation of the smallest fragment, which can be easily implemented in codes for modelling collisional systems. For plausible parameters, our results are consistent with the observed predominance of grains much larger than the blow-out size in Fomalhaut's main belt and in the Herschel cold debris disks.
△ Less
Submitted 27 May, 2014;
originally announced May 2014.
-
Blowing in the wind: The dust wave surrounding $σ$ Ori AB
Authors:
B. B. Ochsendorf,
N. L. J. Cox,
S. Krijt,
F. Salgado,
O. Berné,
J. P. Bernard,
L. Kaper,
A. G. G. M. Tielens
Abstract:
Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $σ$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitativ…
▽ More
Observations with the Spitzer Space Telescope and the WISE satellite have revealed a prominent arc-like structure at 50" ($\simeq0.1$ pc) from the O9.5V/B0.5V system $σ$ Ori AB. We attribute this dust structure to the interaction of radiation pressure from the star with dust carried along by the IC 434 photo-evaporative flow of ionized gas from the dark cloud L1630. We have developed a quantitative model for the interaction of a dusty ionized flow with nearby (massive) stars where radiation pressure stalls dust, piling it up at an appreciable distance (> 0.1 pc), and force it to flow around the star. The model demonstrates that for the conditions in IC 434, the gas will decouple from the dust and will keep its original flow lines. We argue that this dust structure is the first example of a dust wave created by a massive star moving through the interstellar medium. Dust waves (and bow waves) stratify dust grains according to their radiation pressure opacity, which reflects the size distribution and composition of the grain material. Comparison of our model with observations implies that dust-gas coupling through Coulomb interaction is less important than previously thought, challenging our understanding of grain dynamics in hot, ionized regions of space. We describe the difference between dust (and bow) waves and classical bow shocks. We conclude that dust waves and bow waves should be common around stars showing the weak-wind phenomenon and that these structures are best observed at mid-IR to FIR wavelengths. In particular, dust waves and bow waves are most efficiently formed around weak-wind stars moving through a high density medium. Moreover, they provide a unique opportunity to study the direct interaction between a (massive) star and its immediate surroundings.
△ Less
Submitted 8 July, 2014; v1 submitted 28 January, 2014;
originally announced January 2014.
-
Erosion of dust aggregates
Authors:
Alexander Seizinger,
Sebastiaan Krijt,
Wilhelm Kley
Abstract:
Aims: The aim of this work is to gain a deeper insight into how much different aggregate types are affected by erosion. Especially, it is important to study the influence of the velocity of the impacting projectiles. We also want to provide models for dust growth in protoplanetary disks with simple recipes to account for erosion effects.
Methods: To study the erosion of dust aggregates we employ…
▽ More
Aims: The aim of this work is to gain a deeper insight into how much different aggregate types are affected by erosion. Especially, it is important to study the influence of the velocity of the impacting projectiles. We also want to provide models for dust growth in protoplanetary disks with simple recipes to account for erosion effects.
Methods: To study the erosion of dust aggregates we employed a molecular dynamics approach that features a detailed micro-physical model of the interaction of spherical grains. For the first time, the model has been extended by introducing a new visco-elastic damping force which requires a proper calibration. Afterwards, different sample generation methods were used to cover a wide range of aggregate types.
Results: The visco-elastic damping force introduced in this work turns out to be crucial to reproduce results obtained from laboratory experiments. After proper calibration, we find that erosion occurs for impact velocities of 5 m/s and above. Though fractal aggregates as formed during the first growth phase are most susceptible to erosion, we observe erosion of aggregates with rather compact surfaces as well.
Conclusions: We find that bombarding a larger target aggregate with small projectiles results in erosion for impact velocities as low as a few m/s. More compact aggregates suffer less from erosion. With increasing projectile size the transition from accretion to erosion is shifted to higher velocities. This allows larger bodies to grow through high velocity collisions with smaller aggregates.
△ Less
Submitted 16 October, 2013;
originally announced October 2013.
-
Stirring up the dust: A dynamical model for halo-like dust clouds in transitional disks
Authors:
Sebastiaan Krijt,
Carsten Dominik
Abstract:
A small number of young stellar objects show signs of a halo-like structure of optically thin dust. This halo or torus is located within a few AU of the star, but its origin has not yet been understood. A dynamically excited cloud of planetesimals colliding to eventually form dust could produce such a structure. The cause of the dynamical excitation could be one or more planets. This work investig…
▽ More
A small number of young stellar objects show signs of a halo-like structure of optically thin dust. This halo or torus is located within a few AU of the star, but its origin has not yet been understood. A dynamically excited cloud of planetesimals colliding to eventually form dust could produce such a structure. The cause of the dynamical excitation could be one or more planets. This work investigates an inwardly migrating planet that is dynamically scattering planetesimals as a possible cause for the observed structures. If this mechanism is responsible, the observed halo-like structure could be used to infer the existence of planets in these systems. We present analytical estimates on the maximum inclination reached owing to dynamical interactions between planetesimals and a migrating planet. A symplectic integrator is used to simulate the effect of a migrating planet on a population of planetesimals. It is found that an inwardly migrating planet is only able to scatter the material it encounters to highly-inclined orbits if the material is on an eccentric orbit. Such eccentric orbits can be the result of resonance trapping and eccentricity pumping. Simulations show that for a certain range of migration rates and planet masses, resonance capture combined with planetary migration indeed causes the planetesimals to reach eccentric orbits and subsequently get scattered to highly-inclined orbits. The size distribution of the resulting dust is calculated determined to find the total mass and optical depth, which are found to compare reasonably well with the observed structures. Dynamical scattering of planetesimals caused by a planet migrating in, followed by the grinding down of these planetesimals to dust grains, appears to be a promising explanation for the inferred circumstellar dust clouds. Further study is needed to see if the haloes can be used to infer the presence of planets.
△ Less
Submitted 26 May, 2011;
originally announced May 2011.