-
A Dust-Trapping Ring in the Planet-Hosting Disk of Elias 2-24
Authors:
Adolfo S. Carvalho,
Laura M. Perez,
Anibal Sierra,
Maria Jesus Mellado,
Lynne A. Hillenbrand,
Sean Andrews,
Myriam Benisty,
Tilman Birnstiel,
John M. Carpenter,
Viviana V. Guzman,
Jane Huang,
Andrea Isella,
Nicolas Kurtovic,
Luca Ricci,
David J. Wilner
Abstract:
Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid gr…
▽ More
Rings and gaps are among the most widely observed forms of substructure in protoplanetary disks. A gap-ring pair may be formed when a planet carves a gap in the disk, which produces a local pressure maximum following the gap that traps inwardly drifting dust grains and appears as a bright ring due to the enhanced dust density. A dust-trapping ring would provide a promising environment for solid growth and possibly planetesimal production via the streaming instability. We present evidence of dust trapping in the bright ring of the planet-hosting disk Elias 2-24, from the analysis of 1.3 mm and 3 mm ALMA observations at high spatial resolution (0.029 arcsec, 4.0 au). We leverage the high spatial resolution to demonstrate that larger grains are more efficiently trapped and place constraints on the local turbulence ($8 \times 10^{-4} < α_\mathrm{turb} < 0.03$) and the gas-to-dust ratio ($Σ_g / Σ_d < 30$) in the ring. Using a scattering-included marginal probability analysis we measure a total dust disk mass of $M_\mathrm{dust} = 13.8^{+0.7}_{-0.5} \times 10^{-4} \ M_\odot$. We also show that at the orbital radius of the proposed perturber, the gap is cleared of material down to a flux contrast of 10$^{-3}$ of the peak flux in the disk.
△ Less
Submitted 18 June, 2024;
originally announced June 2024.
-
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.
-
Vertical shear instability with dust evolution and consistent cooling times. On the importance of the initial dust distribution
Authors:
Thomas Pfeil,
Til Birnstiel,
Hubert Klahr
Abstract:
Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust.
Aims. We aim…
▽ More
Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust.
Aims. We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales.
Results. We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 $μ$m) generally leads to a less turbulent outcome. While the inner disk regions (within $\sim$ 70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk.
Conclusions. Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially $\gtrsim$ 100 $μ$m and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): PDS 111, an old T Tauri star with a young-looking disk
Authors:
Annelotte Derkink,
Christian Ginski,
Paola Pinilla,
Nicolas Kurtovic,
Lex Kaper,
Alex de Koter,
Per-Gunnar Valegård,
Eric Mamajek,
Frank Backs,
Myriam Benisty,
Til Birnstiel,
Gabriele Columba,
Carsten Dominik,
Antonio Garufi,
Michiel Hogerheijde,
Rob van Holstein,
Jane Huang,
François Ménard,
Christian Rab,
María Claudia Ramírez-Tannus,
Álvaro Ribas,
Jonathan P. Williams,
Alice Zurlo
Abstract:
The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary…
▽ More
The interplay between T Tauri stars and their circumstellar disks, and how this impacts the onset of planet formation has yet to be established. We studied a seemingly old T Tauri star, PDS 111, and its disk. We analyzed optical, infrared, and sub-millimeter observations obtained with VLT/X-shooter, Mercator/HERMES, TESS, VLT/SPHERE, and ALMA, providing a new view on PDS 111 and its protoplanetary disk. The multi-epoch spectroscopy yields photospheric lines to classify the star, and emission lines to study variability in the hot inner disk and to determine the mass-accretion rate. The SPHERE and ALMA observations are used to characterize the dust distribution of the small and large grains, respectively. PDS 111 is a weak-line T Tauri star with spectral type G2, exhibits strong H$α$ variability and with a low mass-accretion rate of $1-5\times10^{-10}$\,M$_{\odot}$\,yr$^{-1}$. We measured an age of the system of 15.9$^{+1.7}_{-3.7}$ Myr using pre-main sequence tracks. The SPHERE observations show a strongly flaring disk with an asymmetric substructure. The ALMA observations reveal a 30 au cavity in the dust continuum emission with a low contrast asymmetry in the South-West of the disk and a dust disk mass of 45.8\,$M_\oplus$. The $^{12}$CO radial extension is at least three times larger than that of the dust emission. Although the measured age is younger than suggested in literature, PDS 111 still seems relatively old; this provides insight into disk properties at an advanced stage of pre-main sequence evolution. The characteristics of this disk are very similar to its younger counterparts: strongly flaring, an average disk mass, a typical radial extent of the disk gas and dust, and the presence of common substructures. This suggests that disk evolution has not significantly changed the disk properties. These results show similarities with the "Peter Pan disks" around M-dwarfs.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
Population Synthesis Models Indicate a Need for Early and Ubiquitous Disk Substructures
Authors:
Luca Delussu,
Tilman Birnstiel,
Anna Miotello,
Paola Pinilla,
Giovanni Rosotti,
Sean M. Andrews
Abstract:
Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust…
▽ More
Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslantα\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
The SPHERE view of the Orion star-forming region
Authors:
P. -G. Valegard,
C. Ginski,
A. Derkink,
A. Garufi,
C. Dominik,
A. Ribas,
J. P. Williams,
M. Benisty,
T. Birnstiel,
S. Facchini,
G. Columba,
M. Hogerheijde,
R. G. Van Holstein,
J. Huang,
M. Kenworthy,
C. F. Manara,
P. Pinilla,
Ch. Rab,
R. Sulaiman,
A. Zurlo
Abstract:
We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analys…
▽ More
We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion Star forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We use polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterize the disk orientation, radius and contrast. We analyse the disks in context of the stellar parameters and the environment of the Orion star-forming region. We use ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore use a combination of new and archival ALMA mm-continuum observations to characterize the dust masses present in the circumstellar disks. Within our sample we detect extended circumstellar disks in 10 of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori and V1012 Ori) and show scattered light asymmetries which may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high resolution imaging observations are also sensitive to close (sub)stellar companions and we detect 9 such objects in our sample of which 5 were previously unknown. We find in particular a possible sub-stellar companion (either a very low mass star or a high mass brown dwarf) 137 au from the star RY Ori. We find a strong anti-correlation between disk detection and multiplicity, with only 2 of our 10 disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimetre flux suggesting that disks that have a high dust content are typically bright in near-infrared scattered light. Conversely we do not find significant correlations between scattered light contrast of the disks and the stellar mass or age.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
The SPHERE view of the Chamaeleon I star-forming region
Authors:
C. Ginski,
A. Garufi,
M. Benisty,
R. Tazaki,
C. Dominik,
A. Ribas,
N. Engler,
T. Birnstiel,
G. Chauvin,
G. Columba,
S. Facchini,
A. Goncharov,
J. Hagelberg,
T. Henning,
M. Hogerheijde,
R. G. van Holstein,
J. Huang,
T. Muto,
P. Pinilla,
K. Kanagawa,
S. Kim,
N. Kurtovic,
M. Langlois,
C. Manara,
J. Milli
, et al. (10 additional authors not shown)
Abstract:
We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems…
▽ More
We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems. In 13 of the 20 systems included in this study we detected resolved scattered light signals from circumstellar dust. For the CR Cha, CT Cha, CV Cha, SY Cha, SZ Cha, and VZ Cha systems we present the first detailed descriptions of the disks in scattered light. The observations found typically smooth or faint disks, often with little substructure, with the notable exceptions of SZ Cha, which shows an extended multiple-ringed disk, and WW Cha, which shows interaction with the cloud environment. New high S/N K- band observations of the HD 97048 system in our survey reveal a significant brightness asymmetry that may point to disk misalignment and subsequent shadowing of outer disk regions, possibly related to the suggested planet candidate in the disk. We resolve for the first time the stellar binary in the CS Cha system. Multiple wavelength observations of the disk around CS Cha have revealed that the system contains small, compact dust grains that may be strongly settled, consistent with numerical studies of circumbinary disks. We find in our sample that there is a strong anti-correlation between the presence of a (close) stellar companion and the detection of circumstellar material with five of our seven nondetections located in binary systems.
△ Less
Submitted 4 March, 2024;
originally announced March 2024.
-
TIPSY: Trajectory of Infalling Particles in Streamers around Young stars. Dynamical analysis of the streamers around S CrA and HL Tau
Authors:
Aashish Gupta,
Anna Miotello,
Jonathan P. Williams,
Til Birnstiel,
Michael Kuffmeier,
Hsi-Wei Yen
Abstract:
Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs.
Aims. In order to ascertain the infalling nature of observed streamer-like structures and…
▽ More
Context. Elongated trails of infalling gas, often referred to as "streamers," have recently been observed around young stellar objects (YSOs) at different evolutionary stages. This asymmetric infall of material can significantly alter star and planet formation processes, especially in the more evolved YSOs.
Aims. In order to ascertain the infalling nature of observed streamer-like structures and then systematically characterize their dynamics, we developed the code TIPSY (Trajectory of Infalling Particles in Streamers around Young stars).
Methods. Using TIPSY, the streamer molecular line emission is first isolated from the disk emission. Then the streamer emission, which is effectively a point cloud in three-dimensional (3D) position-position-velocity space, is simplified to a curve-like representation. The observed streamer curve is then compared to the theoretical trajectories of infalling material. The best-fit trajectories are used to constrain streamer features, such as the specific energy, the specific angular momenta, the infall timescale, and the 3D morphology.
Results. We used TIPSY to fit molecular-line ALMA observations of streamers around a Class II binary system, S CrA, and a Class I/II protostar, HL Tau. Our results indicate that both of the streamers are consistent with infalling motion. TIPSY results and mass estimates suggest that S CrA and HL Tau are accreting material at a rate of $\gtrsim27$ M$_{jupiter}$ Myr$^{-1}$ and $\gtrsim5$ M$_{jupiter}$ Myr$^{-1}$, respectively, which can significantly increase the mass budget available to form planets.
Conclusions. TIPSY can be used to assess whether the morphology and kinematics of observed streamers are consistent with infalling motion and to characterize their dynamics, which is crucial for quantifying their impact on the protostellar systems.
△ Less
Submitted 18 January, 2024;
originally announced January 2024.
-
Dust growth and evolution in protoplanetary disks
Authors:
Tilman Birnstiel
Abstract:
Over the past decade, advancement of observational capabilities, specifically the Atacama Large Millimeter/submillimeter Array (ALMA) and SPHERE instrument, alongside theoretical innovations like pebble accretion, have reshaped our understanding of planet formation and the physics of protoplanetary disks. Despite this progress, mysteries persist along the winded path of micrometer-sized dust, from…
▽ More
Over the past decade, advancement of observational capabilities, specifically the Atacama Large Millimeter/submillimeter Array (ALMA) and SPHERE instrument, alongside theoretical innovations like pebble accretion, have reshaped our understanding of planet formation and the physics of protoplanetary disks. Despite this progress, mysteries persist along the winded path of micrometer-sized dust, from the interstellar medium, through transport and growth in the protoplanetary disk, to becoming gravitationally bound bodies. This review outlines our current knowledge of dust evolution in circumstellar disks, yielding the following insights:
$\bullet$ Theoretical and laboratory studies have accurately predicted the growth of dust particles to sizes that are susceptible to accumulation through transport processes like radial drift and settling.
$\bullet$ Critical uncertainties in that process remain the level of turbulence, the threshold collision velocities at which dust growth stalls, and the evolution of dust porosity.
$\bullet$ Symmetric and asymmetric substructure are widespread. Dust traps appear to be solving several long-standing issues in planet formation models, and they are observationally consistent with being sites of active planetesimal formation.
$\bullet$ In some instances, planets have been identified as the causes behind substructures. This underlines the need to study earlier stages of disks to understand how planets can form so rapidly.
In the future, better probes of the physical conditions in optically thick regions, including densities, turbulence strength, kinematics, and particle properties will be essential for unraveling the physical processes at play.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
A Collection of German Science Interests in the Next Generation Very Large Array
Authors:
M. Kadler,
D. A. Riechers,
J. Agarwal,
A. -K. Baczko,
H. Beuther,
F. Bigiel,
T. Birnstiel,
B. Boccardi,
D. J. Bomans,
L. Boogaard,
T. T. Braun,
S. Britzen,
M. Brüggen,
A. Brunthaler,
P. Caselli,
D. Elsässer,
S. von Fellenberg,
M. Flock,
C. M. Fromm,
L. Fuhrmann,
P. Hartogh,
M. Hoeft,
R. P. Keenan,
Y. Kovalev,
K. Kreckel
, et al. (66 additional authors not shown)
Abstract:
The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com…
▽ More
The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions.
△ Less
Submitted 18 June, 2024; v1 submitted 16 November, 2023;
originally announced November 2023.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): HD 34700 A unveils an inner ring
Authors:
G. Columba,
E. Rigliaco,
R. Gratton,
D. Mesa,
V. D'Orazi,
C. Ginski,
N. Engler,
J. P. Williams,
J. Bae,
M. Benisty,
T. Birnstiel,
P. Delorme,
C. Dominik,
S. Facchini,
F. Menard,
P. Pinilla,
C. Rab,
Á. Ribas,
V. Squicciarini,
R. G. van Holstein,
A. Zurlo
Abstract:
Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms.
Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures.
Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and…
▽ More
Context. The study of protoplanetary disks is fundamental to understand their evolution and interaction with the surrounding environment, and to constrain planet formation mechanisms.
Aims. We aim at characterising the young binary system HD 34700 A, which shows a wealth of structures.
Methods. Taking advantage of the high-contrast imaging instruments SPHERE at the VLT, LMIRCam at the LBT, and of ALMA observations, we analyse this system at multiple wavelengths. We study the rings and spiral arms morphology and the scattering properties of the dust. We discuss the possible causes of all the observed features.
Results. We detect for the first time, in the H$α$ band, a ring extending from $\sim$65 au to ${\sim}$120 au, inside the ring already known from recent studies. These two have different physical and geometrical properties. Based on the scattering properties, the outer ring may consist of grains of typical size $a_{out} > 4 μm$, while the inner ring of smaller grains ($a_{in} <= 0.4 {μm}$). Two extended logarithmic spiral arms stem from opposite sides of the disk. The outer ring appears as a spiral arm itself, with a variable radial distance from the centre and extended substructures. ALMA data confirm the presence of a millimetric dust substructure centred just outside the outer ring, and detect misaligned gas rotation patterns for HD 34700 A and B.
Conclusions. The complexity of HD 34700 A, revealed by the variety of observed features, suggests the existence of one or more disk-shaping physical mechanisms. Possible scenarios, compatible with our findings, involve the presence inside the disk of a yet undetected planet of several Jupiter masses and the system interaction with the surroundings by means of gas cloudlet capture or flybys. Further observations with JWST/MIRI or ALMA (gas kinematics) could shed more light on these.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Dust Coagulation Reconciles Protoplanetary Disk Observations with the Vertical Shear Instability. I. Dust Coagulation and the VSI Dead Zone
Authors:
Thomas Pfeil,
Til Birnstiel,
Hubert Klahr
Abstract:
Protoplanetary disks exhibit a vertical gradient in angular momentum, rendering them susceptible to the Vertical Shear Instability (VSI). The most important condition for the onset of this mechanism is a short timescale of thermal relaxation ($\lesssim 0.1$ orbital timescales). Simulations of fully VSI active disks are characterized by turbulent, vertically extended dust layers. This is in contrad…
▽ More
Protoplanetary disks exhibit a vertical gradient in angular momentum, rendering them susceptible to the Vertical Shear Instability (VSI). The most important condition for the onset of this mechanism is a short timescale of thermal relaxation ($\lesssim 0.1$ orbital timescales). Simulations of fully VSI active disks are characterized by turbulent, vertically extended dust layers. This is in contradiction with recent observations of the outer regions of some protoplanetary disks, which appear highly settled. In this work, we demonstrate that the process of dust coagulation can diminish the cooling rate of the gas in the outer disk and extinct the VSI activity. Our findings indicate that the turbulence strength is especially susceptible to variations in the fragmentation velocity of the grains. A small fragmentation velocity of $\approx 100 \mathrm{\, cm \,s^{-1}}$ results in a fully turbulent simulation, whereas a value of $\approx 400 \mathrm{\, cm \,s^{-1}}$ results in a laminar outer disk, being consistent with observations. We show that VSI turbulence remains relatively unaffected by variations in the maximum particle size in the inner disk regions. However, we find that dust coagulation can significantly suppress the occurrence of VSI turbulence at larger distances from the central star.
△ Less
Submitted 15 December, 2023; v1 submitted 11 October, 2023;
originally announced October 2023.
-
Millimeter emission in photoevaporating disks is determined by early substructures
Authors:
Matías Gárate,
Til Birnstiel,
Paola Pinilla,
Sean M. Andrews,
Raphael Franz,
Sebastian Markus Stammler,
Giovanni Picogna,
Barbara Ercolano,
Anna Miotello,
Nicolás T. Kurtovic
Abstract:
[abridged]Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. We studied how the presence of early substructures affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble…
▽ More
[abridged]Photoevaporation and dust-trapping are individually considered to be important mechanisms in the evolution and morphology of protoplanetary disks. We studied how the presence of early substructures affects the evolution of the dust distribution and flux in the millimeter continuum of disks that are undergoing photoevaporative dispersal. We also tested if the predicted properties resemble those observed in the population of transition disks. We used the numerical code Dustpy to simulate disk evolution considering gas accretion, dust growth, dust-trapping at substructures, and mass loss due to X-ray and EUV (XEUV) photoevaporation and dust entrainment. Then, we compared how the dust mass and millimeter flux evolve for different disk models. We find that, during photoevaporative dispersal, disks with primordial substructures retain more dust and are brighter in the millimeter continuum than disks without early substructures, regardless of the photoevaporative cavity size. Once the photoevaporative cavity opens, the estimated fluxes for the disk models that are initially structured are comparable to those found in the bright transition disk population ($F_\textrm{mm} > 30\, \textrm{mJy}$), while the disk models that are initially smooth have fluxes comparable to the transition disks from the faint population ($F_\textrm{mm} < 30\, \textrm{mJy}$), suggesting a link between each model and population. Our models indicate that the efficiency of the dust trapping determines the millimeter flux of the disk, while the gas loss due to photoevaporation controls the formation and expansion of a cavity, decoupling the mechanisms responsible for each feature. In consequence, even a planet with a mass comparable to Saturn could trap enough dust to reproduce the millimeter emission of a bright transition disk, while its cavity size is independently driven by photoevaporative dispersal.
△ Less
Submitted 15 September, 2023;
originally announced September 2023.
-
Three-Dimensional Dust Stirring by a Giant Planet Embedded in a Protoplanetary Disk
Authors:
Fabian Binkert,
Judit Szulágyi,
Til Birnstiel
Abstract:
The motion of solid particles embedded in gaseous protoplanetary disks is influenced by turbulent fluctuations. Consequently, the dynamics of moderately to weakly coupled solids can be distinctly different from the dynamics of the gas. Additionally, gravitational perturbations from an embedded planet can further impact the dynamics of solids. In this work, we investigate the combined effects of tu…
▽ More
The motion of solid particles embedded in gaseous protoplanetary disks is influenced by turbulent fluctuations. Consequently, the dynamics of moderately to weakly coupled solids can be distinctly different from the dynamics of the gas. Additionally, gravitational perturbations from an embedded planet can further impact the dynamics of solids. In this work, we investigate the combined effects of turbulent fluctuations and planetary dust stirring in a protoplanetary disk on three-dimensional dust morphology and on synthetic ALMA continuum observations. We carry out 3D radiative two-fluid (gas+1-mm-dust) hydrodynamic simulations in which we explicitly model the gravitational perturbation of a Jupiter-mass planet. We derived a new momentum-conserving turbulent diffusion model that introduces a turbulent pressure to the pressureless dust fluid to capture the turbulent transport of dust. The model implicitly captures the effects of orbital oscillations and reproduces the theoretically predicted vertical settling-diffusion equilibrium. We find a Jupiter-mass planet to produce distinct and large-scale three-dimensional flow structures in the mm-size dust, which vary strongly in space. We quantify these effects by locally measuring an effective vertical diffusivity (equivalent alpha) and find azimuthally averaged values in a range $δ_\mathrm{eff}\sim5\cdot 10^{-3} - 2\cdot 10^{-2}$ and local peaks at values of up to $δ_\mathrm{eff}\sim3\cdot 10^{-1}$. In synthetic ALMA continuum observations of inclined disks, we find effects of turbulent diffusion to be observable, especially at disk edges, and effects of planetary dust stirring in edge-on observations.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
The impact of dust evolution on the dead zone outer edge in magnetized protoplanetary disks
Authors:
Timmy N. Delage,
Matías Gárate,
Satoshi Okuzumi,
Chao-Chin Yang,
Paola Pinilla,
Mario Flock,
Sebastian Markus Stammler,
Tilman Birnstiel
Abstract:
[Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the…
▽ More
[Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles. Since the dust content of the disk decreases over million years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced $α$-parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there is enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once it is no longer the case, it is expected to be controlled by gas evolution and occurs on a viscous evolution timescale.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
Carbon Depletion in the Early Solar System
Authors:
Fabian Binkert,
Til Birnstiel
Abstract:
Earth and other rocky objects in the inner Solar System are depleted in carbon compared to objects in the outer Solar System, the Sun, or the ISM. It is believed that this is a result of the selective removal of refractory carbon from primordial circumstellar material. In this work, we study the irreversible release of carbon into the gaseous environment via photolysis and pyrolysis of refractory…
▽ More
Earth and other rocky objects in the inner Solar System are depleted in carbon compared to objects in the outer Solar System, the Sun, or the ISM. It is believed that this is a result of the selective removal of refractory carbon from primordial circumstellar material. In this work, we study the irreversible release of carbon into the gaseous environment via photolysis and pyrolysis of refractory carbonaceous material during the disk phase of the early Solar System. We analytically solve the one-dimensional advection equation and derive an explicit expression that describes the depletion of carbonaceous material in solids under the influence of radial and vertical transport. We find both depletion mechanisms individually fail to reproduce Solar System abundances under typical conditions. While radial transport only marginally restricts photodecomposition, it is the inefficient vertical transport that limits carbon depletion under these conditions. We show explicitly that an increase in the vertical mixing efficiency, and/or an increase in the directly irradiated disk volume, favors carbon depletion. Thermal decomposition requires a hot inner disk (> 500 K) beyond 3 AU to deplete the formation region of Earth and chondrites. We find FU Ori-type outbursts to produce these conditions such that moderately refractory compounds are depleted. However, such outbursts likely do not deplete the most refractory carbonaceous compounds beyond the innermost disk region. Hence, the refractory carbon abundance at 1 AU typically does not reach terrestrial levels. Nevertheless, under specific conditions, we find photolysis and pyrolysis combined to reproduce Solar System abundances.
△ Less
Submitted 13 January, 2023;
originally announced January 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.
-
Reflections on nebulae around young stars: A systematic search for late-stage infall of material onto Class II disks
Authors:
Aashish Gupta,
Anna Miotello,
Carlo F. Manara,
Jonathan P. Williams,
Stefano Facchini,
Giacomo Beccari,
Til Birnstiel,
Christian Ginski,
Alvaro Hacar,
Michael Küffmeier,
Leonardo Testi,
Lukasz Tychoniec,
Hsi-Wei Yen
Abstract:
Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation.
Aims. In order to systematically study late infall of gas onto disks, we i…
▽ More
Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation.
Aims. In order to systematically study late infall of gas onto disks, we identify candidate Class II sources in close vicinity to a reflection nebula (RN) that may be undergoing this process.
Methods. First we targeted Class II sources with known kilo-au scale gas structures - possibly due to late infall of material - and we searched for RNe in their vicinity in optical and near-infrared images. Second, we compiled a catalogue of Class II sources associated with RNe and looked for the large-scale CO structures in archival ALMA data. Using the catalogues of protostars and RNe, we also estimated the probability of Class II sources interacting with surrounding material.
Results. All of the sources with large-scale gas structures also exhibit some reflection nebulosity in their vicinity. Similarly, at least five Class II objects associated with a prominent RNe, and for which adequate ALMA observations are available, were found to have spirals or stream-like structures which may be due to late infall. We report the first detection of these structures around S CrA.
Conclusions. Our results suggest that a non-negligible fraction of Class II disks in nearby star-forming regions may be associated with RNe and could therefore be undergoing late accretion of gas. Surveys of RNe and kilo-au scale gas structures around Class II sources will allow us to better understand the frequency and impact of late-infall phenomena.
△ Less
Submitted 8 January, 2023;
originally announced January 2023.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Characterization of the young star T CrA and its circumstellar environment
Authors:
E. Rigliaco,
R. Gratton,
S. Ceppi,
C. Ginski,
M. Hogerheijde,
M. Benisty,
T. Birnstiel,
M. Dima,
S. Facchini,
A. Garufi,
J. Bae,
M. Langlois,
G. Lodato,
E. Mamajek,
C. F. Manara,
F. Ménard,
Á. Ribas,
A. Zurlo
Abstract:
Birth environments of young stars have strong imprints on the star itself and their surroundings. We present a detailed analysis of the wealthy circumstellar environment around the young Herbig Ae/Be star TCrA. Our aim is to understand the nature of the stellar system and the extended circumstellar structures as seen in scattered light images. We conduct our analysis combining archival data, and n…
▽ More
Birth environments of young stars have strong imprints on the star itself and their surroundings. We present a detailed analysis of the wealthy circumstellar environment around the young Herbig Ae/Be star TCrA. Our aim is to understand the nature of the stellar system and the extended circumstellar structures as seen in scattered light images. We conduct our analysis combining archival data, and new adaptive optics high-contrast and high-resolution images. The scattered light images reveal the presence of a complex environment composed of a bright forward scattering rim of the disk's surface that is seen at very high inclination, a dark lane of the disk midplane, bipolar outflows, and streamer features likely tracing infalling material from the surrounding birth cloud onto the disk. The analysis of the light curve suggests the star is a binary with a period of 29.6yrs. The comparison of the scattered light images with ALMA continuum and 12CO line emission shows the disk is in keplerian rotation, with the northern side of the outflowing material receding, while the southern side approaching the observer. The disk is itself seen edge-on. The direction of the outflows seen in scattered light is in agreement with the direction of the more distant molecular hydrogen emission-line objects (MHOs) associated to the star. Modeling of the SED using a radiative transfer scheme well agrees with the proposed configuration, as well as the hydrodynamical simulation performed using a Smoothed Particle Hydrodynamics code. We find evidence of streamers of accreting material around TCrA. These streamers connect the filament along which TCrA is forming with the outer parts of the disk, suggesting that the strong misalignment between the inner and outer disk is due to a change in the direction of the angular momentum of the material accreting on the disk during the late phase of star formation.
△ Less
Submitted 4 January, 2023;
originally announced January 2023.
-
Comets and Planetesimal Formation
Authors:
Jacob B. Simon,
Jürgen Blum,
Til Birnstiel,
David Nesvorný
Abstract:
In this chapter, we review the processes involved in the formation of planetesimals and comets. We will start with a description of the physics of dust grain growth and how this is mediated by gas-dust interactions in planet-forming disks. We will then delve into the various models of planetesimal formation, describing how these planetesimals form as well as their resulting structure. In doing so,…
▽ More
In this chapter, we review the processes involved in the formation of planetesimals and comets. We will start with a description of the physics of dust grain growth and how this is mediated by gas-dust interactions in planet-forming disks. We will then delve into the various models of planetesimal formation, describing how these planetesimals form as well as their resulting structure. In doing so, we focus on and compare two paradigms for planetesimal formation: the gravitational collapse of particle over-densities (which can be produced by a variety of mechanisms) and the growth of particles into planetesimals via collisional and gravitational coagulation. Finally, we compare the predictions from these models with data collected by the Rosetta and New Horizons missions and that obtained via observations of distant Kuiper Belt Objects.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Constraining the turbulence and the dust disk in IM Lup: onset of planetesimal formation
Authors:
Riccardo Franceschi,
Tilman Birnstiel,
Thomas Henning,
Anirudh Sharma
Abstract:
Observations of protoplanetary disks provide information on planet formation and the reasons for the diversity of planetary systems. The key to understanding planet formation is the study of dust evolution from small grains to pebbles. Smaller grains are well-coupled to the gas dynamics, and their distribution is significantly extended above the disk midplane. Larger grains settle much faster and…
▽ More
Observations of protoplanetary disks provide information on planet formation and the reasons for the diversity of planetary systems. The key to understanding planet formation is the study of dust evolution from small grains to pebbles. Smaller grains are well-coupled to the gas dynamics, and their distribution is significantly extended above the disk midplane. Larger grains settle much faster and are efficiently formed only in the midplane. By combining near-infrared polarized light and millimeter observations, it is possible to constrain the spatial distribution of both the small and large grains. We aim to construct detailed models of the size distribution and vertical/radial structure of the dust particles in protoplanetary disks based on observational data. In particular, we are interested in recovering the dust distribution in the IM Lup protoplanetary disk. We create a physical model for the dust distribution of protoplanetary disks and simulate the radiative transfer of the millimeter continuum and the near-infrared polarized radiation. Using a Markov chain Monte Carlo method, we compare the derived images to the observations available for the IM Lup disk to constrain the best physical model for IM Lup and to recover the vertical grain size distribution. The millimeter and near-infrared emission tightly constrain the dust mass and grain size distribution of our model. We find size segregation in the dust distribution, with millimeter-sized grains in the disk midplane. These grains are efficiently formed in the disk, possibly by sedimentation-driven coagulation, in accord with the short settling timescales predicted by our model. This also suggests a high dust-to-gas ratio at smaller radii in the midplane, possibly triggering streaming instabilities and planetesimal formation in the inner disk.
△ Less
Submitted 2 December, 2022;
originally announced December 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.
-
A Neural Network Subgrid Model of the Early Stages of Planet Formation
Authors:
Thomas Pfeil,
Miles Cranmer,
Shirley Ho,
Philip J. Armitage,
Tilman Birnstiel,
Hubert Klahr
Abstract:
Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation.…
▽ More
Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation. Numerical simulations which include the relevant physical effects are complex and computationally expensive. Here, we present a fast and accurate learned effective model for dust coagulation, trained on data from high resolution numerical coagulation simulations. Our model captures details of the dust coagulation process that were so far not tractable with other dust coagulation prescriptions with similar computational efficiency.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
The impact of dynamic pressure bumps on the observational properties of protoplanetary disks
Authors:
Jochen Stadler,
Matías Gárate,
Paola Pinilla,
Christian Lenz,
Cornelis P. Dullemond,
Til Birnstiel,
Sebastian M. Stammler
Abstract:
Over the last years, large (sub-)millimetre surveys of protoplanetary disks have well constrained the demographics of disks, such as their millimetre luminosities, spectral indices, and disk radii. Additionally, several high-resolution observations have revealed an abundance of substructures in the disks dust continuum. The most prominent are ring like structures, likely due to pressure bumps trap…
▽ More
Over the last years, large (sub-)millimetre surveys of protoplanetary disks have well constrained the demographics of disks, such as their millimetre luminosities, spectral indices, and disk radii. Additionally, several high-resolution observations have revealed an abundance of substructures in the disks dust continuum. The most prominent are ring like structures, likely due to pressure bumps trapping dust particles. The origins and characteristics of these bumps, nevertheless, need to be further investigated. The purpose of this work is to study how dynamic pressure bumps affect observational properties of protoplanetary disks. We further aim to differentiate between the planetary- versus zonal flow-origin of pressure bumps. We perform one-dimensional gas and dust evolution simulations, setting up models with varying pressure bump features. We subsequently run radiative transfer calculations to obtain synthetic images and the different quantities of observations. We find that the outermost pressure bump determines the disks dust size across different millimetre wavelengths. Our modelled dust traps need to form early (< 0.1 Myr), fast (on viscous timescales), and must be long lived (> Myr) to obtain the observed high millimetre luminosities and low spectral indices of disks. While the planetary bump models can reproduce these observables irrespectively of the opacity prescription, the highest opacities are needed for the zonal flow bump model to be in line with observations. Our findings favour the planetary- over the zonal flow-origin of pressure bumps and support the idea that planet formation already occurs in early class 0-1 stages of circumstellar disks. The determination of the disks effective size through its outermost pressure bump also delivers a possible answer to why disks in recent low-resolution surveys appear to have the same sizes across different millimetre wavelengths.
△ Less
Submitted 16 September, 2022;
originally announced September 2022.
-
Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science
Authors:
Ian Stotesbury,
Billy Edwards,
Jean-Francois Lavigne,
Vasco Pesquita,
James J. Veilleux,
Philip Windred,
Ahmed Al-Refaie,
Lawrence Bradley,
Sushuang Ma,
Giorgio Savini,
Giovanna Tinetti,
Til Birnstiel,
Sally Dodson-Robinson,
Barbara Ercolano,
Dax Feliz,
Scott Gaudi,
Nina Hernitschek,
Daniel Holdsworth,
Ing-Guey Jiang,
Matt Griffin,
Nataliea Lowson,
Karan Molaverdikhani,
Hilding Neilson,
Caprice Phillips,
Thomas Preibisch
, et al. (13 additional authors not shown)
Abstract:
With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the tel…
▽ More
With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey.
△ Less
Submitted 7 September, 2022;
originally announced September 2022.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Scattered light detection of a possible disk wind in RY Tau
Authors:
P. -G. Valegård,
C. Ginski,
C. Dominik,
J. Bae,
M. Benisty,
T. Birnstiel,
S. Facchini,
A. Garufi,
M. Hogerheijde,
R. G. van Holstein,
M. Langlois,
C. F. Manara,
P. Pinilla,
Ch. Rab,
Á. Ribas,
L. B. F. M. Waters,
J. Williams
Abstract:
Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-i…
▽ More
Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-infrared but not seen at optical wavelengths. In the near-infrared, dark wedges that separates the horns from the disk, indicating we may see the scattered light from a disk wind. We use archived ALMA and SPHERE/ZIMPOL I-band observations combined with newly acquired SPEHRE/IRDIS H-band observations and available literature to build a simple geometric model of the RY Tau disk and disk wind. We use Monte Carlo radiative transfer modelling \textit{MCMax3D} to create comparable synthetic observations that test the effect of a dusty wind on the optical effect in the observations. We constrain the grain size and dust mass needed in the disk wind to reproduce the effect from the observations. A model geometrically reminiscent of a dusty disk wind with small micron to sub-micron size grains elevated above the disk can reproduce the optical effect seen in the observations. The mass in the obscuring component of the wind has been constrained to $1\times10^{-9} M_{\odot} \leq M \leq 5\times10^{-8} M_{\odot}$ which corresponds to a lower limit mass loss rate in the wind of about $\sim 1\times10^{-8}M_{\odot}\mathrm{yr}^{-1}$. While an illuminate dust cavity cannot be ruled out without measurements of the gas velocity, we argue that a magnetically launched disk wind is the most likely scenario.
△ Less
Submitted 5 October, 2022; v1 submitted 5 September, 2022;
originally announced September 2022.
-
Toward a Population Synthesis of Disks and Planets I. Evolution of Dust with Entrainment in Winds and Radiation Pressure
Authors:
Remo Burn,
Alexandre Emsenhuber,
Jesse Weder,
Oliver Völkel,
Hubert Klahr,
Til Birnstiel,
Barbara Ercolano,
Christoph Mordasini
Abstract:
Millimeter astronomy provides valuable information on the birthplaces of planetary systems. In order to compare theoretical models with observations, the dust component has to be carefully calculated. Here, we aim to study the effects of dust entrainment in photoevaporative winds and the ejection and drag of dust due to effects caused by radiation from the central star. We improved and extended th…
▽ More
Millimeter astronomy provides valuable information on the birthplaces of planetary systems. In order to compare theoretical models with observations, the dust component has to be carefully calculated. Here, we aim to study the effects of dust entrainment in photoevaporative winds and the ejection and drag of dust due to effects caused by radiation from the central star. We improved and extended the existing implementation of a two-population dust and pebble description in the global Bern/Heidelberg planet formation and evolution model. Modern prescriptions for photoevaporative winds were used and we account for settling and advection of dust when calculating entrainment rates. In order to prepare for future population studies with varying conditions, we explore a wide range of disk-, photoevaporation-, and dust-parameters. We find that if dust can grow to pebble sizes, that is, if they are resistant to fragmentation or turbulence is weak, drift dominates and the entrained mass is small but larger than under the assumption of no vertical advection of grains with the gas flow. For the case of fragile dust shattering at velocities of 1 m/s - as indicated in laboratory experiments -, an order of magnitude more dust is entrained which becomes the main dust removal process. Radiation pressure effects disperse massive, dusty disks on timescales of a few 100 Myr. These results highlight the importance of dust entrainment in winds as a solid mass removal process. Furthermore, this model extension lies the basis for future statistical studies of planet formation in their birth environment.
△ Less
Submitted 18 July, 2022;
originally announced July 2022.
-
DustPy: A Python Package for Dust Evolution in Protoplanetary Disks
Authors:
Sebastian Markus Stammler,
Tilman Birnstiel
Abstract:
Many processes during the evolution of protoplanetary disks and during planet formation are highly sensitive to the sizes of dust particles that are present in the disk: The efficiency of dust accretion in the disk and volatile transport on dust particles, gravoturbulent instabilities leading to the formation of planetesimals, or the accretion of pebbles onto large planetary embryos to form giant…
▽ More
Many processes during the evolution of protoplanetary disks and during planet formation are highly sensitive to the sizes of dust particles that are present in the disk: The efficiency of dust accretion in the disk and volatile transport on dust particles, gravoturbulent instabilities leading to the formation of planetesimals, or the accretion of pebbles onto large planetary embryos to form giant planets are typical examples of processes that depend on the sizes of the dust particles involved. Furthermore, radiative properties like absorption or scattering opacities depend on the particle sizes. To interpret observations of dust in protoplanetary disks, a proper estimate of the dust particle sizes is needed.
We present DustPy - A Python package to simulate dust evolution in protoplanetary disks. DustPy solves gas and dust transport including viscous advection and diffusion as well as collisional growth of dust particles. DustPy is written with a modular concept, such that every aspect of the model can be easily modified or extended to allow for a multitude of research opportunities.
△ Less
Submitted 2 August, 2022; v1 submitted 1 July, 2022;
originally announced July 2022.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A Panchromatic View of DO Tau's Complex Kilo-au Environment
Authors:
Jane Huang,
Christian Ginski,
Myriam Benisty,
Bin Ren,
Alexander J. Bohn,
Élodie Choquet,
Karin I. Öberg,
Álvaro Ribas,
Jaehan Bae,
Edwin A. Bergin,
Til Birnstiel,
Yann Boehler,
Stefano Facchini,
Daniel Harsono,
Michiel Hogerheijde,
Feng Long,
Carlo F. Manara,
François Ménard,
Paola Pinilla,
Christophe Pinte,
Christian Rab,
Jonathan P. Williams,
Alice Zurlo
Abstract:
While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT…
▽ More
While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales.
△ Less
Submitted 8 May, 2022; v1 submitted 4 April, 2022;
originally announced April 2022.
-
Setting the Stage for Planet Formation: Measurements and Implications of the Fundamental Disk Properties
Authors:
A. Miotello,
I. Kamp,
T. Birnstiel,
L. I. Cleeves,
A. Kataoka
Abstract:
The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks,…
▽ More
The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks, as well as their spatial distribution, to better understand planet formation. Since PPVI, the field has seen tremendous improvements in observational capabilities, enabling both surveys of large samples of disks and high resolution imaging studies of a few bright disks. Improvements in data quality and sample size have, however, opened up many fundamental questions about properties such as the mass budget of disks, its spatial distribution, and its radial extent. Moreover, the vertical structure of disks has been studied in greater detail with spatially resolved observations, providing new insights on vertical layering and temperature stratification, yet also bringing rise to questions about other properties, such as material transport and viscosity. Each one of these properties - disk mass, surface density distribution, outer radius, vertical extent, temperature structure, and transport - is of fundamental interest as they collectively set the stage for disk evolution and corresponding planet formation theories. In this chapter, we will review our understanding of the fundamental properties of disks including the relevant observational techniques to probe their nature, modelling methods, and the respective caveats. Finally, we discuss the implications for theories of disk evolution and planet formation underlining what new questions have since arisen as our observational facilities have improved.
△ Less
Submitted 18 March, 2022;
originally announced March 2022.
-
Mineral Snowflakes on Exoplanets and Brown Dwarfs: Coagulation and Fragmentation of Cloud Particles with {\sc HyLandS}
Authors:
Dominic Samra,
Christiane Helling,
Tilman Birnstiel
Abstract:
Brown dwarfs and exoplanets provide unique atmospheric regimes that hold information about their formation routes and evolutionary states. Modelling mineral cloud particle formation is key to prepare for missions and instruments like CRIRES+, JWST and ARIEL as well as possible polarimetry missions like {\sc PolStar}. The aim is to support more detailed observations that demand greater understandin…
▽ More
Brown dwarfs and exoplanets provide unique atmospheric regimes that hold information about their formation routes and evolutionary states. Modelling mineral cloud particle formation is key to prepare for missions and instruments like CRIRES+, JWST and ARIEL as well as possible polarimetry missions like {\sc PolStar}. The aim is to support more detailed observations that demand greater understanding of microphysical cloud processes. We extend our kinetic cloud formation model that treats nucleation, condensation, evaporation and settling of mixed material cloud particles to consistently model cloud particle-particle collisions. The new hybrid code, {\sc HyLandS}, is applied to a grid of {\sc Drift-Phoenix} (T, p)-profiles. Effective medium theory and Mie theory are used to investigate the optical properties. Turbulence is the main driving process of collisions, with collisions becoming the dominant process at the cloud base ($p>10^{-4}\,{\rm bar}$). Collisions produce one of three outcomes: fragmenting atmospheres ($\log_{10}(g)=3$), coagulating atmospheres ($\log_{10}(g)=5$, $T_{\rm eff} \leq 1800\, {\rm K}$) and condensational growth dominated atmospheres ($\log_{10}(g\,)=5$, $T_{\rm eff} > 1800\, {\rm K}$). Cloud particle opacity slope at optical wavelengths (HST) is increased with fragmentation, as are the silicate features at mid-infrared wavelengths. The hybrid moment-bin method {\sc HyLandS} demonstrates the feasibility of combining a moment and a bin method whilst assuring element conservation. It provides a powerful and fast tool for capturing general trends of particle collisions, consistently with other microphysical processes. Collisions are important in exoplanet and brown dwarf atmospheres but cannot be assumed to be hit-and-stick only. The spectral effects of collisions complicates inferences of cloud particle size and material composition from observational data.
△ Less
Submitted 14 March, 2022;
originally announced March 2022.
-
A large population study of protoplanetary disks: Explaining the millimeter size-luminosity relation with or without sub-structure
Authors:
Apostolos Zormpas,
Tilman Birnstiel,
Giovanni P. Rosotti,
Sean M. Andrews
Abstract:
Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolut…
▽ More
Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $α$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure.
△ Less
Submitted 23 May, 2022; v1 submitted 2 February, 2022;
originally announced February 2022.
-
Dust entrainment in photoevaporative winds: Synthetic observations of transition disks
Authors:
R. Franz,
G. Picogna,
B. Ercolano,
S. Casassus,
T. Birnstiel,
Ch. Rab,
S. Pérez
Abstract:
X-ray- and extreme-ultraviolet- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may strongly impact disk evolution, affecting both gas and dust distributions. We compute dust densities for the wind regions of XEUV-irradiated transition disks with gap sizes of 20 and 30 AU, and determine whether they can be observed at wavelengths…
▽ More
X-ray- and extreme-ultraviolet- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may strongly impact disk evolution, affecting both gas and dust distributions. We compute dust densities for the wind regions of XEUV-irradiated transition disks with gap sizes of 20 and 30 AU, and determine whether they can be observed at wavelengths $0.7 \lesssim λ_\mathrm{obs} [μ\mathrm{m}] \lesssim 1.8$ in scattered and polarised light with current instrumentation. For an XEUV-driven outflow around a $M_* = 0.7 \mathrm{M}_\odot$ T-Tauri star with $L_X = 2 \cdot 10^{30} \mathrm{erg/s}$, we find dust mass-loss rates $\dot{M}_\mathrm{dust} \lesssim 2.0 \cdot 10^{-3} \dot{M}_\mathrm{gas}$, and if we invoke vertical settling, the outflow is quite collimated. The synthesised images exhibit a distinct chimney-like structure. The relative intensity of these chimneys is low, but under optimal conditions, their detection may still be feasible with current instrumentation such as JWST NIRCam and SPHERE IRDIS.
△ Less
Submitted 23 April, 2022; v1 submitted 28 January, 2022;
originally announced January 2022.
-
Dust entrainment in photoevaporative winds: Densities and imaging
Authors:
R. Franz,
B. Ercolano,
S. Casassus,
G. Picogna,
T. Birnstiel,
S. Pérez,
Ch. Rab,
A. Sharma
Abstract:
X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We constrain the dust densities in a typical XEUV-driven outflow, and determine whether these winds can be observed at $μ\mathrm{m}$-wavelengths in scattered and polarised light. For an XEUV-driven outflow ar…
▽ More
X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We constrain the dust densities in a typical XEUV-driven outflow, and determine whether these winds can be observed at $μ\mathrm{m}$-wavelengths in scattered and polarised light. For an XEUV-driven outflow around a $M_* = 0.7\,\mathrm{M}_\odot$ T-Tauri star with $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$, we find a dust mass-loss rate $\dot{M}_\mathrm{dust} \lesssim 4.1 \cdot 10^{-11}\,\mathrm{M_\odot / yr}$ for an optimistic estimate of dust densities in the wind (compared to $\dot{M}_\mathrm{gas} \approx 3.7 \cdot 10^{-8}\,\mathrm{M_\odot / yr}$). The synthesised scattered-light images suggest a distinct chimney structure emerging at intensities $I/I_{\max} < 10^{-4.5}$ ($10^{-3.5}$) at $λ_\mathrm{obs} = 1.6$ ($0.4$) $μ\mathrm{m}$, while the features in the polarised-light images are even fainter. Observations synthesised from our model do not exhibit clear features for SPHERE IRDIS, but show a faint wind signature for JWST NIRCam under optimal conditions. In conclusion, unambiguous detections of photoevaporative XEUV winds launched from primordial disks are at least challenging with current instrumentation; this provides a possible explanation as to why disk winds are not routinely detected in scattered or polarised light.
△ Less
Submitted 14 December, 2021; v1 submitted 20 October, 2021;
originally announced October 2021.
-
Large gaps and high accretion rates in photoevaporative transition disks with a dead zone
Authors:
Matías Gárate,
Timmy N. Delage,
Jochen Stadler,
Paola Pinilla,
Til Birnstiel,
Sebastian M. Stammler,
Giovanni Picogna,
Barbara Ercolano,
Raphael Franz,
Christian Lenz
Abstract:
Observations of young stars hosting transition disks show that several of them have high accretion rates, despite their disks presenting extended cavities in their dust component. This represents a challenge for theoretical models, which struggle to reproduce both features. We explore if a disk evolution model, including a dead zone and disk dispersal by X-ray photoevaporation, can explain the hig…
▽ More
Observations of young stars hosting transition disks show that several of them have high accretion rates, despite their disks presenting extended cavities in their dust component. This represents a challenge for theoretical models, which struggle to reproduce both features. We explore if a disk evolution model, including a dead zone and disk dispersal by X-ray photoevaporation, can explain the high accretion rates and large gaps (or cavities) measured in transition disks. We implement a dead zone turbulence profile and a photoevaporative mass loss profile into numerical simulations of gas and dust. We perform a population synthesis study of the gas component, and obtain synthetic images and SED of the dust component through radiative transfer calculations. This model results in long lived inner disks and fast dispersing outer disks, that can reproduce both the accretion rates and gap sizes observed in transition disks. For a dead zone of turbulence $α_{dz} = 10^{-4}$ and extent $r_{dz}$ = 10 AU, our population synthesis study shows that $63\%$ of our transition disks are accreting with $\dot{M}_g > 10^{-11} M_\odot/yr$ after opening a gap. Among those accreting transition disks, half display accretion rates higher than $5\times10^{-10} M_\odot/yr$ . The dust component in these disks is distributed in two regions: in a compact inner disk inside the dead zone, and in a ring at the outer edge of the photoevaporative gap, which can be located between 20 AU and 100 AU. Our radiative transfer calculations show that the disk displays an inner disk and an outer ring in the millimeter continuum, a feature observed in some transition disks. A disk model considering X-ray photoevaporative dispersal in combination with dead zones can explain several of the observed properties in transition disks including: the high accretion rates, the large gaps, and long-lived inner disks at mm-emission.
△ Less
Submitted 18 October, 2021;
originally announced October 2021.
-
Mass determination of protoplanetary disks from dust evolution
Authors:
Riccardo Franceschi,
Tilman Birnstiel,
Thomas Henning,
Paola Pinilla,
Dmitry Semenov,
Apostolos Zormpas
Abstract:
The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location…
▽ More
The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location of the "dust lines". This is a new concept in the field, linking the disk size at different continuum wavelengths with the radial distribution of grain populations of different sizes. We aim to use a dust evolution model to test the dependence of the dust line location on disk gas mass. In particular, we are interested in the reliability of the method for disks showing radial substructures, as recent high-resolution observations revealed. Our models show that the determination of the dust line location is a promising approach to the mass estimate of protoplanetay disks, but the exact relation between the dust line location and disk mass depends on the structure of the particular disk. We calibrated the relation for disks without evidence of radial structures, while for more complex structures we ran a simple dust evolution model. However, this method fails when there is evidence of strong dust traps. It is possible to reveal when dust evolution is dominated by traps, providing the necessary information for when the method should be applied with caution.
△ Less
Submitted 18 October, 2021;
originally announced October 2021.
-
The formation of wide exoKuiper belts from migrating dust traps
Authors:
E. Miller,
S. Marino,
S. M. Stammler,
P. Pinilla,
C. Lenz,
T. Birnstiel,
Th. Henning
Abstract:
The question of what determines the width of Kuiper belt analogues (exoKuiper belts) is an open one. If solved, this understanding would provide valuable insights into the architecture, dynamics, and formation of exoplanetary systems. Recent observations by ALMA have revealed an apparent paradox in this field, the presence of radially narrow belts in protoplanetary discs that are likely the birthp…
▽ More
The question of what determines the width of Kuiper belt analogues (exoKuiper belts) is an open one. If solved, this understanding would provide valuable insights into the architecture, dynamics, and formation of exoplanetary systems. Recent observations by ALMA have revealed an apparent paradox in this field, the presence of radially narrow belts in protoplanetary discs that are likely the birthplaces of planetesimals, and exoKuiper belts nearly four times as wide in mature systems. If the parent planetesimals of this type of debris disc indeed form in these narrow protoplanetary rings via streaming instability where dust is trapped, we propose that this width dichotomy could naturally arise if these dust traps form planetesimals whilst migrating radially, e.g. as caused by a migrating planet. Using the dust evolution software DustPy, we find that if the initial protoplanetary disc and trap conditions favour planetesimal formation, dust can still effectively accumulate and form planetesimals as the trap moves. This leads to a positive correlation between the inward radial speed and final planetesimal belt width, forming belts up to $\sim$100 au over 10 Myr of evolution. We show that although planetesimal formation is most efficient in low viscosity ($α= 10^{-4}$) discs with steep dust traps to trigger the streaming instability, the large widths of most observed planetesimal belts constrain $α$ to values $\geq4\times 10^{-4}$ at tens of au, otherwise the traps cannot migrate far enough. Additionally, the large spread in the widths and radii of exoKuiper belts could be due to different trap migration speeds (or protoplanetary disc lifetimes) and different starting locations, respectively. Our work serves as a first step to link exoKuiper belts and rings in protoplanetary discs.
△ Less
Submitted 8 October, 2021;
originally announced October 2021.
-
Thermal Wave Instability as an Origin of Gap and Ring Structures in Protoplanetary Disks
Authors:
Takahiro Ueda,
Mario Flock,
Tilman Birnstiel
Abstract:
Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk…
▽ More
Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk is optically thick enough for stellar light, i.e., small-grain-to-gas mass ratio of $\gtrsim0.0001$. The mid-plane temperature varies as the waves propagate and hence gap and ring structures can be seen in both millimeter and infrared emission. The millimeter substructures can be observed even if the disk is fully optically thick since it is induced by the temperature variation, while density-induced substructures would disappear in the optically thick regime. The fractional separation between TWI-induced ring and gap is $Δr/r \sim$ 0.2-0.4 at $\sim$ 10-50 au, which is comparable to those found by ALMA. Due to the temperature variation, snow lines of volatile species move radially and multiple snow lines are observed even for a single species. The wave propagation velocity is as fast as $\sim$ 0.6 ${\rm au~yr^{-1}}$, which can be potentially detected with a multi-epoch observation with a time separation of a few years.
△ Less
Submitted 28 May, 2021;
originally announced May 2021.
-
First 3-D grid-based gas-dust simulations of circumstellar disks with an embedded planet
Authors:
Fabian Binkert,
Judit Szulágyi,
Til Birnstiel
Abstract:
Substructures are ubiquitous in high resolution (sub-)millimeter continuum observations of circumstellar disks. They are possibly caused by forming planets embedded in the disk. To investigate the relation between observed substructures and young planets, we perform novel three-dimensional two-fluid (gas+1-mm-dust) hydrodynamic simulations of circumstellar disks with embedded planets (Neptune-, Sa…
▽ More
Substructures are ubiquitous in high resolution (sub-)millimeter continuum observations of circumstellar disks. They are possibly caused by forming planets embedded in the disk. To investigate the relation between observed substructures and young planets, we perform novel three-dimensional two-fluid (gas+1-mm-dust) hydrodynamic simulations of circumstellar disks with embedded planets (Neptune-, Saturn-, Jupiter-, 5 Jupiter-mass) at different orbital distances from the star (5.2AU, 30AU, 50AU). We turn these simulations into synthetic (sub-)millimeter ALMA images. We find that all but the Neptune-mass planet open annular gaps in both the gas and the dust component of the disk. We find that the temporal evolution of the dust density distribution is distinctly different of the gas'. For example, the planets cause significant vertical stirring of the dust in the circumstellar disk which opposes the vertical settling. This creates a thicker dust disk than disks without a planet. We find that this effect greatly influences the dust masses derived from the synthetic ALMA images. Comparing the dust disk masses in the 3D simulations and the ones derived from the 2D ALMA synthetic images, we find the former to be a factor of a few (up to 10) larger, pointing to that real disks might be significantly more massive than previously thought based on ALMA continuum images using the optically thin assumption and equation. Finally, we analyze the synthetic ALMA images and provide an empirical relationship between the planet mass and the width of the gap in the ALMA images including the effects of the beam size.
△ Less
Submitted 16 July, 2021; v1 submitted 18 March, 2021;
originally announced March 2021.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Late infall causing disk misalignment and dynamic structures in SU Aur
Authors:
C. Ginski,
S. Facchini,
J. Huang,
M. Benisty,
D. Vaendel,
L. Stapper,
C. Dominik,
J. Bae,
F. Menard,
G. Muro-Arena,
M. Hogerheijde,
M. McClure,
R. G. van Holstein,
T. Birnstiel,
Y. Boehler,
A. Bohn,
M. Flock,
E. E. Mamajek,
C. F. Manara,
P. Pinilla,
C. Pinte,
A. Ribas
Abstract:
Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resol…
▽ More
Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems.
△ Less
Submitted 17 February, 2021;
originally announced February 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.
-
A search for companions via direct imaging in the DSHARP planet-forming disks
Authors:
Sebastián Jorquera,
Laura M. Pérez,
Gaël Chauvin,
Myriam Benisty,
Zhaohuan Zhu,
Andrea Isella,
Jane Huang,
Luca Ricci,
Sean M. Andrews,
Shangjia Zhang,
John Carpenter,
Nicolás T. Kurtovic,
Tilman Birnstiel
Abstract:
The "Disk Substructures at High Angular Resolution Project" (DSHARP) has revealed an abundance and ubiquity of rings and gaps over a large sample of young planet-forming disks, which are hypothesised to be induced by the presence of forming planets. In this context, we present the first attempt to directly image these young companions for 10 of the DSHARP disks, by using NaCo/VLT high contrast obs…
▽ More
The "Disk Substructures at High Angular Resolution Project" (DSHARP) has revealed an abundance and ubiquity of rings and gaps over a large sample of young planet-forming disks, which are hypothesised to be induced by the presence of forming planets. In this context, we present the first attempt to directly image these young companions for 10 of the DSHARP disks, by using NaCo/VLT high contrast observations in L'-band instrument and angular differential imaging techniques. We report the detection of a point-like source candidate at 1.1" (174.9 au) for RU Lup, and at 0.42" (55 AU) for Elias 24. In the case of RU Lup, the proper motion of the candidate is consistent with a stationary background contaminant, based on the astrometry derived from our observations and available archival data. For Elias 24 the point-like source candidate is located in one of the disk gaps at 55 AU. Assuming it is a planetary companion, our analysis suggest a mass ranging from $0.5 M_J$ up to $5 M_J$, depending on the presence of a circumplanetary disk and its contribution to the luminosity of the system. However, no clear confirmation is obtained at this stage, and follow-up observations are mandatory to verify if the proposed source is physical, comoving with the stellar host, and associated with a young massive planet sculpting the gap observed at 55\,AU. For all the remaining systems, the lack of detections suggests the presence of planetary companions with masses lower than $5M_J$, based on our derived mass detection limits. This is consistent with predictions of both hydrodynamical simulations and kinematical signatures on the disk, and allows us to set upper limits on the presence of massive planets in these young disks.
△ Less
Submitted 31 March, 2021; v1 submitted 18 December, 2020;
originally announced December 2020.
-
On the diversity of asymmetries in gapped protoplanetary disks
Authors:
Nienke van der Marel,
Til Birnstiel,
Antonio Garufi,
Enrico Ragusa,
Valentin Christiaens,
Daniel Price,
Steph Sallum,
Dhruv Muley,
Logan Francis,
Ruobing Dong
Abstract:
Protoplanetary disks with large inner dust cavities are thought to host massive planetary or substellar companions. These disks show asymmetries and rings in the millimeter continuum, caused by dust trapping in pressure bumps, and potentially vortices or horseshoes. The origin of the asymmetries and their diversity remains unclear. We present a comprehensive study of 16 disks for which the gas sur…
▽ More
Protoplanetary disks with large inner dust cavities are thought to host massive planetary or substellar companions. These disks show asymmetries and rings in the millimeter continuum, caused by dust trapping in pressure bumps, and potentially vortices or horseshoes. The origin of the asymmetries and their diversity remains unclear. We present a comprehensive study of 16 disks for which the gas surface density profile has been constrained by CO isotopologue data. We compare the azimuthal extents of the dust continuum profiles with the local gas surface density in each disk, and find that the asymmetries correspond to higher Stokes numbers or low gas surface density. We discuss which asymmetric structures can be explained by a horseshoe, a vortex or spiral density waves. Second, we reassess the gas gap radii from the $^{13}$CO maps, which are about a factor 2 smaller than the dust ring radii, suggesting that companions in these disks are in the brown dwarf mass regime ($\sim 15-50 M_{\rm Jup}$) or in the Super-Jovian mass regime ($\sim 3-15 M_{\rm Jup}$) on eccentric orbits. This is consistent with the estimates from contrast curves on companion mass limits. These curves rule out (sub)stellar companions ($q>$0.05) for the majority of the sample at the gap location, but it remains possible at even smaller radii. Third, we find that spiral arms in scattered light images are primarily detected around high luminosity stars with disks with wide gaps, which can be understood by the dependence of the spiral arm pitch angle on disk temperature and companion mass.
△ Less
Submitted 22 October, 2020; v1 submitted 20 October, 2020;
originally announced October 2020.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): A close low mass companion to ET Cha
Authors:
C. Ginski,
F. Ménard,
Ch. Rab,
E. E. Mamajek,
R. G. van Holstein,
M. Benisty,
C. F. Manara,
R. Asensio Torres,
A. Bohn,
T. Birnstiel,
P. Delorme,
S. Facchini,
A. Garufi,
R. Gratton,
M. Hogerheijde,
J. Huang,
M. Kenworthy,
M. Langlois,
P. Pinilla,
C. Pinte,
Á. Ribas,
G. Rosotti,
T. O. B. Schmidt,
M. van den Ancker,
Z. Wahhaj
, et al. (3 additional authors not shown)
Abstract:
To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed.…
▽ More
To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Constraining the parameter space for the Solar Nebula
Authors:
Christian T. Lenz,
Hubert Klahr,
Tilman Birnstiel,
Katherine Kretke,
Sebastian Stammler
Abstract:
If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)pl…
▽ More
If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)planets, once the gas disk is gone. Protoplanetary disks provide a plethora of different parameters to explore. The question is whether this parameter space can be constrained, allowing simulations to reproduce the Solar System. Models and observations of planet formation provide constraints on the initial planetesimal mass in certain regions of the Solar Nebula. By making use of pebble flux-regulated planetesimal formation, we perform a parameter study with nine different disk parameters like the initial disk mass, initial disk size, initial dust-to-gas ratio, turbulence level, and more. We find that the distribution of mass in planetesimals in the disk depends on the planetesimal formation timescale and the pebbles' drift timescale. Multiple disk parameters can influence pebble properties and thus planetesimal formation. However, it is still possible to draw some conclusions on potential parameter ranges. Pebble flux-regulated planetesimal formation seems to be very robust, allowing simulations with a wide range of parameters to meet the initial planetesimal constraints for the Solar Nebula. I.e., it does not require a lot of fine tuning.
△ Less
Submitted 15 June, 2020;
originally announced June 2020.
-
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.
-
A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk
Authors:
Jane Huang,
Sean M. Andrews,
Cornelis P. Dullemond,
Karin I. Oberg,
Chunhua Qi,
Zhaohuan Zhu,
Tilman Birnstiel,
John M. Carpenter,
Andrea Isella,
Enrique Macias,
Melissa K. McClure,
Laura M. Perez,
Richard Teague,
David J. Wilner,
Shangjia Zhang
Abstract:
The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observatio…
▽ More
The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $α$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."
△ Less
Submitted 29 January, 2020;
originally announced January 2020.
-
Dust entrainment in photoevaporative winds: The impact of X-rays
Authors:
R. Franz,
G. Picogna,
B. Ercolano,
T. Birnstiel
Abstract:
X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We investigate the dust entrainment in XEUV-driven photoevaporative winds and compare our results to existing MHD and EUV-only models. For an X-ray luminosity of $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$ emitte…
▽ More
X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We investigate the dust entrainment in XEUV-driven photoevaporative winds and compare our results to existing MHD and EUV-only models. For an X-ray luminosity of $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$ emitted by a $M_* = 0.7\,\mathrm{M}_\odot$ star, corresponding to a wind mass-loss rate of $\dot{M}_\mathrm{w} \simeq 2.6 \cdot 10^{-8} \,\mathrm{M_\odot/yr}$, we find dust entrainment for sizes $a_0 \lesssim 11\,μ$m ($9\,μ$m) from the inner $25\,$AU ($120\,$AU). This is an enhancement over dust entrainment in less vigorous EUV-driven winds with $\dot{M}_\mathrm{w} \simeq 10^{-10}\,\mathrm{M_\odot/yr}$. Our numerical model also shows deviations of dust grain trajectories from the gas streamlines even for $μ$m-sized particles. In addition, we find a correlation between the size of the entrained grains and the maximum height they reach in the outflow.
△ Less
Submitted 28 January, 2020;
originally announced January 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.