-
Spirals, rings, and vortices shaped by shadows in protoplanetary disks: from radiative hydrodynamical simulations to observable signatures
Authors:
Alexandros Ziampras,
Cornelis P. Dullemond,
Tilman Birnstiel,
Myriam Benisty,
Richard P. Nelson
Abstract:
Numerous protoplanetary disks exhibit shadows in scattered light observations. These shadows are typically cast by misaligned inner disks and are associated with observable structures in the outer disk such as bright arcs and spirals. Investigating the dynamics of the shadowed outer disk is therefore essential in understanding the formation and evolution of these structures. We carry out twodimens…
▽ More
Numerous protoplanetary disks exhibit shadows in scattered light observations. These shadows are typically cast by misaligned inner disks and are associated with observable structures in the outer disk such as bright arcs and spirals. Investigating the dynamics of the shadowed outer disk is therefore essential in understanding the formation and evolution of these structures. We carry out twodimensional radiation hydrodynamics simulations that include radiative diffusion and dust-gas dynamics to study the formation of substructure in shadowed disks. We find that spiral arms are launched at the edge of each shadow, permeating the entire disk. The local dissipation of these spirals results in an angular momentum flux, opening multiple gaps and leading to a series of concentric, regularly-spaced rings We find that ring formation is favored in weakly turbulent disks where dust growth is taking place. These conditions are met for typical class-II disks, in which bright rings should form well within a fraction of their lifetime (0.1-0.2 Myr). For hotter disks gap opening is more efficient, such that the gap edges quickly collapse into vortices that can appear as bright arcs in continuum emission before decaying into rings or merging into massive, long-lived structures. Synthetic observations show that these structures should be observable in scattered light and millimeter continuum emission, providing a new way to probe the presence of substructure in protoplanetary disks. Our results suggest that the formation of rings and gaps is a common process in shadowed disks, and can explain the rich radial substructure observed in several protoplanetary disks.
△ Less
Submitted 17 October, 2024;
originally announced October 2024.
-
Dusty substructures induced by planets in ALMA disks: how dust growth and dynamics changes the picture
Authors:
Alexandros Ziampras,
Prakruti Sudarshan,
Cornelis P. Dullemond,
Mario Flock,
Vittoria Berta,
Richard P. Nelson,
Andrea Mignone
Abstract:
Protoplanetary disks exhibit a rich variety of substructure in millimeter continuum emission, often attributed to unseen planets. As these planets carve gaps in the gas, dust particles can accumulate in the resulting pressure bumps, forming bright features in the dust continuum. We investigate the role of dust dynamics in the gap-opening process with 2D radiation hydrodynamics simulations of plane…
▽ More
Protoplanetary disks exhibit a rich variety of substructure in millimeter continuum emission, often attributed to unseen planets. As these planets carve gaps in the gas, dust particles can accumulate in the resulting pressure bumps, forming bright features in the dust continuum. We investigate the role of dust dynamics in the gap-opening process with 2D radiation hydrodynamics simulations of planet--disk interaction and a two-population dust component modeled as a pressureless fluid. We consider the opacity feedback and backreaction due to drag forces as mm grains accumulate in pressure bumps at different stages of dust growth. We find that dust dynamics can significantly affect the resulting substructure driven by the quasi-thermal-mass planet with $M_p/M_\star=10^{-4}$. Opacity feedback causes nonaxisymmetric features to become more compact in azimuth, whereas the drag-induced backreaction tends to dissolve nonaxisymmetries. For our fiducial model, this results in multiple concentric rings of dust rather than the expected vortices and corotating dust clumps found in models without dust feedback. A higher coagulation fraction disproportionately enhances the effect of dust opacity feedback, favoring the formation of crescents rather than rings. Our results suggest that turbulent diffusion is not always necessary to explain the rarity of observed nonaxisymmetric features, and that incorporating dust dynamics is vital for interpreting the observed substructure in protoplanetary disks. We also describe and test the implementation of the publicly-available dust fluid module in the PLUTO code.
△ Less
Submitted 23 September, 2024;
originally announced September 2024.
-
Formation and evolution of a protoplanetary disk: combining observations, simulations and cosmochemical constraints
Authors:
Alessandro Morbidelli,
Yves Marrocchi,
Adnan Ali Ahmad,
Asmita Bhandare,
Sebastien Charnoz,
Benoit Commercon,
Cornellis P. Dullemond,
Tristan Guillot,
Patrick Hennebelle,
Yueh-Ning Lee,
Francesco Lovascio,
Raphael Marschall,
Bernard Marty,
Anaelle Maury,
Okamoto Tamami
Abstract:
We present a plausible and coherent view of the evolution of the protosolar disk that is consistent with the cosmochemical constraints and compatible with observations of other protoplanetary disks and sophisticated numerical simulations. The evidence that high-temperature condensates, CAIs and AOAs, formed near the protosun before being transported to the outer disk can be explained by either an…
▽ More
We present a plausible and coherent view of the evolution of the protosolar disk that is consistent with the cosmochemical constraints and compatible with observations of other protoplanetary disks and sophisticated numerical simulations. The evidence that high-temperature condensates, CAIs and AOAs, formed near the protosun before being transported to the outer disk can be explained by either an early phase of vigorous radial spreading of the disk, or fast transport of these condensates from the vicinity of the protosun towards large disk radii via the protostellar outflow. The assumption that the material accreted towards the end of the infall phase was isotopically distinct allows us to explain the observed dichotomy in nucleosynthetic isotopic anomalies of meteorites and leads to intriguing predictions on the isotopic composition of refractory elements in comets. When the infall of material waned, the disk started to evolve as an accretion disk. Initially, dust drifted inwards, shrinking the radius of the dust component to ~ 45 au, probably about 1/2 of the width of the gas component. Then structures must have emerged, producing a series of pressure maxima in the disk which trapped the dust on My timescales. This allowed planetesimals to form at radically distinct times without changing significantly of isotopic properties. There was no late accretion of material onto the disk via streamers. The disk disappeared in ~5 Myr, as indicated by paleomagnetic data in meteorites. In conclusion, the evolution of the protosolar disk seems to have been quite typical in terms of size, lifetime, and dust behavior, suggesting that the peculiarities of the Solar system with respect to extrasolar planetary system probably originate from the chaotic nature of planet formation and not at the level of the parental disk.
△ Less
Submitted 10 September, 2024;
originally announced September 2024.
-
Warped Disk Evolution in Grid-Based Simulations
Authors:
C. N. Kimmig,
C. P. Dullemond
Abstract:
Observations show evidence that a significant fraction of protoplanetary disks contain warps. A warp in a disk evolves in time affecting the appearance of shadows and greatly influencing kinematic signatures. So far, many theoretical studies of warped disks have been conducted using Smoothed Particle Hydrodynamics (SPH) methods. In our approach, we use a grid-based method in spherical coordinates…
▽ More
Observations show evidence that a significant fraction of protoplanetary disks contain warps. A warp in a disk evolves in time affecting the appearance of shadows and greatly influencing kinematic signatures. So far, many theoretical studies of warped disks have been conducted using Smoothed Particle Hydrodynamics (SPH) methods. In our approach, we use a grid-based method in spherical coordinates which has notable advantages: the method allows for accurate modelling of low viscosity values and the resolution does not depend on density or mass of the disk, which allows surface structures to be resolved. We perform 3D simulations using FARGO3D to simulate the evolution of a warped disk and compare the results to one-dimensional models using a ring code. Additionally, we extensively investigate the applicability of grid-based methods to misaligned disks and test their dependency on grid resolution as well as disk viscosity. We find that grid-based simulations are capable of simulating disks not aligned to the grid geometry. Our three-dimensional simulation of a warped disk compares well to one-dimensional models in evolution of inclination. However, we find a twist which is not captured in 1D models. After thorough analysis we suspect this to be a physical effect possibly caused by non-linear effects neglected in the one-dimensional equations. Evaluating the internal dynamics, we find sloshing and breathing motions as predicted in local shearing box analysis. They can become supersonic, which may have consequences on kinematic observations of warped disks. Warped disks can be accurately modelled in 3D grid-based simulations when using reasonably good resolution, especially in the $θ$-direction. We find good agreement with the linear approximation of the sloshing motion which highlights the reliability of 1D models.
△ Less
Submitted 4 June, 2024;
originally announced June 2024.
-
Spatial distribution of crystalline silicates in protoplanetary disks: How to interpret mid-infrared observations
Authors:
Hyerin Jang,
L. B. F. M. Waters,
I. Kamp,
C. P. Dullemond
Abstract:
Crystalline silicates are an important tracer to the dust evolution in protoplanetary disks. In the inner disk, amorphous silicates are annealed by the high temperatures. These crystalline silicates are radially and vertically distributed in the disk. We aim to model the spatial distribution of crystalline silicate in the disk and its mid-IR spectra to study the effect on dust spectral features an…
▽ More
Crystalline silicates are an important tracer to the dust evolution in protoplanetary disks. In the inner disk, amorphous silicates are annealed by the high temperatures. These crystalline silicates are radially and vertically distributed in the disk. We aim to model the spatial distribution of crystalline silicate in the disk and its mid-IR spectra to study the effect on dust spectral features and to compare these to observations. We modeled a T-Tauri protoplanetary disk and defined the crystallization region from the crystallization and residence timescales. Radial mixing and drift were compared to find a vertically mixed region. We used the DISKLAB code to obtain the spatial distribution of the crystalline silicates, and MCMax code to model the mid-infrared spectrum. In our modeled disk, different grain sizes get crystallized in different regions in the disk. Crystallized dust in the disk surface is well mixed with the midplane due to vertical mixing and gets distributed to the outer disk by radial transport. Our model shows different contributions of the disk zones to the dust spectral features. Feature strengths change when varying the spatial distribution of crystalline dust. Our modeled spectra qualitatively agree with observations, but the modeled 10 $μ$m feature is strongly dominated by crystalline dust. Models with reduced crystallinity and depletion of small crystalline dust in the inner disk show a better match with observations. Mid-IR observations of the disk surface represent the radial distribution of small dust in the midplane and provide us with dust properties in the inner disk. The inner and outer disks contribute more to shorter and longer wavelength features, respectively. Amorphization, sublimation, and dust evolution have to be considered to match observations. This study could interpret the spectra of protoplanetary disks taken with the MIRI on board the JWST.
△ Less
Submitted 1 May, 2024;
originally announced May 2024.
-
Feasibility of detecting shadows in disks induced by infall
Authors:
A. Krieger,
M. Kuffmeier,
S. Reissl,
C. P. Dullemond,
C. Ginski,
S. Wolf
Abstract:
Observations performed with high-resolution imaging techniques revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner with respect to an outer disk. The cause of misalignment, however, is still debated. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which inv…
▽ More
Observations performed with high-resolution imaging techniques revealed the existence of shadows in circumstellar disks that can be explained by the misalignment of an inner with respect to an outer disk. The cause of misalignment, however, is still debated. In this study, we investigate the feasibility of observing shadows induced by one prominent scenario that may lead to misalignment, which involves the late infall of material onto a protostellar system. In particular, we use previously performed hydrodynamical simulations of such events, and generate flux maps in the visible, near-infrared, submillimeter, and millimeter wavelength range using Monte Carlo radiative transfer. Based on that, we derive synthetic observations of these systems performed with the instruments SPHERE/VLT and ALMA, which we use as a basis for our subsequent analysis. We find that near-infrared observations with SPHERE are particularly well suited for detecting shadows via direct imaging alongside other features such as gaps, arcs, and streamers. On the contrary, performing a shadow detection based on reconstructed ALMA observations is very challenging due to the high sensitivity that is required for this task. Thus, in cases that allow for a detection, sophisticated analyses may be needed, for instance by the utilization of carefully constructed azimuthal profiles, aiding the search for potentially shallow shadows. Lastly, we conclude that late infall-induced disk misalignment offers a plausible explanation for the emergence of shadows that are observed in various systems.
△ Less
Submitted 13 March, 2024;
originally announced March 2024.
-
FargoCPT: A 2D Multi-Physics Code for Simulating the Interaction of Disks with Stars, Planets and Particles
Authors:
Thomas Rometsch,
Lucas M. Jordan,
Tobias W. Moldenhauer,
Dennis Wehner,
Steven Rendon Restrepo,
Tobias W. A. Müller,
Giovanni Picogna,
Wilhelm Kley,
Cornelis P. Dullemond
Abstract:
Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, a…
▽ More
Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, an updated version of FARGO incorporating other previous enhancements to the code, to provide a simulation environment tailored to study interactions between stars, planets, and disks, ensuring accurate representation of planet systems, hydrodynamics, and dust dynamics with a focus on usability. Methods: The radiation-hydrodynamics part of FargoCPT uses a second-order upwind scheme in 2D polar coordinates supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. Integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapting to drag laws for the Stokes and Epstein regimes. Moreover, FargoCPT provides mechanisms to simulate accretion onto the stars and planets. Results: The code has been tested in practice by its use in various publications. Additionally, it comes with an automated test suite to test the physics modules. Conclusions: FargoCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in the ongoing exploration of planet formation.
△ Less
Submitted 29 January, 2024;
originally announced January 2024.
-
The SPHERE view of three interacting twin disc systems in polarised light
Authors:
Philipp Weber,
Sebastián Pérez,
Greta Guidi,
Nicolás T. Kurtovic,
Alice Zurlo,
Antonio Garufi,
Paola Pinilla,
Satoshi Mayama,
Rob G. van Holstein,
Cornelis P. Dullemond,
Nicolás Cuello,
David Principe,
Lucas Cieza,
Camilo González-Ruilova,
Julien Girard
Abstract:
Dense stellar environments as hosts of ongoing star formation increase the probability of gravitational encounters among stellar systems during the early stages of evolution. Stellar interaction may occur through non-recurring, hyperbolic or parabolic passages (a so-called 'fly-by'), through secular binary evolution, or through binary capture. In all three scenarios, the strong gravitational pertu…
▽ More
Dense stellar environments as hosts of ongoing star formation increase the probability of gravitational encounters among stellar systems during the early stages of evolution. Stellar interaction may occur through non-recurring, hyperbolic or parabolic passages (a so-called 'fly-by'), through secular binary evolution, or through binary capture. In all three scenarios, the strong gravitational perturbation is expected to manifest itself in the disc structures around the individual stars. Here, we present near-infrared polarised light observations that were taken with the SPHERE/IRDIS instrument of three known interacting twin-disc systems: AS 205, EM* SR 24, and FU Orionis. The scattered light exposes spirals likely caused by the gravitational interaction. On a larger scale, we observe connecting filaments between the stars. We analyse their very complex polarised intensity and put particular attention to the presence of multiple light sources in these systems. The local angle of linear polarisation indicates the source whose light dominates the scattering process from the bridging region between the two stars. Further, we show that the polarised intensity from scattering with multiple relevant light sources results from an incoherent summation of the individuals' contribution. This can produce nulls of polarised intensity in an image, as potentially observed in AS 205. We discuss the geometry and content of the systems by comparing the polarised light observations with other data at similar resolution, namely with ALMA continuum and gas emission. Collective observational data can constrain the systems' geometry and stellar trajectories, with the important potential to differentiate between dynamical scenarios of stellar interaction.
△ Less
Submitted 12 December, 2022; v1 submitted 25 November, 2022;
originally announced November 2022.
-
Rapid Formation of Massive Planetary Cores in a Pressure Bump
Authors:
Tommy Chi Ho Lau,
Joanna Drążkowska,
Sebastian M. Stammler,
Tilman Birnstiel,
Cornelis P. Dullemond
Abstract:
Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular s…
▽ More
Models of planetary core growth by either planetesimal or pebble accretion are traditionally disconnected from the models of dust evolution and formation of the first gravitationally-bound planetesimals. The state-of-the-art models typically start with massive planetary cores already present. We aim to study the formation and growth of planetary cores in a pressure bump, motivated by the annular structures observed in protoplanetary disks, starting with sub-micron-sized dust grains. We connect the models of dust coagulation and drift, planetesimal formation in the streaming instability, gravitational interactions between planetesimals, pebble accretion, and planet migration, into one uniform framework. We find that planetesimals forming early at the massive end of the size distribution grow quickly dominantly by pebble accretion. These few massive bodies grow on the timescales of ~100 000 years and stir the planetesimals formed later preventing the emergence of further planetary cores. Additionally, a migration trap occurs allowing for retention of the growing cores. Pressure bumps are favourable locations for the emergence and rapid growth of planetary cores by pebble accretion as the dust density and grain size are increased and the pebble accretion onset mass is reduced compared to a smooth-disk model.
△ Less
Submitted 21 December, 2022; v1 submitted 8 November, 2022;
originally announced November 2022.
-
Razor-thin dust layers in protoplanetary disks: Limits on the vertical shear instability
Authors:
C. P. Dullemond,
A. Ziampras,
D. Ostertag,
C. Dominik
Abstract:
Context: Recent observations with the Atacama Large Millimeter Array (ALMA) have shown that the large dust aggregates observed at millimeter wavelengths settle to the midplane into a remarkably thin layer. Aims: We intend to find out if the geometric thinness of these layers is evidence against the vertical shear instability (VSI) operating in these disks. Methods: We performed hydrodynamic simula…
▽ More
Context: Recent observations with the Atacama Large Millimeter Array (ALMA) have shown that the large dust aggregates observed at millimeter wavelengths settle to the midplane into a remarkably thin layer. Aims: We intend to find out if the geometric thinness of these layers is evidence against the vertical shear instability (VSI) operating in these disks. Methods: We performed hydrodynamic simulations of a protoplanetary disk with a locally isothermal equation of state, and let the VSI fully develop. We sprinkled dust particles and followed their motion as they got stirred up by the VSI. We determined for which grain size the layer becomes geometrically thin enough to be consistent with ALMA observations. We then verified if, with these grain sizes, it is still possible to generate a moderately optically thick layer at millimeter wavelengths, as observations appear to indicate. Results: We found that even very large dust aggregates with Stokes numbers close to unity get stirred up to relatively large heights above the midplane by the VSI, which is in conflict with the observed geometric thinness. For grains so large that the Stokes number exceeds unity, the layer can be made to remain thin, but we show that it is hard to make dust layers optically thick at ALMA wavelengths (e.g., tau(1.3mm)>=1) with such large dust aggregates. Conclusions: We conclude that protoplanetary disks with geometrically thin midplane dust layers cannot be VSI unstable, at least not down to the disk midplane. Explanations for the inhibition of the VSI include a reduced dust-to-gas ratio of the small dust grains that are responsible for the radiative cooling of the disk. A reduction of small grains by a factor of between 10 and 100 is sufficient to quench the VSI. Such a reduction is plausible in dust growth models, and still consistent with observations at optical and infrared wavelengths.
△ Less
Submitted 24 October, 2022;
originally announced October 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.
-
Dust entrainment in magnetically and thermally driven disk winds
Authors:
Peter J. Rodenkirch,
Cornelis P. Dullemond
Abstract:
Magnetically and thermally driven disk winds have gained popularity in the light of the current paradigm of low viscosities in protoplanetary disks that nevertheless present large accretion rates even in the presence of inner cavities. The possibility of dust entrainment in these winds may explain recent scattered light observations and constitutes a way of dust transport towards outer regions of…
▽ More
Magnetically and thermally driven disk winds have gained popularity in the light of the current paradigm of low viscosities in protoplanetary disks that nevertheless present large accretion rates even in the presence of inner cavities. The possibility of dust entrainment in these winds may explain recent scattered light observations and constitutes a way of dust transport towards outer regions of the disk. We aim to study the dust dynamics in these winds and explore the differences between photoevaporation and magnetically driven disk winds in this regard. We quantify maximum entrainable grain sizes, the flow angle, and the general detectability. We used the FARGO3D code to perform global, 2.5D axisymmetric, nonideal MHD simulations including ohmic and ambipolar diffusion. Dust was treated as a pressureless fluid. Synthetic observations were created with the radiative transfer code RADMC-3D. We find a significant difference in the dust entrainment efficiency of warm, ionized winds such as photoevaporation and magnetic winds including X-ray and extreme ultraviolet (XEUV) heating compared to cold magnetic winds. The maximum entrainable grain size varies from $3\,μ\mathrm{m}$ to $6\,μ\mathrm{m}$ for ionized winds to $1\,μ\mathrm{m}$ for cold magnetic winds. Dust grains in cold magnetic winds tend to flow along a shallower angle compared to the warm winds. With increasing distance to the central star, the dust entrainment efficiency decreases. Larger values of the turbulent viscosity increase the maximum grain size radius of possible dust entrainment. Our simulations indicate that diminishing dust content in the outer regions of the wind can be mainly attributed to the dust settling in the disk. In the synthetic images, the dusty wind appears as a faint, conical emission region which is brighter for a cold magnetic wind.
△ Less
Submitted 5 January, 2022;
originally announced January 2022.
-
Misaligned disks induced by infall
Authors:
M. Kuffmeier,
C. P. Dullemond,
S. Reissl,
F. G. Goicovic
Abstract:
Arc- and tail-like structures associated with disks around Herbig stars can be a consequence of infall events occurring after the initial collapse phase of a forming star. An encounter event of gas with a star can lead to the formation of a second-generation disk after the initial protostellar collapse phase. Additionally, observations of shadows in disks can be well described by a configuration o…
▽ More
Arc- and tail-like structures associated with disks around Herbig stars can be a consequence of infall events occurring after the initial collapse phase of a forming star. An encounter event of gas with a star can lead to the formation of a second-generation disk after the initial protostellar collapse phase. Additionally, observations of shadows in disks can be well described by a configuration of misaligned inner and outer disk, such that the inner disk casts a shadow on the outer disk. Carrying out altogether eleven 3D hydrodynamical models with the moving mesh code AREPO, we test whether a late encounter of an existing star-disk system with a cloudlet of gas can lead to the formation of an outer disk that is misaligned with respect to the primordial inner disk. Our models demonstrate that a second-generation disk with large misalignment with respect to an existing primordial disk can easily form if the infall angle is large. The second-generation outer disk is more eccentric, though the asymmetric infall also triggers eccentricity of the inner disk of $e\approx 0.05$ to $0.1$. Retrograde infall can lead to the formation of counter-rotating disks and enhanced accretion. As the angular momentum of the inner disk is reduced, the inner disk shrinks and a gap forms between the two disks. The resulting misaligned disk system can survive for $\sim 100$ kyr or longer without aligning each other even for low primordial disk masses given an infall mass of $\sim 10^{-4}$ M$_{\odot}$. A synthetic image reveals shadows in the outer disk similar to the ones observed in multiple transition disks that are caused by the misaligned inner disk. We conclude that late inclined infall onto a star-disk system leads to the formation of a misaligned outer disk. Infall might therefore be responsible for observations of shadows in at least some transition disks.
△ Less
Submitted 8 October, 2021;
originally announced October 2021.
-
On the equations of warped disc dynamics
Authors:
C. P. Dullemond,
C. N. Kimmig,
J. J. Zanazzi
Abstract:
The 1-D evolution equations for warped discs come in two flavors: For very viscous discs the internal torque vector G is uniquely determined by the local conditions in the disc, and warps tend to damp out rapidly if they are not continuously driven. For very inviscid discs, on the other hand, G becomes a dynamic quantity, and a warp will propagate through the disc as a wave. The equations governin…
▽ More
The 1-D evolution equations for warped discs come in two flavors: For very viscous discs the internal torque vector G is uniquely determined by the local conditions in the disc, and warps tend to damp out rapidly if they are not continuously driven. For very inviscid discs, on the other hand, G becomes a dynamic quantity, and a warp will propagate through the disc as a wave. The equations governing both regimes are usually treated separately. A unified set of equations was postulated recently by Martin et al. (2019), but not yet derived from the underlying physics. The standard method for deriving these equations is based on a perturbation series expansion, which is a powerful, but somewhat abstract technique. A more straightforward method is to employ the warped shearing box framework of Ogilvie and Latter (2013), which so far has not yet been used to derive the equations for the wavelike regime. The goal of this paper is to analyze the warped disc equations in both regimes using the warped shearing box framework, to derive a unified set of equations, valid for small warps, and to discuss how our results can be interpreted in terms of the affine tilted-slab approach of Ogilvie (2018).
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
Self-Sustaining Vortices in Protoplanetary Disks: Setting the Stage for Planetary System Formation
Authors:
Zsolt Regaly,
Kundan Kadam,
Cornelis P. Dullemond
Abstract:
The core accretion scenario of planet formation assumes that planetesimals and planetary embryos are formed during the primordial, gaseous phases of the protoplanetary disk. However, how the dust particles overcome the traditional growth barriers is not well understood. The recently proposed viscous ring-instability may explain the concentric rings observed in protoplanetary disks by assuming that…
▽ More
The core accretion scenario of planet formation assumes that planetesimals and planetary embryos are formed during the primordial, gaseous phases of the protoplanetary disk. However, how the dust particles overcome the traditional growth barriers is not well understood. The recently proposed viscous ring-instability may explain the concentric rings observed in protoplanetary disks by assuming that the dust grains can reduce the gas conductivity, which can weaken the magneto-rotational instability. We present an analysis of this model with the help of GPU-based numerical hydrodynamic simulations of coupled gas and dust in the thin-disk limit. During the evolution of the disk the dusty rings become Rossby unstable and break up into a cascade of small-scale vortices. The vortices form secularly stable dusty structures, which could be sites of planetesimal formation by the streaming instability as well as direct gravitational collapse. The phenomenon of self-sustaining vortices is consistent with observational constraints of exoplanets and sets a favorable environment for planetary system formation.
△ Less
Submitted 19 July, 2021; v1 submitted 26 June, 2021;
originally announced June 2021.
-
Modeling the nonaxisymmetric structure in the HD 163296 disk with planet-disk interaction
Authors:
Peter J. Rodenkirch,
Thomas Rometsch,
Cornelis P. Dullemond,
Philipp Weber,
Wilhelm Kley
Abstract:
High resolution ALMA observations revealed a variety of rich substructures in numerous protoplanetary disks. These structures consist of rings, gaps and asymmetric features. It is debated whether planets can be accounted for these substructures in the dust continuum. Characterizing the origin of asymmetries as seen in HD 163296 might lead to a better understanding of planet formation and the under…
▽ More
High resolution ALMA observations revealed a variety of rich substructures in numerous protoplanetary disks. These structures consist of rings, gaps and asymmetric features. It is debated whether planets can be accounted for these substructures in the dust continuum. Characterizing the origin of asymmetries as seen in HD 163296 might lead to a better understanding of planet formation and the underlying physical parameters of the system. We test the possibility of the formation of the crescent-shaped asymmetry in the HD 163296 disk through planet-disk interaction. The goal is to obtain constraints on planet masses and eccentricities and disk viscosities. Two dimensional, multi-fluid, hydrodynamical simulations are performed with the FARGO3D code including three embedded planets. Dust is described with the pressureless fluid approach and is distributed over eight size bins. Resulting grids are post-processed with the radiative transfer code RADMC-3D and the CASA software to model synthetic observations. We find that the crescent-shaped asymmetry can be qualitatively modeled with a Jupiter mass planet at a radial distance of 48 au. Dust is trapped preferably in the trailing Lagrange point L5 with a mass of 10 to 15 earth masses. Increased values of eccentricity of the innermost Jupiter mass planet damages the stability of the crescent-shaped feature and does not reproduce the observed radial proximity to the first prominent ring in the system. Generally, a low level of viscosity ($α\leq 2\cdot10^{-3}$) is necessary to allow the existence of such a feature. Including dust feedback the leading point L4 can dominantly capture dust for dust grains with an initial Stokes number $\leq 3.6\cdot 10^{-2}$. The observational results suggest a negligible effect of dust feedback since only one such feature has been detected so far.
△ Less
Submitted 16 December, 2020;
originally announced December 2020.
-
Migration jumps of planets in transition disks
Authors:
Thomas Rometsch,
Peter J. Rodenkirch,
Wilhelm Kley,
Cornelis P. Dullemond
Abstract:
Transition disks form a special class of protoplanetary disks that are characterized by a deficiency of disk material close to the star. In a subgroup, inner holes in these disks can stretch out to a few tens of au while there is still mass accretion onto the central star observed. We analyse the proposition that this type of wide transition disks is generated by the interaction of the disk with a…
▽ More
Transition disks form a special class of protoplanetary disks that are characterized by a deficiency of disk material close to the star. In a subgroup, inner holes in these disks can stretch out to a few tens of au while there is still mass accretion onto the central star observed. We analyse the proposition that this type of wide transition disks is generated by the interaction of the disk with a system of embedded planets. We performed 2D hydrodynamics simulations of a flat disk using either a locally isothermal equation of state or considering also radiative effects. Two 3 to 9 Jupiter mass planets were embedded in the disk and their dynamical evolution due to disk-planet interaction was followed for over 100 000 years. The simulations account for mass accretion onto the star and planets. We included models with parameters geared to the system PDS 70. To assess the observability of features in our models we performed synthetic ALMA observations.
For systems with a more massive inner planet there are phases where both planets migrate outward engaged in a 2:1 mean motion resonance via the Masset-Snellgrove mechanism. In sufficiently massive disks the formation of a vortex in the outer disk can trigger rapid outward migration of the outer planet where its distance increases by tens of au within a few thousand years. Later, the outer planet migrates back inwards settling again into resonance with the inner planet. We call this emerging composite phenomenon a 'migration jump'. Outward migration and the migration jumps are accompanied by a high mass accretion rate onto the star. The synthetic images reveal numerous substructures depending on the type of dynamical behaviour. Our results suggest that the outward migration of two embedded planets is a prime candidate for the explanation of the observed high stellar mass accretion rate in wide transition disks.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
The efficiency of dust trapping in ringed proto-planetary discs
Authors:
Giovanni P. Rosotti,
Richard Teague,
Cornelis Dullemond,
Richard A. Booth,
Cathie Clarke
Abstract:
When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $α/St$ and the width of the gas press…
▽ More
When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $α/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($α/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $α$ turbulent parameter ($α\sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.
△ Less
Submitted 23 April, 2020;
originally announced April 2020.
-
Importance of radiative effects in gap opening by planets in protoplanetary disks
Authors:
Alexandros Ziampras,
Wilhelm Kley,
Cornelis P. Dullemond
Abstract:
Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100…
▽ More
Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100 au). We test this assumption by conducting hydrodynamics simulations of embedded planets in thin locally isothermal and radiative disks that mimic the systems HD 163296 and AS 209 in order to examine the effect of including the energy equation in a seemingly locally isothermal environment as far as planet-disk interaction is concerned. We find that modeling such disks with an ideal equation of state makes a difference in terms of the number of produced rings and the spiral arm contrast in the disk. Locally isothermal disks produce sharper annular or azimuthal features and overestimate a single planet's gap-opening capabilities by producing multiple gaps. In contrast, planets in radiative disks carve a single gap for typical disk parameters. Consequently, for accurate modeling of planets with semimajor axes up to about 100 au, radiative effects should be taken into account even in seemingly locally isothermal disks. In addition, for the case of AS 209, we find that the primary gap is significantly different between locally isothermal and radiative models. Our results suggest that multiple planets are required to explain the ring-rich structures in such systems.
△ Less
Submitted 4 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.
-
Midplane temperature and outer edge of the protoplanetary disk around HD 163296
Authors:
Cornelis Dullemond,
Andrea Isella,
Sean Andrews,
Iuliia Skobleva,
Natalia Dzyurkevich
Abstract:
Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the…
▽ More
Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, also in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2-1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be an indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.
△ Less
Submitted 27 November, 2019;
originally announced November 2019.
-
Cloudlet capture by Transitional Disk and FU Orionis stars
Authors:
Cornelis Petrus Dullemond,
Michael Küffmeier,
Felipe Goicovic,
Misato Fukagawa,
Veronika Oehl,
Manuel Kramer
Abstract:
After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete ont…
▽ More
After its formation, a young star spends some time traversing the molecular cloud complex in which it was born. It is therefore not unlikely that, well after the initial cloud collapse event which produced the star, it will encounter one or more low mass cloud fragments, which we call "cloudlets" to distinguish them from full-fledged molecular clouds. Some of this cloudlet material may accrete onto the star+disk system, while other material may fly by in a hyperbolic orbit. In contrast to the original cloud collapse event, this process will be a "cloudlet flyby" and/or "cloudlet capture" event: A Bondi-Hoyle-Lyttleton type accretion event, driven by the relative velocity between the star and the cloudlet. As we will show in this paper, if the cloudlet is small enough and has an impact parameter similar or less than $GM_{*}/v^2_\infty$ (with $v_\infty$ being the approach velocity), such a flyby and/or capture event would lead to arc-shaped or tail-shaped reflection nebulosity near the star. Those shapes of reflection nebulosity can be seen around several transitional disks and FU Orionis stars. Although the masses in the those arcs appears to be much less than the disk masses in these sources, we speculate that higher-mass cloudlet capture events may also happen occasionally. If so, they may lead to the tilting of the outer disk, because the newly infalling matter will have an angular momentum orientation entirely unrelated to that of the disk. This may be one possible explanation for the highly warped/tilted inner/outer disk geometries found in several transitional disks. We also speculate that such events, if massive enough, may lead to FU Orionis outbursts.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Late encounter-events as a source of disks and spiral structures -- Forming second generation disks
Authors:
M. Kuffmeier,
F. G. Goicovic,
C. P. Dullemond
Abstract:
Both observations of arc-like structures and luminosity bursts of stars > 1 Myr in age indicate that at least some stars undergo late infall events. We investigate scenarios of replenishing the mass reservoir around a star via capturing and infalling events of cloudlets. We carry out altogether 24 three-dimensional hydrodynamical simulations of cloudlet encounters with a Herbig star of mass 2.5 so…
▽ More
Both observations of arc-like structures and luminosity bursts of stars > 1 Myr in age indicate that at least some stars undergo late infall events. We investigate scenarios of replenishing the mass reservoir around a star via capturing and infalling events of cloudlets. We carry out altogether 24 three-dimensional hydrodynamical simulations of cloudlet encounters with a Herbig star of mass 2.5 solar mass using the moving mesh code AREPO. To account for the two possibilities of a star or a cloudlet traveling through the interstellar medium (ISM), we put either the star or the cloudlet at rest with respect to the background gas. For absent cooling in the adiabatic runs, almost none of the cloudlet gas is captured due to high thermal pressure. However, second-generation disks easily form when accounting for cooling of the gas. The disk radii range from several 100 au to about 1000 au and associated arc-like structures up to 10 000 au in length form around the star for runs with and without stellar irradiation. Consistent with angular momentum conservation, the arcs and disks are larger for larger impact parameters. Accounting for turbulence in the cloudlet only mildly changes the model outcome. In the case of the star being at rest with the background gas, the disk formation and mass replenishment process is more pronounced and the associated arc-shaped streamers are longer-lived. The results of our models confirm that late encounter events lead to the formation of transitional disks associated with arc-shaped structures such as observed for AB Aurigae or HD 100546. In addition, we find that second-generation disks and their associated filamentary arms are longer lived (>100 000 yrs) in infall events, when the star is at rest with the background gas.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Global axisymmetric simulations of photoevaporation and magnetically driven protoplanetary disk winds
Authors:
Peter J. Rodenkirch,
Hubert Klahr,
Christian Fendt,
Cornelis P. Dullemond
Abstract:
Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present…
▽ More
Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $β\geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $β< 10^7$.
△ Less
Submitted 15 January, 2020; v1 submitted 11 November, 2019;
originally announced November 2019.
-
Spiral arms in the proto-planetary disc HD100453 detected with ALMA: evidence for binary-disc interaction and a vertical temperature gradient
Authors:
G. P. Rosotti,
M. Benisty,
A. Juhász,
R. Teague,
C. Clarke,
C. Dominik,
C. P. Dullemond,
P. D. Klaassen,
L. Matrà,
T. Stolker
Abstract:
Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J…
▽ More
Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.
△ Less
Submitted 1 November, 2019;
originally announced November 2019.
-
Effect of wind-driven accretion on planetary migration
Authors:
C. N. Kimmig,
C. P. Dullemond,
W. Kley
Abstract:
Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a r…
▽ More
Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a renaissance, as evidence is mounting that PPDs may be less turbulent than previously thought, and 3-D non-ideal MHD modeling of the wind-launching process is maturing. Aim: We wish to investigate how wind-driven accretion may affect migration. We aim for a qualitative exploration of the main effects, rather than a quantitative prediction. Methods: We perform 2-D hydrodynamic simulations with the FARGO3D code in the $(r,φ)$-plane. The vertical coordinate and the launching of the wind are not treated explicitly. Instead, the torque of the wind onto the disk is treated using a simple 2-parameter formula treating the wind mass loss rate and the lever arm. Results: We find that the wind-driven accretion process has a different way of replenishing the co-orbital region than the viscous accretion. The former always injects mass from the outer edge of the co-orbital region and removes mass from the inner edge, while the latter injects or removes mass from the co-orbital region depending on the radial density gradients in the disk. The migration behavior can differ very much and under certain conditions it can drive rapid type-III-like outward migration. We derive an analytic expression for the parameters under which this outward migration occurs. Conclusion: If wind-driven accretion plays a role in the secular evolution of PPDs, migration studies have to include this process as well, because it can strongly affect the resulting migration rate and direction.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
The impact of planet wakes on the location and shape of the water iceline in a protoplanetary disk
Authors:
Alexandros Ziampras,
Sareh Ataiee,
Wilhelm Kley,
Cornelis P. Dullemond,
Clément Baruteau
Abstract:
Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure.
Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric natu…
▽ More
Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure.
Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets.
Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters.
Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context.
Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.
△ Less
Submitted 18 October, 2019;
originally announced October 2019.
-
An Ideal Testbed for Planet-disk Interaction: Two Giant Protoplanets in Resonance Shaping the PDS 70 Protoplanetary Disk
Authors:
Jaehan Bae,
Zhaohuan Zhu,
Clément Baruteau,
Myriam Benisty,
Cornelis P. Dullemond,
Stefano Facchini,
Andrea Isella,
Miriam Keppler,
Laura M. Pérez,
Richard Teague
Abstract:
While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features…
▽ More
While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$α$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$μ$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$μ$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$μ$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems.
△ Less
Submitted 20 September, 2019;
originally announced September 2019.
-
The DSHARP Rings: Evidence of Ongoing Planetesimal Formation?
Authors:
Sebastian M. Stammler,
Joanna Drazkowska,
Til Birnstiel,
Hubert Klahr,
Cornelis P. Dullemond,
Sean M. Andrews
Abstract:
Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of…
▽ More
Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process.
We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations.
We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Observability of Forming Planets and their Circumplanetary Disks II. -- SEDs and Near-Infrared Fluxes
Authors:
J. Szulágyi,
C. P. Dullemond,
A. Pohl,
S. P. Quanz
Abstract:
Detection of forming planets means detection of the circumplanetary disk (CPD) in reality, since the planet is still surrounded by a disk at this evolutionary stage. Yet, no comprehensive CPD modeling was done in near-infrared wavelengths, where high contrast imaging is a powerful tool to detect these objects. We combined 3D radiative hydrodynamic simulations of various embedded planets with RADMC…
▽ More
Detection of forming planets means detection of the circumplanetary disk (CPD) in reality, since the planet is still surrounded by a disk at this evolutionary stage. Yet, no comprehensive CPD modeling was done in near-infrared wavelengths, where high contrast imaging is a powerful tool to detect these objects. We combined 3D radiative hydrodynamic simulations of various embedded planets with RADMC-3D radiative transfer post-processing that includes scattering of photons on dust particles. We made synthetic images for VLT NaCo/ERIS in the Ks, L', M' bands as well as examined the spectral energy distributions (SEDs) of disks between 1 $μm$ and 10 cm. We found that the observed magnitudes from the planet's vicinity will mostly depend on the CPD parameters, not on the planet's. The CPD is 20-100x brighter than the embedded planet in near-IR. We also show how the CPD parameters, e.g. the dust-to-gas ratio will affect the resulting CPD magnitudes. According to the SEDs, the best contrast ratio between the CPD and circumstellar disks is in sub-mm/radio wavelengths and between 8-33 microns in case if the planet opened a resolvable, deep gap ($\ge 5 \rm{M_{Jup}}$), while the contrast is particularly poor in the near-IR. Hence, to detect the forming planet and its CPD, the best chance today is targeting the sub-mm/radio wavelengths and the 10-micron silicate feature vicinity. In order to estimate the forming planet's mass from the observed brightness, it is necessary to run system specific disk modeling.
△ Less
Submitted 9 May, 2019;
originally announced May 2019.
-
One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering
Authors:
Zhaohuan Zhu,
Shangjia Zhang,
Yan-Fei Jiang,
Akimasa Kataoka,
Tilman Birnstiel,
Cornelis P. Dullemond,
Sean M. Andrews,
Jane Huang,
Laura M. Perez,
John M. Carpenter,
Xue-Ning Bai,
David J. Wilner,
Luca Ricci
Abstract:
ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out tha…
▽ More
ALMA surveys have suggested that the dust in Class II disks may not be enough to explain the averaged solid mass in exoplanets, under the assumption that the mm disk continuum emission is optically thin. This optically thin assumption seems to be supported by recent DSHARP observations where the measured optical depths of spatially resolved disks are mostly less than one. However, we point out that dust scattering can considerably reduce the emission from an optically thick region. If that scattering is ignored, the optical depth will be considerably underestimated. An optically thick disk with scattering can be misidentified as an optically thin disk. Dust scattering in more inclined disks can reduce the intensity even further, making the disk look even fainter. The measured optical depth of $\sim$0.6 in several DSHARP disks can be naturally explained by optically thick dust with an albedo of $\sim$0.9 at 1.25 mm. Using the DSHARP opacity, this albedo corresponds to a dust population with the maximum grain size ($s_{max}$) of 0.1-1 mm. For optically thick scattering disks, the measured spectral index $α$ can be either larger or smaller than 2 depending on if the dust albedo increases or decreases with wavelength. Using the DSHARP opacity, $α<2$ corresponds to $s_{max}$ of 0.03-0.3 mm. We describe how this optically thick scattering scenario could explain the observed scaling between submm continuum sizes and luminosities, and might help ease the tension between the dust size constraints from polarization and dust continuum measurements. We suggest that a significant amount of disk mass can be hidden from ALMA observations at short millimeter wavelengths. For compact disks smaller than 30 au, we can easily underestimate the dust mass by more than a factor of 10. Longer wavelength observations (e.g. VLA or SKA) are desired to probe the dust mass in disks.
△ Less
Submitted 7 May, 2019; v1 submitted 3 April, 2019;
originally announced April 2019.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): IV. Characterizing substructures and interactions in disks around multiple star systems
Authors:
Nicolás Kurtovic,
Laura Pérez,
Myriam Benisty,
Zhaohuan Zhu,
Shangjia Zhang,
Jane Huang,
Sean M. Andrews,
Cornellis P. Dullemond,
Andrea Isella,
Xue-Ning Bai,
John M. Carpenter,
Viviana V. Guzmán,
Luca Ricci,
David J. Wilner
Abstract:
To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symme…
▽ More
To characterize the substructures induced in protoplanetary disks by the interaction between stars in multiple systems, we study the $1.25\,$mm continuum and the $^{12}$CO$(J=2-1)$ spectral line emission of the triple systems HT Lup and AS 205, at scales of $\approx 5\,$au, as part of the 'Disk Substructures at High Angular Resolution Project' (DSHARP). In the continuum emission, we find two symmetric spiral arms in the disk around AS 205 N, with pitch angle of $14^\circ$, while the southern component AS 205 S, itself a spectroscopic binary, is surrounded by a compact inner disk and a bright ring at a radius of $34\,$au. The $^{12}$CO line exhibits clear signatures of tidal interactions, with spiral arms, extended arc-like emission, and high velocity gas, possible evidence of a recent close encounter between the disks in the AS 205 system, as these features are predicted by hydrodynamic simulations of fly-by encounters. In the HT Lup system, we detect continuum emission from all three components. The primary disk, HT Lup A, also shows two-armed symmetric spiral structure with a pitch angle of $4^\circ$, while HT Lup B and C, located at $25$ and $434\,$au in projected separation from HT Lup A, are barely resolved with $\sim5$ and $\sim10\,$au in diameter, respectively. The gas kinematics for the closest pair indicates a different sense of rotation for each disk, which could be explained by either a counter rotation of the two disks in different, close to parallel, planes, or by a projection effect of these disks with a close to $90^\circ$ misalignment between them.
△ Less
Submitted 11 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): III. Spiral Structures in the Millimeter Continuum of the Elias 27, IM Lup, and WaOph 6 Disks
Authors:
Jane Huang,
Sean M. Andrews,
Laura M. Pérez,
Zhaohuan Zhu,
Cornelis P. Dullemond,
Andrea Isella,
Myriam Benisty,
Xue-Ning Bai,
Tilman Birnstiel,
John M. Carpenter,
Viviana V. Guzmán,
A. Meredith Hughes,
Karin I. Öberg,
Luca Ricci,
David J. Wilner,
Shangjia Zhang
Abstract:
We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics…
▽ More
We present an analysis of ALMA 1.25 millimeter continuum observations of spiral structures in three protoplanetary disks from the Disk Substructures at High Angular Resolution Project. The disks around Elias 27, IM Lup, and WaOph 6 were observed at a resolution of $\sim40-60$ mas ($\sim6-7$ au). All three disks feature $m=2$ spiral patterns in conjunction with annular substructures. Gas kinematics established by $^{12}$CO $J=2-1$ observations indicate that the continuum spiral arms are trailing. The arm-interarm intensity contrasts are modest, typically less than 3. The Elias 27 spiral pattern extends throughout much of the disk, and the arms intersect the gap at $R\sim69$ au. The spiral pattern in the IM Lup disk is particularly complex-it extends about halfway radially through the disk, exhibiting pitch angle variations with radius and interarm features that may be part of ring substructures or spiral arm branches. Spiral arms also extend most of the way through the WaOph 6 disk, but the source overall is much more compact than the other two disks. We discuss possible origins for the spiral structures, including gravitational instability and density waves induced by a stellar or planetary companion. Unlike the millimeter continuum counterparts of many of the disks with spiral arms detected in scattered light, these three sources do not feature high-contrast crescent-like asymmetries or large ($R>20$ au) emission cavities. This difference may point to multiple spiral formation mechanisms operating in disks.
△ Less
Submitted 14 January, 2019; v1 submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): X. Multiple rings, a misaligned inner disk, and a bright arc in the disk around the T Tauri star HD 143006
Authors:
Laura M. Pérez,
Myriam Benisty,
Sean M. Andrews,
Andrea Isella,
Cornelis P. Dullemond,
Jane Huang,
Nicolás T. Kurtovic,
Viviana V. Guzmán,
Zhaohuan Zhu,
Tilman Birnstiel,
Shangjia Zhang,
John M. Carpenter,
David J. Wilner,
Luca Ricci,
Xue-Ning Bai,
Erik Weaver,
Karin I. Öberg
Abstract:
We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting…
▽ More
We present a detailed analysis of new ALMA observations of the disk around the T-Tauri star HD 143006, which at 46 mas (7.6 au) resolution reveal new substructures in the 1.25 mm continuum emission. The disk resolves into a series of concentric rings and gaps together with a bright arc exterior to the rings that resembles hydrodynamics simulations of a vortex, and a bridge-like feature connecting the two innermost rings. Although our $^{12}$CO observations at similar spatial resolution do not show obvious substructure, they reveal an inner disk depleted of CO emission. From the continuum emission and the CO velocity field we find that the innermost ring has a higher inclination than the outermost rings and the arc. This is evidence for either a small ($\sim8^{\circ}$) or moderate ($\sim41^{\circ}$) misalignment between the inner and outer disk, depending on the specific orientation of the near/far sides of the inner/outer disk. We compare the observed substructures in the ALMA observations with recent scattered light data from VLT/SPHERE of this object. In particular, the location of narrow shadow lanes in the SPHERE image combined with pressure scale height estimates, favor a large misalignment of about $41^{\circ}$. We discuss our findings in the context of a dust-trapping vortex, planet-carved gaps, and a misaligned inner disk due to the presence of an inclined companion to HD 143006.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP) - IX. A high definition study of the HD 163296 planet forming disk
Authors:
Andrea Isella,
Jane Huang,
Sean M. Andrews,
Cornelis P. Dullemond,
Tilman Birnstiel,
Shangjia Zhang,
Zhaohuan Zhu,
Viviana V. Guzmán,
Laura M. Pérez,
Xue-Ning Bai,
Myriam Benisty,
John M. Carpenter,
Luca Ricci,
David J. Wilner
Abstract:
ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm d…
▽ More
ALMA observations of protoplanetary disks acquired by the Disk Substructure at High Angular Resolution Project (DSHARP) resolve the dust and gas emission on angular scales as small as 3 astronomical units, offering an unprecedented detailed view of the environment where planets form. In this article, we present and discuss observations of the HD 163296 protoplanetary disk that imaged the 1.25 mm dust continuum and $^{12}$CO J=2-1 rotational line emission at a spatial resolution of 4 and 10 au, respectively. The continuum observations resolve and allow us to characterize the previously discovered dust rings at radii of 67 and 100 au. They also reveal new small scale structures, such as a dark gap at 10 au, a bright ring at 15 au, a dust crescent at a radius of 55 au, and several fainter azimuthal asymmetries. The observations of the CO and dust emission inform about the vertical structure of the disk and allow us to directly constrain the dust extinction optical depth at the dust rings. Furthermore, the observed asymmetries in the dust continuum emission corroborate to the hypothesis that the complex structure of the HD 163296 disk is the result of the gravitational interaction with yet unseen planets.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
Disk Substructures at High Angular Resolution Program (DSHARP): VIII. The Rich Ringed Substructures in the AS 209 Disk
Authors:
Viviana V. Guzmán,
Jane Huang,
Sean M. Andrews,
Andrea Isella,
Laura M. Pérez,
John M. Carpenter,
Cornelis P. Dullemond,
Luca Ricci,
Tilman Birnstiel,
Shangjia Zhang,
Zhaohuan Zhu,
Xue-Ning Bai,
Myriam Benisty,
Karin I. Öberg,
David J. Wilner
Abstract:
We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained…
▽ More
We present a detailed analysis of the high-angular resolution (0.''037, corresponding to 5 au) observations of the 1.25 mm continuum and $^{12}$CO $2-1$ emission from the disk around the T Tauri star AS 209. AS 209 hosts one of the most unusual disks from the DSHARP sample, the first high angular resolution ALMA survey of disks (Andrews et al. 2018), as nearly all of the emission can be explained with concentric Gaussian rings. In particular, the dust emission consists of a series of narrow and closely spaced rings in the inner $\sim$60 au, two well-separated bright rings in the outer disk, centered at 74 and 120 au, and at least two fainter emission features at 90 and 130 au. We model the visibilities with a parametric representation of the radial surface brightness profile, consisting of a central core and 7 concentric Gaussian rings. Recent hydro-dynamical simulations of low viscosity disks show that super-Earth planets can produce the multiple gaps seen in AS 209 millimeter continuum emission. The $^{12}$CO line emission is centrally peaked and extends out to $\sim$300 au, much farther than the millimeter dust emission. We find axisymmetric, localized deficits of CO emission around four distinct radii, near 45, 75, 120 and 210 au. The outermost gap is located well beyond the edge of the millimeter dust emission, and therefore cannot be due to dust opacity and must be caused by a genuine CO surface density reduction, due either to chemical effects or depletion of the overall gas content.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): VII. The Planet-Disk Interactions Interpretation
Authors:
Shangjia Zhang,
Zhaohuan Zhu,
Jane Huang,
Viviana V. Guzmán,
Sean M. Andrews,
Tilman Birnstiel,
Cornelis P. Dullemond,
John M. Carpenter,
Andrea Isella,
Laura M. Pérez,
Myriam Benisty,
David J. Wilner,
Clément Baruteau,
Xue-Ning Bai,
Luca Ricci
Abstract:
The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resul…
▽ More
The Disk Substructures at High Angular Resolution Project (DSHARP) provides a large sample of protoplanetary disks having substructures which could be induced by young forming planets. To explore the properties of planets that may be responsible for these substructures, we systematically carry out a grid of 2-D hydrodynamical simulations including both gas and dust components. We present the resulting gas structures, including the relationship between the planet mass and 1) the gaseous gap depth/width, and 2) the sub/super-Keplerian motion across the gap. We then compute dust continuum intensity maps at the frequency of the DSHARP observations. We provide the relationship between the planet mass and 1) the depth/width of the gaps at millimeter intensity maps, 2) the gap edge ellipticity and asymmetry, and 3) the position of secondary gaps induced by the planet. With these relationships, we lay out the procedure to constrain the planet mass using gap properties, and study the potential planets in the DSHARP disks. We highlight the excellent agreement between observations and simulations for AS 209 and the detectability of the young Solar System analog. Finally, under the assumption that the detected gaps are induced by young planets, we characterize the young planet population in the planet mass-semimajor axis diagram. We find that the occurrence rate for $>$ 5 $M_J$ planets beyond 5-10 au is consistent with direct imaging constraints. Disk substructures allow us probe a wide-orbit planet population (Neptune to Jupiter mass planets beyond 10 au) that is not accessible to other planet searching techniques.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP) VI: Dust trapping in thin-ringed protoplanetary disks
Authors:
Cornelis P. Dullemond,
Tilman Birnstiel,
Jane Huang,
Nicolás T. Kurtovic,
Sean M. Andrews,
Viviana V. Guzmán,
Laura M. Pérez,
Andrea Isella,
Zhaohuan Zhu,
Myriam Benisty,
David J. Wilner,
Xue-Ning Bai,
John M. Carpenter,
Shangjia Zhang,
Luca Ricci
Abstract:
A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bu…
▽ More
A large fraction of the protoplanetary disks observed with ALMA display multiple well-defined and nearly perfectly circular rings in the continuum, in many cases with substantial peak-to-valley contrast. The DSHARP campaign shows that several of these rings are very narrow in radial extent. In this paper we test the hypothesis that these dust rings are caused by dust trapping in radial pressure bumps, and if confirmed, put constraints on the physics of the dust trapping mechanism. We model this process analytically in 1D, assuming axisymmetry. By comparing this model to the data, we find that all rings are consistent with dust trapping. Based on a plausible model of the dust temperature we find that several rings are narrower than the pressure scale height, providing strong evidence for dust trapping. The rings have peak absorption optical depth in the range between 0.2 and 0.5. The dust masses stored in each of these rings is of the order of tens of Earth masses, though much ambiguity remains due to the uncertainty of the dust opacities. The dust rings are dense enough to potentially trigger the streaming instability, but our analysis cannot give proof of this mechanism actually operating. Our results show, however, that the combination of very low alpha_turb << 5e-4 and very large grains a_grain >> 0.1 cm can be excluded by the data for all the rings studied in this paper.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): V. Interpreting ALMA maps of protoplanetary disks in terms of a dust model
Authors:
Tilman Birnstiel,
Cornelis P. Dullemond,
Zhaohuan Zhu,
Sean M. Andrews,
Xue-Ning Bai,
David J. Wilner,
John M. Carpenter,
Jane Huang,
Andrea Isella,
Myriam Benisty,
Laura M. Pérez,
Shangjia Zhang
Abstract:
The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very n…
▽ More
The Disk Substructures at High Angular Resolution Project (DSHARP) is the largest homogeneous high-resolution ($\sim 0.035$ arcsec, or $\sim$ 5 au) disk continuum imaging survey with ALMA so far. In the coming years, many more disks will be mapped with ALMA at similar resolution. Interpreting the results in terms of the properties and quantities of the emitting dusty material is, however, a very non-trivial task. This is in part due to the uncertainty in the dust opacities, an uncertainty which is not likely to be resolved any time soon. It is also partly due to the fact that, as the DSHARP survey has shown, these disk often contain regions of intermediate to high optical depth, even at millimeter wavelengths and at relatively large radius in the disk. This makes the interpretation challenging, in particular if the grains are large and have a large albedo. On the other hand, the highly structured features seen in the DSHARP survey, of which strong indications were already seen in earlier observations, provide a unique opportunity to study the dust growth and dynamics. To provide continuity within the DSHARP project, its follow-up projects, and projects by other teams interested in these data, we present here the methods and opacity choices used within the DSHARP collaboration to link the measured intensity $I_ν$ to dust surface density $Σ_d$.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures
Authors:
Jane Huang,
Sean M. Andrews,
Cornelis P. Dullemond,
Andrea Isella,
Laura M. Pérez,
Viviana V. Guzmán,
Karin I. Öberg,
Zhaohuan Zhu,
Shangjia Zhang,
Xue-Ning Bai,
Myriam Benisty,
Tilman Birnstiel,
John M. Carpenter,
A. Meredith Hughes,
Luca Ricci,
Erik Weaver,
David J. Wilner
Abstract:
The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructu…
▽ More
The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.
△ Less
Submitted 12 August, 2020; v1 submitted 10 December, 2018;
originally announced December 2018.
-
The Disk Substructures at High Angular Resolution Project (DSHARP): I. Motivation, Sample, Calibration, and Overview
Authors:
Sean M. Andrews,
Jane Huang,
Laura M. Pérez,
Andrea Isella,
Cornelis P. Dullemond,
Nicolás T. Kurtovic,
Viviana V. Guzmán,
John M. Carpenter,
David J. Wilner,
Shangjia Zhang,
Zhaohuan Zhu,
Tilman Birnstiel,
Xue-Ning Bai,
Myriam Benisty,
A. Meredith Hughes,
Karin I. Öberg,
Luca Ricci
Abstract:
We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 a…
▽ More
We introduce the Disk Substructures at High Angular Resolution Project (DSHARP), one of the initial Large Programs conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal of DSHARP is to find and characterize substructures in the spatial distributions of solid particles for a sample of 20 nearby protoplanetary disks, using very high resolution (0.035 arcsec, or 5 au FWHM) observations of their 240 GHz (1.25 mm) continuum emission. These data provide a first homogeneous look at the small-scale features in disks that are directly relevant to the planet formation process, quantifying their prevalence, morphologies, spatial scales, spacings, symmetry, and amplitudes, for targets with a variety of disk and stellar host properties. We find that these substructures are ubiquitous in this sample of large, bright disks. They are most frequently manifested as concentric, narrow emission rings and depleted gaps, although large-scale spiral patterns and small arc-shaped azimuthal asymmetries are also present in some cases. These substructures are found at a wide range of disk radii (from a few au to more than 100 au), are usually compact ($<$10 au), and show a wide range of amplitudes (brightness contrasts). Here we discuss the motivation for the project, describe the survey design and the sample properties, detail the observations and data calibration, highlight some basic results, and provide a general overview of the key conclusions that are presented in more detail in a series of accompanying articles. The DSHARP data -- including visibilities, images, calibration scripts, and more -- are released for community use at https://almascience.org/alma-data/lp/DSHARP.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
Surface waves in protoplanetary disks induced by outbursts: Concentric rings in scattered light
Authors:
A. D. Schneider,
C. P. Dullemond,
B. Bitsch
Abstract:
Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical…
▽ More
Context: Vertically hydrostatic protoplanetary disk models are based on the assumption that the main heating source, stellar irradiation, does not vary much with time. However, it is known that accreting young stars are variable sources of radiation. This is particularly evident for outbursting sources such as EX Lupi and FU Orionis stars. Aim: We investigate how such outbursts affect the vertical structure of the outer regions of the protoplanetary disk, in particular their appearance in scattered light at optical and near-infrared wavelengths. Methods: We employ the 3D FARGOCA radiation-hydrodynamics code, in polar coordinates, to compute the time-dependent behavior of the axisymmetric disk structure. The outbursting inner disk region is not included explicitly. Instead, its luminosity is added to the stellar luminosity and is thus included in the irradiation of the outer disk regions. For time snapshots of interest we insert the density structure into the RADMC-3D radiative transfer code and compute the appearance of the disk at optical/near-infrared wavelengths. Results: We find that, depending on the amplitude of the outbursts, the vertical structure of the disk can become highly dynamic, featuring circular surface waves of considerable amplitude. These "hills" and "valleys" on the disk's surface show up in the scattered light images as bright and dark concentric rings. Initially these rings are small and act as standing waves, but they subsequently lead to outward propagating waves, like the waves produced by a stone thrown into a pond. These waves continue long after the actual outburst has died out. Conclusions: We propose that some of the multi-ringed structures seen in optical/infrared images of several protoplanetary disks may have their origin in outbursts that occurred decades or centuries ago.
△ Less
Submitted 8 September, 2018;
originally announced September 2018.
-
Shadows and asymmetries in the T Tauri disk HD 143006: Evidence for a misaligned inner disk
Authors:
M. Benisty,
A. Juhasz,
S. Facchini,
P. Pinilla,
J. de Boer,
L. M. Perez,
M. Keppler,
G. Muro-Arena,
M. Villenave,
S. Andrews,
C. Dominik,
C. P. Dullemond,
A. Gallenne,
A. Garufi,
C. Ginski,
A. Isella
Abstract:
While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, r…
▽ More
While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, reaching an angular resolution of ~0.037" (~6 au). We obtained two datasets, one with a 145 mas diameter coronagraph, and the other without, enabling us to probe the disk structure down to an angular separation of ~0.06" (~10 au). In our observations, the disk of HD 143006 is clearly resolved up to ~0.5" and shows a clear large-scale asymmetry with the eastern side brighter than the western side. We detect a number of additional features, including two gaps and a ring. The ring shows an overbrightness at a position angle (PA) of ~140 deg, extending over a range in position angle of ~60 deg, and two narrow dark regions. The two narrow dark lanes and the overall large-scale asymmetry are indicative of shadowing effects, likely due to a misaligned inner disk. We demonstrate the remarkable resemblance between the scattered light image of HD 143006 and a model prediction of a warped disk due to an inclined binary companion. The warped disk model, based on the hydrodynamic simulations combined with 3D radiative transfer calculations, reproduces all major morphological features. However, it does not account for the observed overbrightness at PA~140 deg. Shadows have been detected in several protoplanetary disks, suggesting that misalignment in disks is not uncommon. However, the origin of the misalignment is not clear. As-yet-undetected stellar or massive planetary companions could be responsible for them, and naturally account for the presence of depleted inner cavities.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
The Millimeter Continuum Size-Frequency Relationship in the UZ Tau E Disk
Authors:
Anjali Tripathi,
Sean M. Andrews,
Tilman Birnstiel,
Claire J. Chandler,
Andrea Isella,
Laura M. Perez,
Robert J. Harris,
Luca Ricci,
David J. Wilner,
John M. Carpenter,
Nuria Calvet,
Stuartt. A. Corder,
Adam T. Deller,
Cornelis P. Dullemond,
Jane S. Greaves,
Thomas Henning,
Woojin Kwon,
T. Joseph W. Lazio,
Hendrik Linz,
Leonardo Testi
Abstract:
We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence…
▽ More
We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto ν^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.
△ Less
Submitted 16 May, 2018;
originally announced May 2018.
-
Planetesimal formation during protoplanetary disk buildup
Authors:
Joanna Drazkowska,
Cornelis P. Dullemond
Abstract:
Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during…
▽ More
Models of dust coagulation and subsequent planetesimal formation are usually computed on the backdrop of an already fully formed protoplanetary disk model. At the same time, observational studies suggest that planetesimal formation should start early, possibly even before the protoplanetary disk is fully formed. In this paper, we investigate under which conditions planetesimals already form during the disk buildup stage, in which gas and dust fall onto the disk from its parent molecular cloud. We couple our earlier planetesimal formation model at the water snow line to a simple model of disk formation and evolution. We find that under most conditions planetesimals only form after the buildup stage when the disk becomes less massive and less hot. However, there are parameters for which planetesimals already form during the disk buildup. This occurs when the viscosity driving the disk evolution is intermediate ($α_v \sim 10^{-3}-10^{-2}$) while the turbulent mixing of the dust is reduced compared to that ($α_t \le 0.03 \cdot α_v$), and with the assumption that water vapor is vertically well-mixed with the gas. Such $α_t \ll α_v$ scenario could be expected for layered accretion, where the gas flow is mostly driven by the active surface layers, while the midplane layers, where most of the dust resides, are quiescent.
△ Less
Submitted 23 February, 2023; v1 submitted 1 March, 2018;
originally announced March 2018.
-
Impact splash chondrule formation during planetesimal recycling
Authors:
Tim Lichtenberg,
Gregor J. Golabek,
Cornelis P. Dullemond,
Maria Schönbächler,
Taras V. Gerya,
Michael R. Meyer
Abstract:
Chondrules are the dominant bulk silicate constituent of chondritic meteorites and originate from highly energetic, local processes during the first million years after the birth of the Sun. So far, an astrophysically consistent chondrule formation scenario, explaining major chemical, isotopic and textural features, remains elusive. Here, we examine the prospect of forming chondrules from planetes…
▽ More
Chondrules are the dominant bulk silicate constituent of chondritic meteorites and originate from highly energetic, local processes during the first million years after the birth of the Sun. So far, an astrophysically consistent chondrule formation scenario, explaining major chemical, isotopic and textural features, remains elusive. Here, we examine the prospect of forming chondrules from planetesimal collisions. We show that intensely melted bodies with interior magma oceans became rapidly chemically equilibrated and physically differentiated. Therefore, collisional interactions among such bodies would have resulted in chondrule-like but basaltic spherules, which are not observed in the meteoritic record. This inconsistency with the expected dynamical interactions hints at an incomplete understanding of the planetary growth regime during the protoplanetary disk phase. To resolve this conundrum, we examine how the observed chemical and isotopic features of chondrules constrain the dynamical environment of accreting chondrite parent bodies by interpreting the meteoritic record as an impact-generated proxy of planetesimals that underwent repeated collision and reaccretion cycles. Using a coupled evolution-collision model we demonstrate that the vast majority of collisional debris feeding the asteroid main belt must be derived from planetesimals which were partially molten at maximum. Therefore, the precursors of chondrite parent bodies either formed primarily small, from sub-canonical aluminum-26 reservoirs, or collisional destruction mechanisms were efficient enough to shatter planetesimals before they reached the magma ocean phase. Finally, we outline the window in parameter space for which chondrule formation from planetesimal collisions can be reconciled with the meteoritic record and how our results can be used to further constrain early solar system dynamics.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
Dust-driven viscous ring-instability in protoplanetary disks
Authors:
C. P. Dullemond,
A. B. T. Penzlin
Abstract:
Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. Th…
▽ More
Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. The presence of copious amounts of dust tends to strongly reduce the conductivity of the gas, thereby inhibiting the magneto-rotational instability, and thus reducing the turbulence in the disk. From viscous disk theory it is known that a disk tends to increase its surface density in regions where the viscosity (i.e. turbulence) is low. Local maxima in the gas pressure tend to attract dust through radial drift, increasing the dust content even more. We investigate mathematically if this could potentially lead to a feedback loop in which a perturbation in the dust surface density could perturb the gas surface density, leading to increased dust drift and thus amplification of the dust perturbation and, as a consequence, the gas perturbation. We find that this is indeed possible, even for moderately small dust grain sizes, which drift less efficiently, but which are more likely to affect the gas ionization degree. We speculate that this instability could be triggered by the small dust population initially, and when the local pressure maxima are strong enough, the larger dust grains get trapped and lead to the familiar ring-like shapes. We also discuss the many uncertainties and limitations of this model.
△ Less
Submitted 27 November, 2017; v1 submitted 1 September, 2017;
originally announced September 2017.
-
Efficiency of thermal relaxation by radiative processes in protoplanetary discs: constraints on hydrodynamic turbulence
Authors:
M. G. Malygin,
H. Klahr,
D. Semenov,
Th. Henning,
C. P. Dullemond
Abstract:
Hydrodynamic, non-magnetic instabilities can provide turbulent stress in the regions of protoplanetary discs, where the MRI can not develop. The induced motions influence the grain growth, from which formation of planetesimals begins. Thermal relaxation of the gas constrains origins of the identified hydrodynamic sources of turbulence in discs.
We estimate the radiative relaxation timescale of t…
▽ More
Hydrodynamic, non-magnetic instabilities can provide turbulent stress in the regions of protoplanetary discs, where the MRI can not develop. The induced motions influence the grain growth, from which formation of planetesimals begins. Thermal relaxation of the gas constrains origins of the identified hydrodynamic sources of turbulence in discs.
We estimate the radiative relaxation timescale of temperature perturbations and study the dependence of this timescale on the perturbation wavelength, the location within the disc, the disc mass, and the dust-to-gas mass ratio. We then apply thermal relaxation criteria to localise modes of the convective overstability, the vertical shear instability, and the zombie vortex instability.
Our calculations employed the latest tabulated dust and gas mean opacities and we account for the collisional coupling to the emitting species.
The relaxation criterion defines the bulk of a typical T Tauri disc as unstable to the development of linear hydrodynamic instabilities. The midplane is unstable to the convective overstability from at most $2\mbox{ au}$ and up to $40\mbox{ au}$, as well as beyond $140\mbox{ au}$. The vertical shear instability can develop between $15\mbox{ au}$ and $180\mbox{ au}$. The successive generation of (zombie) vortices from a seeded noise can work within the inner $0{.}8\mbox{ au}$.
Dynamic disc modelling with the evolution of dust and gas opacities is required to clearly localise the hydrodynamic turbulence, and especially its non-linear phase.
△ Less
Submitted 22 April, 2017;
originally announced April 2017.
-
Mid-infrared interferometric variability of DG Tau: implications for the inner-disk structure
Authors:
J. Varga,
K. É. Gabányi,
P. Ábrahám,
L. Chen,
Á. Kóspál,
J. Menu,
Th. Ratzka,
R. van Boekel,
C. P. Dullemond,
Th. Henning,
W. Jaffe,
A. Juhász,
A. Moór,
L. Mosoni,
N. Sipos
Abstract:
Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $μ$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. W…
▽ More
Context. DG Tau is a low-mass pre-main sequence star, whose strongly accreting protoplanetary disk exhibits a so-far enigmatic behavior: its mid-infrared thermal emission is strongly time-variable, even turning the 10 $μ$m silicate feature from emission to absorption temporarily. Aims. We look for the reason for the spectral variability at high spatial resolution and at multiple epochs. Methods. We study the temporal variability of the mid-infrared interferometric signal, observed with the VLTI/MIDI instrument at six epochs between 2011 and 2014. We fit a geometric disk model to the observed interferometric signal to obtain spatial information about the disk. We also model the mid-infrared spectra by template fitting to characterize the profile and time dependence of the silicate emission. We use physically motivated radiative transfer modeling to interpret the mid-infrared interferometric spectra. Results. The inner disk (r<1-3 au) spectra exhibit a 10 $μ$m absorption feature related to amorphous silicate grains. The outer disk (r>1-3 au) spectra show a crystalline silicate feature in emission, similar to the spectra of comet Hale-Bopp. The striking difference between the inner and outer disk spectral feature is highly unusual among T Tauri stars. The mid-infrared variability is dominated by the outer disk. The strength of the silicate feature changed by more than a factor of two. Between 2011 and 2014 the half-light radius of the mid-infrared-emitting region decreased from 1.15 to 0.7 au. Conclusions. For the origin of the absorption we discuss four possible explanations: a cold obscuring envelope, an accretion heated inner disk, a temperature inversion on the disk surface and a misaligned inner geometry. The silicate emission in the outer disk can be explained by dusty material high above the disk plane, whose mass can change with time, possibly due to turbulence in the disk.
△ Less
Submitted 19 April, 2017;
originally announced April 2017.
-
Redistribution of CO at the Location of the CO Ice Line in evolving Gas and Dust Disks
Authors:
Sebastian Markus Stammler,
Tilman Birnstiel,
Olja Panić,
Cornelis Petrus Dullemond,
Carsten Dominik
Abstract:
Context. Ice lines are suggested to play a significant role in grain growth and planetesimal formation in protoplanetary disks. Evaporation fronts directly influence the gas and ice abundances of volatile species in the disk and therefore the coagulation physics and efficiency and the chemical composition of the resulting planetesimals.
Aims. In this work we investigate the influence of the exis…
▽ More
Context. Ice lines are suggested to play a significant role in grain growth and planetesimal formation in protoplanetary disks. Evaporation fronts directly influence the gas and ice abundances of volatile species in the disk and therefore the coagulation physics and efficiency and the chemical composition of the resulting planetesimals.
Aims. In this work we investigate the influence of the existence of the CO ice line on the particle growth and on the distribution of CO in the disk.
Methods. We include the possibility of tracking the CO content and/or other volatiles in particles and in the gas in our existing dust coagulation and disk evolution model and developed a method for evaporation and condensation of CO using the Hertz-Knudsen equation. Our model does not include fragmentation, yet, which will be part of further investigations.
Results. We find no enhanced grain growth just outside the ice line where the particle size is limited by radial drift. Instead we find a depletion of solid material inside the ice line which is solely due to evaporation of the CO. Such a depression inside the ice line may be observable and may help to quantify the processes described in this work. Furthermore, we find that the viscosity and diffusivity of the gas heavily influence the re-distribution of vaporized CO at the ice line and can lead to an increase in the CO abundance by up to a factors of a few in the region just inside the ice line. Depending on the strength of the gaseous transport mechanisms the position of the ice line in our model can change by up to 10 AU and consequently, the temperature at that location can range from 21 K to 23 K.
△ Less
Submitted 31 January, 2017; v1 submitted 9 January, 2017;
originally announced January 2017.