Showing 1–50 of 80 results for author: Flock, M

From traffic jams to roadblocks: The outer regions of TW Hya with ALMA Band 8

Authors: S. Das, N. T. Kurtovic, M. Flock

Abstract: We present Atacama Large Millimeter Array (ALMA) observations of TW Hya at 0.65 mm with 0.5 arcsecond angular resolution, together with high angular resolution archival observations at 0.87 mm, 1.3 mm, 2.1 mm and 3.1 mm. We constrain the outer disk emission with both image-plane retrieval, and visibility-plane modeling with non-parametric and parametric fitting tools. Our results confirm emission… ▽ More We present Atacama Large Millimeter Array (ALMA) observations of TW Hya at 0.65 mm with 0.5 arcsecond angular resolution, together with high angular resolution archival observations at 0.87 mm, 1.3 mm, 2.1 mm and 3.1 mm. We constrain the outer disk emission with both image-plane retrieval, and visibility-plane modeling with non-parametric and parametric fitting tools. Our results confirm emission in the outer disk regions of TW Hya (60 au < R < 110 au) at 0.65 mm, 0.87 mm and 1.3 mm. With image-plane retrieval, we resolve the new continuum gap and ring, namely D79 and B86, at 0.87 mm and 1.3 mm. With visibility-plane modeling, we also detect this substructure at 0.65 mm, and it is consistent in location with the outer ring proposed by Ilee et al. 2022 at the 1-sigma level. Furthermore, it has a high spectral index of 3.7, which may indicate dust grain sizes << 1 mm. It may be a dust trap or a traffic jam, that has a flux density of 60 mJy and a mass (1.59 earth masses) that accounts for up to 2% of the dust disk at 0.65 mm. In conclusion, we confirm the existence of a faint ring in the outer regions of TW Hya at multiple millimeter wavelengths. With visibility-plane modeling, we are able to set constrains that are 3 times better than the resolution of our Band 8 observations. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. The associated data can be found in: https://doi.org/10.5281/zenodo.12677551

Support for fragile porous dust in a gravitationally self-regulated disk around IM Lup

Authors: Takahiro Ueda, Ryo Tazaki, Satoshi Okuzumi, Mario Flock, Prakruti Sudarshan

Abstract: Protoplanetary disks, the birthplace of planets, are expected to be gravitationally unstable in their early phase of evolution. IM Lup, a well-known T-Tauri star, is surrounded by a protoplanetary disk with spiral arms likely caused by gravitational instability. The IM Lup disk has been observed using various methods, but developing a unified explanatory model is challenging. Here we present a phy… ▽ More Protoplanetary disks, the birthplace of planets, are expected to be gravitationally unstable in their early phase of evolution. IM Lup, a well-known T-Tauri star, is surrounded by a protoplanetary disk with spiral arms likely caused by gravitational instability. The IM Lup disk has been observed using various methods, but developing a unified explanatory model is challenging. Here we present a physical model of the IM Lup disk that offers a comprehensive explanation for diverse observations spanning from near-infrared to millimeter wavelengths. Our findings underscore the importance of dust fragility in retaining the observed millimeter emission and reveal the preference for moderately porous dust to explain observed millimeter polarization. We also find that the inner disk region is likely heated by gas accretion, providing a natural explanation for bright millimeter emission within 20 au. The actively heated inner region in the model casts a 100-au-scale shadow, aligning seamlessly with the near-infrared scattered light observation. The presence of accretion heating also supports the fragile dust scenario in which accretion efficiently heat the disk midplane. Due to the fragility of dust, it is unlikely that a potential embedded planet at 100 au formed via pebble accretion in a smooth disk, pointing to local dust enhancement boosting pebble accretion or alternative pathways such as outward migration or gravitational fragmentation. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted for publication. Author version

Mixing is easy: New insights for cosmochemical evolution from pre-stellar core collapse

Authors: Asmita Bhandare, Benoît Commerçon, Guillaume Laibe, Mario Flock, Rolf Kuiper, Thomas Henning, Andrea Mignone, Gabriel-Dominique Marleau

Abstract: Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. W… ▽ More Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot and high-pressure regions formed during gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas-dust interaction and resolve zones of dust concentration during star-disc formation. We explore if thermal and dynamical conditions developed during disc formation can generate gas flows that efficiently mix and transport well-coupled gas and dust components. We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a 2D axisymmetric geometry and follow the azimuthal component of velocity. Dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10 and 100 micron-sized neutral spherical dust. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 Msun and explore effects of initial rotation rates and cloud core sizes. Our study underlines mechanisms for early transport of dust from inner hot disc regions via the occurrence of meridional flows and outflow. The vortical flow fosters dynamical mixing and retention of dust while thermal pressure driven outflow replenishes dust in the outer disc. Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse. [Abridged] △ Less

Submitted 12 May, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: Language edited version

The inner disk rim of HD 163296: linking radiative hydrostatic models with infrared interferometry

Authors: Ondřej Chrenko, Mario Flock, Takahiro Ueda, Antoine Mérand, Myriam Benisty, Raúl O. Chametla

Abstract: Previous studies of the protoplanetary disk HD 163296 revealed that the morphology of its sub-au infrared emission encompasses the terminal sublimation front of dust grains, referred to as the inner rim, but also extends into the (supposedly) dust-free region within it. Here, we present a set of radiative hydrostatic simulations of the inner rim in order to assess how much the rim alone can contri… ▽ More Previous studies of the protoplanetary disk HD 163296 revealed that the morphology of its sub-au infrared emission encompasses the terminal sublimation front of dust grains, referred to as the inner rim, but also extends into the (supposedly) dust-free region within it. Here, we present a set of radiative hydrostatic simulations of the inner rim in order to assess how much the rim alone can contribute to the observed interferometric visibilities $V$, half-light radii $R_{\mathrm{hl}}$, and fractional disk fluxes $\mathcal{F}$ in the wavelength range $1.5$--$13\,μ\mathrm{m}$. In our set of models, we regulate the cooling efficiency of the disk via the boundary condition for radiation diffusion and we also modify the shape of the sublimation front. We find that when the cooling efficiency is reduced, the infrared photosphere at the rim becomes hotter, leading to an increase of $R_{\mathrm{hl}}$ sufficient to match the observations. However, the near-infrared disk flux is typically too low ($\mathcal{F}\simeq0.25$ at $1.5\,μ\mathrm{m}$), resulting in H-band visibility curves located above the observed data. We show that the match to the H-band observations up to moderate baselines can be improved when a wall-shaped rather than curved sublimation front is considered. Nevertheless, our model visibilities always exhibit a bounce at long baselines, which is not observed, confirming the need for additional emission interior to the rim. In summary, our study illustrates how the temperature structure and geometry of the inner rim needs to change in order to boost the rim's infrared emission. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted for publication in AJ, 20 pages, 10 figures

The GRAVITY young stellar object survey: XI. Imaging the hot gas emission around the Herbig Ae star HD 58647

Authors: Y. -I. Bouarour, R. Garcia Lopez, J. Sanchez-Bermudez, A. Caratti o Garatti, K. Perraut, N. Aimar, A. Amorim, J. -P. Berger, G. Bourdarot, W. Brandner, Y. Clénet, P. T. de Zeeuw, C. Dougados, A. Drescher, A. Eckart, F. Eisenhauer, M. Flock, P. Garcia, E. Gendron, R. Genzel, S. Gillessen, S. Grant, G. Heißel, Th. Henning, L. Jocou , et al. (23 additional authors not shown)

Abstract: We aim to investigate the origin of the HI Br$γ$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$γ$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of… ▽ More We aim to investigate the origin of the HI Br$γ$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$γ$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of HD58647 with both high spectral and high angular resolution. Thanks to the extensive $uv$ coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Br$γ$-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Br$γ$-emitting region. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14\degr. The size of the continuum emission at 5-sigma levels is $\sim$ 1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas $\pm$ 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD58647 traced by the Br$γ$ ,has been imaged. This allowed us to constrain the size of the Br$γ$-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas $\pm$ 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: accepted in A&A on 20/11/2023

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks II. Secondary instabilities and stability regions

Authors: Julio David Melon Fuksman, Mario Flock, Hubert Klahr

Abstract: The vertical shear instability (VSI) is a hydrodynamical instability likely to produce turbulence in the dead zones of protoplanetary disks. Various aspects of this instability remain to be understood, including the disk regions where it can operate and the physical phenomena leading to its saturation. In this work, we studied the growth and evolution of secondary instabilities parasitic to the VS… ▽ More The vertical shear instability (VSI) is a hydrodynamical instability likely to produce turbulence in the dead zones of protoplanetary disks. Various aspects of this instability remain to be understood, including the disk regions where it can operate and the physical phenomena leading to its saturation. In this work, we studied the growth and evolution of secondary instabilities parasitic to the VSI, examining their relation with its saturation in axisymmetric radiation-hydrodynamical simulations of protoplanetary disks. We also constructed stability maps for our disk models, considering temperature stratifications enforced by stellar irradiation and radiative cooling and incorporating the effects of dust-gas collisions and molecular line emission. We found that the flow pattern produced by the interplay of the axisymmetric VSI modes and the baroclinic torque forms bands of nearly uniform specific angular momentum. In the high-shear regions in between these bands, the Kelvin-Helmholtz instability (KHI) is triggered. The significant transfer of kinetic energy to small-scale eddies produced by the KHI and possibly even the baroclinic acceleration of eddies limit the maximum energy of the VSI modes, likely leading to the saturation of the VSI. A third instability mechanism, consisting of an amplification of eddies by baroclinic torques, forms meridional vortices with Mach numbers up to $\sim 0.4$. Our stability analysis suggests that protoplanetary disks can be VSI-unstable in surface layers up to tens of au for reasonably high gas emissivities, even in regions where the midplane is stable. This picture is consistent with current observations of disks showing thin midplane millimeter-sized dust layers while appearing vertically extended in optical and near-infrared wavelengths. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Vertical shear instability in two-moment radiation-hydrodynamical simulations of irradiated protoplanetary disks I. Angular momentum transport and turbulent heating

Authors: Julio David Melon Fuksman, Mario Flock, Hubert Klahr

Abstract: We studied the linear and nonlinear evolution of the Vertical Shear Instability (VSI) in axisymmetric models of protoplanetary disks, focusing on the transport of angular momentum, the produced temperature perturbations, and the applicability of local stability conditions. We modeled the gas-dust mixture via high-resolution two-moment (M1) radiation-hydrodynamical simulations including stellar irr… ▽ More We studied the linear and nonlinear evolution of the Vertical Shear Instability (VSI) in axisymmetric models of protoplanetary disks, focusing on the transport of angular momentum, the produced temperature perturbations, and the applicability of local stability conditions. We modeled the gas-dust mixture via high-resolution two-moment (M1) radiation-hydrodynamical simulations including stellar irradiation with frequency-dependent opacities. We found that, given sufficient depletion of small grains (with a dust-to-gas mass ratio of $10\%$ of our nominal value of $10^{-3}$ for $<0.25$ $μ$m grains), the VSI can operate in surface disk layers while being inactive close to the midplane, resulting in a suppression of the VSI body modes. The VSI reduces the initial vertical shear in bands of approximately uniform specific angular momentum, whose formation is likely favored by the enforced axisymmetry. Similarities with Reynolds stresses and angular momentum distributions in 3D simulations suggest that the VSI-induced angular momentum mixing in the radial direction may be predominantly axisymmetric. The stability regions in our models are well explained by local stability criteria, while the employment of global criteria is still justifiable up to a few scale heights above the midplane, at least as long as VSI modes are radially optically thin. Turbulent heating produces only marginal temperature increases of at most $0.1\%$ and $0.01\%$ in the nominal and dust-depleted models, respectively, peaking at a few (approximately three) scale heights above the midplane. We conclude that it is unlikely that the VSI can, in general, lead to any significant temperature increase since that would either require it to efficiently operate in largely optically thick disk regions or to produce larger levels of turbulence than predicted by models of passive irradiated disks. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Kinematic signatures of planet-disk interactions in VSI-turbulent protoplanetary disks

Authors: Marcelo Barraza-Alfaro, Mario Flock, Thomas Henning

Abstract: Context. Planets are thought to form inside weakly ionized regions of protoplanetary disks, where turbulence creates ideal conditions for solid growth. However, the nature of this turbulence is still uncertain. In this zone, vertical shear instability (VSI) can operate, inducing a low level of gas turbulence and large-scale motions. Resolving kinematic signatures of VSI may reveal the origin of tu… ▽ More Context. Planets are thought to form inside weakly ionized regions of protoplanetary disks, where turbulence creates ideal conditions for solid growth. However, the nature of this turbulence is still uncertain. In this zone, vertical shear instability (VSI) can operate, inducing a low level of gas turbulence and large-scale motions. Resolving kinematic signatures of VSI may reveal the origin of turbulence in planet-forming disks. However, an exploration of kinematic signatures of the interplay between VSI and forming planets is needed for a correct interpretation of radio interferometric observations. Robust detection of VSI would lead to a deeper understanding of the impact of gas turbulence on planet formation. Aims. The goal of this study is to explore the effect of VSI on the disk substructures triggered by an embedded massive planet. We focus on the impact of this interplay on CO kinematic observations with ALMA. Methods. We conduct global 3D hydrodynamical simulations of VSI-unstable disks with and without embedded massive planets, exploring Saturn- and Jupiter-mass cases. We study the effect of planets on the VSI gas dynamics, comparing with viscous disks. Post-processing the simulations with a radiative transfer code, we examine the kinematic signatures expected in CO molecular line emission, varying disk inclination. Further, we simulate ALMA high-resolution observations to test the observability of VSI and planetary signatures. Results. The embedded planet dampens the VSI along a radial region, most effective at the disk midplane. For the Saturn case, the VSI modes are distorted by the planet's spirals producing mixed kinematic signatures. For the Jupiter case, the planet's influence dominates the disk gas kinematics. Conclusions. The presence of massive embedded planets can weaken the VSI large-scale gas flows, limiting its observability in CO kinematic observations. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 27 pages, 17 figures and 2 tables

Probing the Temperature Structure of the Inner Region of a Protoplanetary Disk

Authors: Takahiro Ueda, Satoshi Okuzumi, Akimasa Kataoka, Mario Flock

Abstract: Midplane heating induced by disk accretion plays a key role in determining the disk temperature particularly at the inner disk midplane where planets form. However, the efficiency of accretion heating has been not well constrained by observations. We construct two-dimensional models of the Class II disk around CW Tau, taking into account the midplane heating. The models are compared with the ALMA… ▽ More Midplane heating induced by disk accretion plays a key role in determining the disk temperature particularly at the inner disk midplane where planets form. However, the efficiency of accretion heating has been not well constrained by observations. We construct two-dimensional models of the Class II disk around CW Tau, taking into account the midplane heating. The models are compared with the ALMA dust continuum observations at Bands 4, 6, 7 and 8, with an angular resolution of 0.1 arcsec. The observed brightness temperatures are almost wavelength-indenpendent at $\lesssim$10 au. We find that if the maximum dust size $a_{\rm max}$ is $\lesssim100~{\rm μm}$, the brightness temperatures predicted by the model exceed the observed values, regardless of the efficiency of accretion heating. The low observed brightness temperatures can be explained if millimeter scattering reduces the intensity. If the disk is passive, $a_{\rm max}$ needs to be either $\sim150~{\rm μm}$ or $\gtrsim$ few ${\rm cm}$. The accretion heating significantly increases the brightness temperature particularly when $a_{\rm max}\lesssim300~{\rm μm}$, and hence $a_{\rm max}$ needs to be either $\sim300~{\rm μm}$ or $\gtrsim$ few ${\rm cm}$. The midplane temperature is expected to be $\sim$1.5-3 times higher than the observed brightness temperatures, depending on the models. The dust settling effectively increases the temperature of the dust responsible for the millimeter emission in the active disk, which makes the model with $300~{\rm μm}$-sized dust overpredicts the brightness temperatures when strong turbulence is absent. Porous dust (porosity of 0.9) makes the accretion heating more efficient so that some sort of reduction in accretion heating is required. Future longer wavelength and higher angular resolution observations will help us constrain the heating mechanisms of the inner protoplanetary disks. △ Less

Submitted 21 May, 2023; originally announced May 2023.

Comments: 19 pages, 25 figures, accepted for publication in A&A

Journal ref: A&A 675, A176 (2023)

Meridional Circulation driven by Planetary Spiral Wakes in Radiative and Magnetized Protoplanetary Discs

Authors: Marco Cilibrasi, Mario Flock, Judit Szulágyi

Abstract: We study a Jupiter-mass planet formation for the first time in radiative magneto-hydrodynamics (MHD) simulations and compare it with pure hydrodynamical simulations, as well as to different isothermal configurations. We found that the meridional circulation is the same in every setup. The planetary spiral wakes drive a vertical stirring inside the protoplanetary disc and the encounter with these s… ▽ More We study a Jupiter-mass planet formation for the first time in radiative magneto-hydrodynamics (MHD) simulations and compare it with pure hydrodynamical simulations, as well as to different isothermal configurations. We found that the meridional circulation is the same in every setup. The planetary spiral wakes drive a vertical stirring inside the protoplanetary disc and the encounter with these shock fronts also helps in delivering gas vertically onto the Hill-sphere. The accretion dynamics are unchanged: the planet accretes vertically, and there is outflow in the midplane regions inside the Hill-sphere. We determined the effective $α$-viscosity generated in the disc by the various angular momentum loss mechanisms, which showed that magnetic fields produce high turbulence in the ideal MHD limit, that grows from $α\sim 10^{-2.5}$ up to $\sim 10^{-1.5}$ after the planet spirals develop. In the HD simulations, the planetary spirals contribute to $α\sim 10^{-3}$, making this a very important angular momentum transport mechanism. Due to the various $α$ values in the different setups, the gap opening is different in each case. In the radiative MHD setups, the high turbulent viscosity prevents gap opening, leading to a higher Hill mass, and no clear dust trapping regions. While the Hill accretion rate is $10^{-6} \rm{M_{Jup}/yr}$ in all setups, the accretion variability is orders of magnitude higher in radiative runs than in isothermal ones. Finally, with higher-resolution runs, the magneto-rotational instability started to be resolved, changing the effective viscosity and increasing the heating in the disc. △ Less

Submitted 13 May, 2023; originally announced May 2023.

The impact of dust evolution on the dead zone outer edge in magnetized protoplanetary disks

Authors: Timmy N. Delage, Matías Gárate, Satoshi Okuzumi, Chao-Chin Yang, Paola Pinilla, Mario Flock, Sebastian Markus Stammler, Tilman Birnstiel

Abstract: [Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the… ▽ More [Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles. Since the dust content of the disk decreases over million years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced $α$-parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there is enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once it is no longer the case, it is expected to be controlled by gas evolution and occurs on a viscous evolution timescale. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 23 pages, 13 figures, Accepted for publication in A&A

Journal ref: A&A 674, A190 (2023)

Forbidden emission lines in protostellar outflows and jets with MUSE

Authors: Lizxandra Flores-Rivera, Mario Flock, Nicolás Kurtovic, Bernd Husemann, Andrea Banzatti, Simon C. Ringqvist, Sebastian Kamann, André Müller, Christian Fendt, Rebeca Garcia Lopez, Gabriel-Dominique Marleau, Thomas Henning, Carlos Carrasco-Gonzalez, Roy van Boekel, Miriam Keppler, Ralf Launhardt, Yuhiko Aoyama

Abstract: Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral fiel… ▽ More Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral field spectrograph of the Multi Unit Spectroscopic Explorer (MUSE), at the Very Large Telescope (VLT). With a resolution of 0.025 X 0.025 arcsec$^{2}$, we aim to derive the position angle of the outflow/jet (PA$_{outflow/jet}$) that is connected with the inner disk. The forbidden emission lines analyzed here have their origin at the inner parts of the protoplanetary disk. From the maximum intensity emission along the outflow/jet in DL Tau, CI Tau, DS Tau, IP Tau, and IM Lup, we were able to reliably measure the PA$_{outflow/jet}$ for most of the identified lines. We found that our estimates agree with PA$_{dust}$ for most of the disks. These estimates depend on the signal-to-noise level and the collimation of the outflow (jet). The outflows/jets in CIDA 9, GO Tau, and GW Lup are too compact for a PA$_{outflow/jet}$ to be estimated. Based on our kinematics analysis, we confirm that DL Tau and CI Tau host a strong outflow/jet with line-of-sight velocities much greater than 100 km s$^{-1}$, whereas DS Tau, IP Tau, and IM Lup velocities are lower and their structures encompass low-velocity components to be more associated with winds. Our estimates for the mass-loss rate, $\dot{M}_{loss}$, range between (1.1-6.5)X10$^{-7}$-10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for the disk-outflow/jet systems analyzed here. The outflow/jet systems analyzed here are aligned within around 1 degree between the inner and outer disk. Further observations are needed to confirm a potential misalignment in IM Lup. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted on December 19 2022

On the underestimation of dust mass in protoplanetary disks: Effects of disk structure and dust properties

Authors: Yao Liu, Hendrik Linz, Min Fang, Thomas Henning, Sebastian Wolf, Mario Flock, Giovanni P. Rosotti, Hongchi Wang, Dafa Li

Abstract: The total amount of dust grains in protoplanetary disks is one of the key properties that characterize the potential for planet formation. With (sub-)millimeter flux measurements, literature studies usually derive the dust mass using an analytic form under the assumption of optically thin emission, which may lead to substantial underestimation. In this work, we conduct a parameter study with the g… ▽ More The total amount of dust grains in protoplanetary disks is one of the key properties that characterize the potential for planet formation. With (sub-)millimeter flux measurements, literature studies usually derive the dust mass using an analytic form under the assumption of optically thin emission, which may lead to substantial underestimation. In this work, we conduct a parameter study with the goal of investigating the effects of disk structure and dust properties on the underestimation through self-consistent radiative transfer models. Different dust models, scattering modes and approaches for dust settling are considered and compared. The influences of disk substructures, such as rings and crescents, on the mass derivation are investigated as well. The results indicate that the traditional analytic method can underestimate the mass by a factor of a few to hundreds, depending on the optical depth along the line of sight set mainly by the true dust mass, disk size and inclination. As an application, we perform a detailed radiative transfer modeling of the spectral energy distribution of DoAr 33, one of the observed DSHARP disks. When the DSHARP dust opacities are adopted, the most probable dust mass returned from the Bayesian analysis is roughly 7 times higher than the value given by the analytic calculation. Our study demonstrates that estimating disk dust masses from radiative transfer modeling is one solution for alleviating the problem of insufficient mass for planet formation raised in the ALMA era. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: 16 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 668, A175 (2022)

Trapping (sub-)Neptunes similar to TOI-216b at the inner disk rim: Implications for the disk viscosity and the Neptunian desert

Authors: Ondřej Chrenko, Raúl O. Chametla, David Nesvorný, Mario Flock

Abstract: [Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive co-rotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, lose the support of the co-rotation torque as their co-rotation regions… ▽ More [Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive co-rotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, lose the support of the co-rotation torque as their co-rotation regions become depleted, and the trapping efficiency then becomes uncertain. We study what it takes to trap such gap-opening planets at the inner disk rim. We performed 2D locally isothermal and non-isothermal hydrodynamic simulations of planet migration. A viscosity transition was introduced in the disk to (i) create a density drop and (ii) mimic the viscosity increase as the planet migrated from a dead zone towards a region with active magneto-rotational instability (MRI). We chose TOI-216b as a Neptune-like upper-limit test case, but we also explored different planetary masses, both on fixed and evolving orbits. For planet-to-star mass ratios $q\simeq(4$-$8)\times10^{-5}$, the density drop at the disk rim becomes reshaped due to a gap opening and is often replaced with a small density bump centered on the planet's corotation. Trapping is possible only if the bump retains enough gas mass and if the co-rotation region becomes azimuthally asymmetric, with an island of librating streamlines that accumulate a gas overdensity ahead of the planet. The overdensity exerts a positive torque that can counteract the negative torque of spiral arms. In our model, efficient trapping depends on the $α$ viscosity and its contrast across the viscosity transition. In order to trap TOI-216b, $α_{\mathrm{DZ}}=10^{-3}$ in the dead zone requires $α_{\mathrm{MRI}}\gtrsim5\times10^{-2}$ in the MRI-active zone. If $α_{\mathrm{DZ}}=5\times10^{-4}$, $α_{\mathrm{MRI}}\gtrsim7.5\times10^{-2}$ is needed. △ Less

Submitted 30 August, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A, 21 pages, 16 figures; v2 with minor text corrections

Journal ref: A&A 666, A63 (2022)

Millimeter Gap Contrast as a Probe for Turbulence Level in Protoplanetary Disks

Authors: Yao Liu, Gesa H. -M. Bertrang, Mario Flock, Giovanni P. Rosotti, Ewine F. van Dishoeck, Yann Boehler, Stefano Facchini, Can Cui, Sebastian Wolf, Min Fang

Abstract: Turbulent motions are believed to regulate angular momentum transport and influence dust evolution in protoplanetary disks. Measuring the strength of turbulence is challenging through gas line observations because of the requirement for high spatial and spectral resolution data, and an exquisite determination of the temperature. In this work, taking the well-known HD 163296 disk as an example, we… ▽ More Turbulent motions are believed to regulate angular momentum transport and influence dust evolution in protoplanetary disks. Measuring the strength of turbulence is challenging through gas line observations because of the requirement for high spatial and spectral resolution data, and an exquisite determination of the temperature. In this work, taking the well-known HD 163296 disk as an example, we investigated the contrast of gaps identified in high angular resolution continuum images as a probe for the level of turbulence. With self-consistent radiative transfer models, we simultaneously analyzed the radial brightness profiles along the disk major and minor axes, and the azimuthal brightness profiles of the B67 and B100 rings. By fitting all the gap contrasts measured from these profiles, we constrained the gas-to-dust scale height ratio $Λ$ to be $3.0_{-0.8}^{+0.3}$, $1.2_{-0.1}^{+0.1}$ and ${\ge}\,6.5$ for the D48, B67 and B100 regions, respectively. The varying gas-to-dust scale height ratios indicate that the degree of dust settling changes with radius. The inferred values for $Λ$ translate into a turbulence level of $α_{\rm turb}\,{<}\,3\times10^{-3}$ in the D48 and B100 regions, which is consistent with previous upper limits set by gas line observations. However, turbulent motions in the B67 ring are strong with $α_{\rm turb}\,{\sim}1.2\,{\times}\,10^{-2}$. Due to the degeneracy between $Λ$ and the depth of dust surface density drops, the turbulence strength in the D86 gap region is not constrained. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 13 pages, 7 figures, accepted for publication in Science China Physics, Mechanics & Astronomy

Journal ref: Sci. China-Phys. Mech. Astron. 65, 129511 (2022)

Modeling the arc and ring structures in the HD 143006 disk

Authors: Yao Liu, Mario Flock, Min Fang

Abstract: Rings and asymmetries in protoplanetary disks are considered as signposts of ongoing planet formation. In this work, we conduct three-dimensional radiative transfer simulations to model the intriguing disk around HD 143006 that has three dust rings and a bright arc. A complex geometric configuration, with a misaligned inner disk, is assumed to account for the asymmetric structures. The two-dimensi… ▽ More Rings and asymmetries in protoplanetary disks are considered as signposts of ongoing planet formation. In this work, we conduct three-dimensional radiative transfer simulations to model the intriguing disk around HD 143006 that has three dust rings and a bright arc. A complex geometric configuration, with a misaligned inner disk, is assumed to account for the asymmetric structures. The two-dimensional surface density is constructed by iteratively fitting the ALMA data. We find that the dust temperature displays a notable discontinuity at the boundary of the misalignment. The ring masses range from 0.6 to 16Mearth that are systematically lower than those inferred in the younger HL Tau disk. The arc occupies nearly 20% of the total dust mass. Such a high mass fraction of dust grains concentrated in a local region is consistent with the mechanism of dust trapping into vortices. Assuming a gas-to-dust mass ratio of 30 that is constant throughout the disk, the dense and cold arc is close to the threshold of being gravitationally unstable, with the Toomre parameter Q~1.3. Nevertheless, our estimate of Q relies on the assumption for the unknown gas-to-dust mass ratio. Adopting a lower gas-to-dust mass ratio would increase the inferred Q value. Follow-up high resolution observations of dust and gas lines are needed to clarify the origin of the substructures. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 12 Pages, 8 figures, Accepted for publication in Science China Physics, Mechanics & Astronomy

Journal ref: Sci. China-Phys. Mech. Astron. 65, 269511 (2022)

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

Authors: G. Lesur, B. Ercolano, M. Flock, M. -K. Lin, C. -C. Yang, J. A. Barranco, P. Benitez-Llambay, J. Goodman, A. Johansen, H. Klahr, G. Laibe, W. Lyra, P. Marcus, R. P. Nelson, J. Squire, J. B. Simon, N. Turner, O. M. Umurhan, A. N. Youdin

Abstract: The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t… ▽ More The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Probing Inner-Edge of Dead Zones in Protoplanetary Disks with ALMA and Next Generation Very Large Array

Authors: Takahiro Ueda, Luca Ricci, Mario Flock, Zachary Castro

Abstract: The discovery of substructures in protoplanetary disks with ALMA has provided us key insights on the formation of planets. However, observational constraints on the formation of rocky planets have been still sparse, especially because of the limited spatial resolution. The inner edge of so-called dead zone is one of the preferential sites of rocky planet formation. We investigate the capabilities… ▽ More The discovery of substructures in protoplanetary disks with ALMA has provided us key insights on the formation of planets. However, observational constraints on the formation of rocky planets have been still sparse, especially because of the limited spatial resolution. The inner edge of so-called dead zone is one of the preferential sites of rocky planet formation. We investigate the capabilities of ALMA and ngVLA for observing a dust concentration expected at the inner edge of the dead-zone around a Herbig star. Herbig Ae/Be stars are useful laboratories for exploring the evolution of rocky grains in protoplanetary disks because of their high luminosity which pushes the dead-zone inner edge outward. We find that, thanks to its unprecedented angular resolution and sensitivity, ngVLA can detect the dust concentration at the dead-zone inner edge, with a reasonable integration time of 10 hrs at $λ=3,7$ mm and 1 cm. The dust concentration is expected to be optically thick at the ALMA wavelengths and cannot be spatially resolved due to its limited resolution. On the other hand, the flux density from the inner disk regions ($\sim$3--4 au) observed with current VLA is higher for disks with a dust ring, and hence would be a useful indicator that help us choose potential candidates of disks having a dust concentration at the inner most region. With these observations we can characterize the process of dust concentration in the innermost disk regions, where rocky planets can form. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 9 pages, 7 figures, accepted for publication in ApJ

TOI-216: Resonant Constraints on Planet Migration

Authors: David Nesvorny, Ondrej Chrenko, Mario Flock

Abstract: TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_ψ\simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit.… ▽ More TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_ψ\simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit. We carry out $N$-body simulations of planet migration in a protoplanetary gas disk to explain the orbital configuration of TOI-216 planets. We find that TOI-216b's migration must have been halted near its current orbital radius to allow for a convergent migration of the two planets into the resonance. For the inferred damping-to-migration timescale ratio $τ_e/τ_a \simeq 0.02$, overstable librations in the resonance lead to a limit cycle with $A_ψ\simeq 80^\circ$ and $e_{\rm b}<0.1$. The system could have remained in this configuration for the greater part of the protoplanetary disk lifetime. If the gas disk was removed from inside out, this would have reduced the libration amplitude to $A_ψ\simeq 60^\circ$ and boosted $e_{\rm b}$ via the resonant interaction with TOI-216c. Our results suggest a relatively fast inner disk removal ($\sim 10^5$ yr). Another means of explaining the large libration amplitude is stochastic stirring from a (turbulent) gas disk. For that to work, overstable librations would need to be suppressed, $τ_e/τ_a \simeq 0.05$, and very strong turbulent stirring (or some other source of large stochastic forcing) would need to overcome the damping effects of gas. Hydrodynamical simulations can be performed to test these models. △ Less

Submitted 18 October, 2021; originally announced October 2021.

Comments: ApJ

Steady-state accretion in magnetized protoplanetary disks

Authors: Timmy N. Delage, Satoshi Okuzumi, Mario Flock, Paola Pinilla, Natalia Dzyurkevich

Abstract: [abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $α$ is determined self-consistently… ▽ More [abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $α$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Comments: Accepted for publication in A&A

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

The Substructures in Disks undergoing Vertical Shear Instability: II. Observational Predictions for the Dust Continuum

Authors: Diana Blanco, Luca Ricci, Mario Flock, Neal Turner

Abstract: High-angular resolution observations at sub-millimeter/millimeter wavelengths of disks surrounding young stars have shown that their morphology is made of azimuthally-symmetric or point-symmetric substructures, in some cases with spiral arms, localized spur- or crescent-shaped features. The majority of theoretical studies with the aim of interpreting the observational results have focused on disk… ▽ More High-angular resolution observations at sub-millimeter/millimeter wavelengths of disks surrounding young stars have shown that their morphology is made of azimuthally-symmetric or point-symmetric substructures, in some cases with spiral arms, localized spur- or crescent-shaped features. The majority of theoretical studies with the aim of interpreting the observational results have focused on disk models with planets, under the assumption that the disk substructures are due to the disk-planet interaction. However, so far only in very few cases exoplanets have been detected in these systems. Furthermore, some substructures are expected to appear \textit{before} planets form, as they are necessary to drive the concentration of small solids which can lead to the formation of planetesimals. In this work we present observational predictions from high-resolution 3D radiative hydrodynamical models which follow the evolution of gas and solids in a protoplanetary disk. We focus on substructures in the distribution of millimeter-sized and smaller solid particles produced by the vertical shear instability. We show that their characteristics are compatible with some of the shallow gaps detected in recent observations at sub-mm/mm wavelengths, and present predictions for future observations with better sensitivity and angular resolution with ALMA and a Next Generation Very Large Array. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 22 pages, 14 figures, accepted for publication in The Astrophysical Journal

Observability of the Vertical Shear Instability in protoplanetary disk CO kinematics

Authors: Marcelo Barraza-Alfaro, Mario Flock, Sebastian Marino, Sebastián Pérez

Abstract: Context. Dynamical and turbulent motions of gas in a protoplanetary disk are crucial for their evolution and affect planet formation. Recent observations suggest weak turbulence in the disk's outer regions. However, the physical mechanism of turbulence in these outer regions remains uncertain. The vertical shear instability (VSI) is a promising mechanism to produce turbulence in disks. Aims. Our a… ▽ More Context. Dynamical and turbulent motions of gas in a protoplanetary disk are crucial for their evolution and affect planet formation. Recent observations suggest weak turbulence in the disk's outer regions. However, the physical mechanism of turbulence in these outer regions remains uncertain. The vertical shear instability (VSI) is a promising mechanism to produce turbulence in disks. Aims. Our aim is to study the observability of the gas velocity structure produced by the VSI via CO kinematics with ALMA. Methods. We perform global 3D hydrodynamical simulations of a VSI-unstable disk. We post-process the simulation results with radiative transfer calculations, and produce synthetic predictions of CO rotational emission lines. Following, we compute the line of sight velocity map, and its deviations from a sub-Keplerian equilibrium solution. We explore the detectability of the VSI by identifying kinematic signatures using realistic simulated observations. Results. Our 3D simulations of the VSI show the steady state dynamics of the gas in great detail. From the velocity structure we infer a turbulent stress value of $α_{rφ}=1.4 \times 10^{-4}$. On large scales, we observe velocity deviations of 50 m s$^{-1}$ as axisymmetric rings. We find optimal conditions at $i \lesssim 20^{\circ}$ to trace for the kinematic structures of the VSI. We found that current diagnostics to constrain gas turbulence from non-thermal broadening of the line emission are not applicable to anisotropic VSI turbulence. Conclusions. The detection of kinematic signatures produced by the VSI is possible with ALMA. Observations including an extended antenna configuration combined with the highest spectral resolution available are needed for robust detection. The characterization of the large-scale velocity perturbations is required to constrain the turbulence level produced by the VSI from gas observations. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A. 17 pages with 10 figures and 1 table

Journal ref: A&A 653, A113 (2021)

Thermal Wave Instability as an Origin of Gap and Ring Structures in Protoplanetary Disks

Authors: Takahiro Ueda, Mario Flock, Tilman Birnstiel

Abstract: Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk… ▽ More Recent millimeter and infrared observations have shown that gap and ring-like structures are common in both dust thermal emission and scattered-light of protoplanetary disks. We investigate the impact of the so-called Thermal Wave Instability (TWI) on the millimeter and infrared scattered-light images of disks. We perform 1+1D simulations of the TWI and confirm that the TWI operates when the disk is optically thick enough for stellar light, i.e., small-grain-to-gas mass ratio of $\gtrsim0.0001$. The mid-plane temperature varies as the waves propagate and hence gap and ring structures can be seen in both millimeter and infrared emission. The millimeter substructures can be observed even if the disk is fully optically thick since it is induced by the temperature variation, while density-induced substructures would disappear in the optically thick regime. The fractional separation between TWI-induced ring and gap is $Δr/r \sim$ 0.2-0.4 at $\sim$ 10-50 au, which is comparable to those found by ALMA. Due to the temperature variation, snow lines of volatile species move radially and multiple snow lines are observed even for a single species. The wave propagation velocity is as fast as $\sim$ 0.6 ${\rm au~yr^{-1}}$, which can be potentially detected with a multi-epoch observation with a time separation of a few years. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 10 pages, 8 figures, accepted for publication in ApJL

Observable scattered light features from inclined and non-inclined planets embedded in protoplanetary discs

Authors: D. Kloster, M. Flock

Abstract: Over the last few years instruments such as VLT/SPHERE and Subaru/HiCIAO have been able to take detailed scattered light images of protoplanetary discs. Many of the features observed in these discs are generally suspected to be caused by an embedded planet, and understanding the cause of these features requires detailed theoretical models. In this work we investigate disc-planet interactions using… ▽ More Over the last few years instruments such as VLT/SPHERE and Subaru/HiCIAO have been able to take detailed scattered light images of protoplanetary discs. Many of the features observed in these discs are generally suspected to be caused by an embedded planet, and understanding the cause of these features requires detailed theoretical models. In this work we investigate disc-planet interactions using the PLUTO code to run 2D and 3D hydrodynamic simulations of protoplanetary discs with embedded 30 M$_{\oplus}$ and 300 M$_{\oplus}$ planets on both an inclined ($i = 2.86^{\circ}$) and non-inclined orbit, using an $α$-viscosity of $4 \times 10^{-3}$. We produce synthetic scattered-light images of these discs at \emph{H-band} wavelengths using the radiative transfer code RADMC3D. We find that while the surface density evolution in 2D and 3D simulations of inclined and non-inclined planets remain fairly similar, their observational appearance is remarkably different. Most of the features seen in the synthetic \emph{H-band} images are connected to density variations of the disc at around 3.3 scale heights above and below the midplane, which emphasizes the need for 3D simulations. Planets on sustained orbital inclinations disrupt the disc's upper-atmosphere and produce radically different observable features and intensity profiles, including shadowing effects and intensity variation in the order of 10-20 times the surrounding background. The vertical optical depth to the disc midplane for \emph{H-band} wavelengths is $τ\approx 20$ in the disc gap created by the high-mass planet. We conclude that direct imaging of planets embedded in the disc remains difficult to observe, even for massive planets in the gap. △ Less

Submitted 24 May, 2021; originally announced May 2021.

Journal ref: MNRAS, 487, 5372 (2019)

Streaming instability in a global patch simulation of protoplanetary disks

Authors: Mario Flock, Andrea Mignone

Abstract: In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two comp… ▽ More In the recent years, sub/mm observations of protoplanetary disks have discovered an incredible diversity of substructures in the dust emission. An important result was the finding that dust grains of mm size are embedded in very thin dusty disks. This implies that the dust mass fraction in the midplane becomes comparable to the gas, increasing the importance of the interaction between the two components there. We address this problem by means of numerical 2.5D simulations in order to study the gas and dust interaction in fully global stratified disks. To this purpose, we employ the recently developed dust grain module in the PLUTO code. Our model focuses on a typical T Tauri disk model, simulating a short patch of the disk at 10 au which includes grains of constant Stokes number of $St=0.01$ and $St=0.1$, corresponding to grains with sizes of 0.9 cm and 0.9 mm, respectively, for the given disk model. By injecting a constant pebble flux at the outer domain, the system reaches a quasi steady state of turbulence and dust concentrations driven by the streaming instability. For our given setup and using resolutions up to 2500 cells per scale height we resolve the streaming instability, leading to local dust clumping and concentrations. Our results show dust density values of around 10-100 times the gas density with a steady state pebble flux between $3.5 \times 10^{-4}$ and $2.5 \times 10^{-3} M_{\rm Earth}/\mathit{year}$ for the models with $\mathit{St}=0.01$ and $\mathit{St}=0.1$. The grain size and pebble flux for model $\mathit{St}=0.01$ compares well with dust evolution models of the first million years of disk evolution. For those grains the scatter opacity dominates the extinction coefficient at mm wavelengths. These types of global dust and gas simulations are a promising tool for studies of the gas and dust evolution at pressure bumps in protoplanetary disks. △ Less

Submitted 28 March, 2021; originally announced March 2021.

Comments: accepted at A&A

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Late infall causing disk misalignment and dynamic structures in SU Aur

Authors: C. Ginski, S. Facchini, J. Huang, M. Benisty, D. Vaendel, L. Stapper, C. Dominik, J. Bae, F. Menard, G. Muro-Arena, M. Hogerheijde, M. McClure, R. G. van Holstein, T. Birnstiel, Y. Boehler, A. Bohn, M. Flock, E. E. Mamajek, C. F. Manara, P. Pinilla, C. Pinte, A. Ribas

Abstract: Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resol… ▽ More Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 18 pages, 12 figures, published in ApJL on 18-02-2021

Journal ref: ApJL, 908, L25 (2021)

Physical and Chemical Structure of the Disk and Envelope of the Class 0/I Protostar L1527

Authors: Lizxandra Flores-Rivera, Susan Terebey, Karen Willacy, Andrea Isella, Neal Turner, Mario Flock

Abstract: Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that… ▽ More Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that allow us to distinguish the disk, the infalling envelope and outflow of this Class 0/I protostar. RadChemT combines dynamics, radiative transfer, gas chemistry and gas-grain reactions to generate models which can be directly compared with observations for individual protostars. Rather than individually fit abundances to a large number of free parameters, we aim to best match the spectral line maps by (i) adopting a physical model based on density structure and luminosity derived primarily from previous work that fit SED and 2D imaging data, updating it to include a narrow jet detected in CARMA and ALMA data near ($\leq 75$au) the protostar, and then (ii) computing the resulting astrochemical abundances for 292 chemical species. Our model reproduces the C$^{18}$O and N$_{2}$H$^{+}$ line strengths within a factor of 3.0; this is encouraging considering the pronounced abundance variation (factor $> 10^3$) between the outflow shell and CO snowline region near the midplane. Further, our modeling confirms suggestions regarding the anti-correlation between N$_{2}$H$^{+}$ and the CO snowline between 400 au to 2,000 au from the central star. Our modeling tools represent a new and powerful capability with which to exploit the richness of spectral line imaging provided by modern submillimeter interferometers. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 25 pages, 9 figures, and 7 tables. Accepted by ApJ

The Ophiuchus DIsc Survey Employing ALMA (ODISEA)-III: the evolution of substructures in massive discs at 3-5 au resolution

Authors: Lucas A. Cieza, Camilo González-Ruilova, Antonio S. Hales, Paola Pinilla, Dary Ruíz-Rodríguez, Alice Zurlo, Simón Casassus, Sebastián Pérez, Hector Cánovas, Carla Arce-Tord, Mario Flock, Nicolas Kurtovic, Sebastian Marino, Pedro H. Nogueira, Laura Perez, Daniel J. Price, David A. Principe, Jonathan P. Williams

Abstract: We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in… ▽ More We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $<$10 au). We find that two discs around embedded protostars lack the clear gaps and rings that are ubiquitous in more evolved sources with Class II SEDs. Our sample includes 5 objects with previously known large dust cavities (r $>$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques. △ Less

Submitted 5 December, 2020; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

Hydrodynamical simulations of protoplanetary disks including irradiation of stellar photons. I. Resolution study for Vertical Shear Instability (VSI)

Authors: Lizxandra Flores-Rivera, Mario Flock, Ryohei Nakatani

Abstract: In recent years hydrodynamical (HD) models have become important to describe the gas kinematics in protoplanetary disks, especially in combination with models of photoevaporation and/or magnetic-driven winds. We focus on diagnosing the the vertical extent of the VSI at 203 cells per scale height and allude at what resolution per scale height we obtain convergence. Finally, we determine the regions… ▽ More In recent years hydrodynamical (HD) models have become important to describe the gas kinematics in protoplanetary disks, especially in combination with models of photoevaporation and/or magnetic-driven winds. We focus on diagnosing the the vertical extent of the VSI at 203 cells per scale height and allude at what resolution per scale height we obtain convergence. Finally, we determine the regions where EUV, FUV and X-Rays are dominant in the disk. We perform global HD simulations using the PLUTO code. We adopt a global isothermal accretion disk setup, 2.5D (2 dimensions, 3 components) which covers a radial domain from 0.5 to 5.0 and an approximately full meridional extension. We determine the 50 cells per scale height to be the lower limit to resolve the VSI. For higher resolutions, greater than 50 cells per scale height, we observe the convergence for the saturation level of the kinetic energy. We are also able to identify the growth of the `body' modes, with higher growth rate for higher resolution. Full energy saturation and a turbulent steady state is reached after 70 local orbits. We determine the location of the EUV-heated region defined by the radial column density to be 10$^{19}$ cm$^{-2}$ located at $H_\mathrm{R}\sim9.7$, and the FUV/X-Rays-heated boundary layer defined by 10$^{22}$ cm$^{-2}$ located at $H_\mathrm{R}\sim6.2$, making it necessary to introduce the need of a hot atmosphere. For the first time, we report the presence of small scale vortices in the r-Z plane, between the characteristic layers of large scale vertical velocity motions. Such vortices could lead to dust concentration, promoting grain growth. Our results highlight the importance to combine photoevaporation processes in the future high-resolution studies of the turbulence and accretion processes in disks. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Journal ref: A&A 644, A50 (2020)

A Tale of Two Transition Disks: ALMA long-baseline observations of ISO-Oph 2 reveal two closely packed non-axisymmetric rings and a $\sim$2 au cavity

Authors: Camilo González-Ruilova, Lucas A. Cieza, Antonio S. Hales, Sebastián Pérez, Alice Zurlo, Carla Arce-Tord, Simón Casassus, Hector Cánovas, Mario Flock, Gregory J. Herczeg, Paola Pinilla, Daniel J. Price, David A. Principe, Dary Ruíz-Rodríguez, Jonathan P. Williams

Abstract: ISO-Oph 2 is a wide-separation (240 au) binary system where the primary star harbors a massive (M$_{dust}$ $\sim$40 M$_{\oplus}$) ring-like disk with a dust cavity $\sim$50 au in radius and the secondary hosts a much lighter (M$_{dust}$ $\sim$0.8 M$_{\oplus}$) disk. As part of the high-resolution follow-up of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA) project, we present 1.3 mm continuum… ▽ More ISO-Oph 2 is a wide-separation (240 au) binary system where the primary star harbors a massive (M$_{dust}$ $\sim$40 M$_{\oplus}$) ring-like disk with a dust cavity $\sim$50 au in radius and the secondary hosts a much lighter (M$_{dust}$ $\sim$0.8 M$_{\oplus}$) disk. As part of the high-resolution follow-up of the "Ophiuchus Disk Survey Employing ALMA" (ODISEA) project, we present 1.3 mm continuum and $^{12}$CO molecular line observations of the system at 0''02 (3 au) resolution. We resolve the disk around the primary into two non-axisymmetric rings and find that the disk around the secondary is only $\sim$7 au across and also has a dust cavity (r $\sim$2.2 au). Based on the infrared flux ratio of the system and the M0 spectral type of the primary, we estimate the mass of the companion to be close to the brown dwarf limit. Hence, we conclude that the ISO-Oph 2 system contains the largest and smallest cavities, the smallest measured disk size, and the resolved cavity around the lowest mass object (M$_{\star}$ $\sim$0.08 M$_\odot$) in Ophiuchus. From the $^{12}$CO data, we find a bridge of gas connecting both disks. While the morphology of the rings around the primary might be due to an unseen disturber within the cavity, we speculate that the bridge might indicate an alternative scenario in which the secondary has recently flown by the primary star causing the azimuthal asymmetries in its disk. The ISO-Oph 2 system is therefore a remarkable laboratory to study disk evolution, planet formation, and companion-disk interactions. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: 10 pages, 4 figures, 1 table

MSC Class: 85-11 ACM Class: J.2

ALMA and VLA Observations of EX Lupi in its Quiescent State

Authors: Jacob Aaron White, Á. Kóspál, A. G. Hughes, P. Ábrahám, V. Akimkin, A. Banzatti, L. Chen, F. Cruz-Sáenz de Miera, A. Dutrey, M. Flock, S. Guilloteau, A. S. Hales, T. Henning, K. Kadam, D. Semenov, A. Sicilia-Aguilar, R. Teague, E. I. Vorobyov

Abstract: Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with ar… ▽ More Extreme outbursts in young stars may be a common stage of pre-main-sequence stellar evolution. These outbursts, caused by enhanced accretion and accompanied by increased luminosity, can also strongly impact the evolution of the circumstellar environment. We present ALMA and VLA observations of EX Lupi, a prototypical outburst system, at 100 GHz, 45 GHz, and 15 GHz. We use these data, along with archival ALMA 232 GHz data, to fit radiative transfer models to EX Lupi's circumstellar disk in its quiescent state following the extreme outburst in 2008. The best fit models show a compact disk with a characteristic dust radius of 45 au and a total mass of 0.01 M$_{\odot}$. Our modeling suggests grain growth to sizes of at least 3 mm in the disk, possibly spurred by the recent outburst, and an ice line that has migrated inward to $0.2-0.3$ au post-outburst. At 15 GHz, we detected significant emission over the expected thermal disk emission which we attribute primarily to stellar (gyro)synchrotron and free-free disk emission. Altogether, these results highlight what may be a common impact of outbursts on the circumstellar dust. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: Accepted to ApJ, 15 pages, 8 figures

High Resolution Parameter Study of the Vertical Shear Instability

Authors: Natascha Manger, Hubert Klahr, Wilhelm Kley, Mario Flock

Abstract: Theoretical models of protoplanetary disks have shown the Vertical Shear Instability (VSI) to be a prime candidate to explain turbulence in the dead zone of the disk. However, simulations of the VSI have yet to show consistent levels of key disk turbulence parameters like the stress-to-pressure ratio $α$. We aim to reconcile these different values by performing a parameter study on the VSI with fo… ▽ More Theoretical models of protoplanetary disks have shown the Vertical Shear Instability (VSI) to be a prime candidate to explain turbulence in the dead zone of the disk. However, simulations of the VSI have yet to show consistent levels of key disk turbulence parameters like the stress-to-pressure ratio $α$. We aim to reconcile these different values by performing a parameter study on the VSI with focus on the disk density gradient $p$ and aspect ratio $h := H/R$. We use full 2$π$ 3D simulations of the disk for chosen set of both parameters. All simulations are evolved for 1000 reference orbits, at a resolution of 18 cells per h. We find that the saturated stress-to-pressure ratio in our simulations is dependent on the disk aspect ratio with a \review{strong} scaling of $α\propto h^{2.6}$, in contrast to the traditional $α$ model, where viscosity scales as $ν\propto αh^2$ with a constant $α$. We also observe consistent formation of large scale vortices across all investigated parameters. The vortices show uniformly aspect ratios of $χ\approx 10$ and radial widths of approximately 1.5 $H$. With our findings we can reconcile the different values reported for the stress-to-pressure ratio from both isothermal and radiation hydrodynamics models, and show long-term evolution effects of the VSI that could aide in the formation of planetesimals. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: 15 pages, 18 figures, submitted to MNRAS

On moving shadows and pressure bumps in HD 169142

Authors: Gesa H. -M. Bertrang, Mario Flock, Miriam Keppler, Trifon Trifonov, Anna B. T. Penzlin, Henning Avenhaus, Thomas Henning, Matias Montesinos

Abstract: The search for young planets had its first breakthrough with the detection of the accreting planet PDS70b. In this study, we aim to broaden our understanding towards the formation of multi-planet systems such as HR8799 or the Solar System. Our previous study on HD169142, one of the closest Herbig stars, points towards a shadow-casting protoplanetary candidate. Here, we present follow-up observatio… ▽ More The search for young planets had its first breakthrough with the detection of the accreting planet PDS70b. In this study, we aim to broaden our understanding towards the formation of multi-planet systems such as HR8799 or the Solar System. Our previous study on HD169142, one of the closest Herbig stars, points towards a shadow-casting protoplanetary candidate. Here, we present follow-up observations to test our previously proposed hypothesis. We set our new data into context with previous observations to follow structural changes in the disk over the course of 6 years. We find spatially resolved systematic changes in the position of the previously described surface brightness dip in the inner ring. We further find changes in the brightness structure in azimuthal direction along the ring. And finally, a comparison of our SPHERE data with recent ALMA observations reveals a wavelength dependent radial profile of the bright ring. The time-scale on which the changes in the ring's surface brightness occur suggest that they are caused by a shadow cast by a 1-10Mj planet surrounded by dust, an orbit comparable to those of the giant planets in our own Solar System. Additionally, we find the first indications for temperature-induced instabilities in the ring. And finally, we trace a pressure maxima, for the first time spatially resolved, with a width of 4.5au. The density distribution of the ring at mm wavelengths around the pressure maxima could further indicate effects from snow lines or even the dynamics and feedback of the larger grains. △ Less

Submitted 29 April, 2021; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: Submitted to ApJ, in rev

Gas and dust dynamics in starlight-heated protoplanetary disks

Authors: Mario Flock, Neal J. Turner, Richard P. Nelson, Wladimir Lyra, Natascha Manger, Hubert Klahr

Abstract: Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI… ▽ More Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$μ$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $α\simeq 10^{-4}$. The value of $α$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $α$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: ApJ accepted

A two-moment radiation hydrodynamics scheme applicable to simulations of planet formation in circumstellar disks

Authors: Julio David Melon Fuksman, Hubert Klahr, Mario Flock, Andrea Mignone

Abstract: We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly… ▽ More We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly anisotropic radiation transport as can be expected in the vicinity of an accreting planet in a protoplanetary disk. To reduce the computational cost caused by the timescale disparity between radiation and matter fields, we integrate their evolution equations separately in an operator-split way, using substepping to evolve the radiation equations. We further increase the code's efficiency by adopting the reduced speed of light approximation (RSLA). Our integration scheme for the evolution equations of radiation fields relies on implicit-explicit schemes, in which radiation-matter interaction terms are integrated implicitly while fluxes are integrated via Godunov-type solvers. The module is suitable for general astrophysical computations in 1, 2, and 3 dimensions in Cartesian, spherical and cylindrical coordinates, and can be implemented on rotating frames. We demonstrate the algorithm performance on different numerical benchmarks, paying particular attention to the applicability of the RSLA for computations of physical processes in protoplanetary disks. We show 2D simulations of vertical convection in disks and 3D simulations of gas accretion by planetary cores, which are the first of their kind to be solved with a two-moment approach. △ Less

Submitted 18 January, 2021; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 27 pages, 16 figures

Journal ref: The Astrophysical Journal, 906:78, 2021

Birth of convective low-mass to high-mass second Larson cores

Authors: Asmita Bhandare, Rolf Kuiper, Thomas Henning, Christian Fendt, Mario Flock, Gabriel-Dominique Marleau

Abstract: Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the depen… ▽ More Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 22 pages, 20 figures, and 2 tables. Accepted for publication in A&A. Collapse movie: https://keeper.mpdl.mpg.de/f/f04abdeabdf3472fb56d/

Journal ref: A&A 638, A86 (2020)

Population synthesis of exocometary gas around A stars

Authors: S. Marino, M. Flock, Th. Henning, Q. Kral, L. Matrà, M. C. Wyatt

Abstract: The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon prod… ▽ More The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long timescales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded onto an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity ($α\sim0.1$), indicating that gas is lost on timescales $\sim1-10$ Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesise that the observed carbon cavities could be due to radiation pressure or accreting planets. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS, 23 pages, 13 figures

Planet formation and migration near the silicate sublimation front in protoplanetary disks

Authors: Mario Flock, Neal J. Turner, Gijs D. Mulders, Yasuhiro Hasegawa, Richard P. Nelson, Bertram Bitsch

Abstract: The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet for… ▽ More The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems. △ Less

Submitted 9 October, 2019; originally announced October 2019.

Comments: A&A in press, MPIA press release 10.10.2019

A Particle Module for the PLUTO Code: III -- Dust

Authors: A. Mignone, M. Flock, B. Vaidya

Abstract: The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust p… ▽ More The implementation of a new particle module describing the physics of dust grains coupled to the gas via drag forces is the subject of this work. The proposed particle-gas hybrid scheme has been designed to work in Cartesian as well as in cylindrical and spherical geometries. The numerical method relies on a Godunov-type second-order scheme for the fluid and an exponential midpoint rule for dust particles which overcomes the stiffness introduced by the linear coupling term. Besides being time-reversible and globally second-order accurate in time, the exponential integrator provides energy errors which are always bounded and it remains stable in the limit of arbitrarily small particle stopping times yielding the correct asymptotic solution. Such properties make this method preferable to the more widely used semi-implicit or fully implicit schemes at a very modest increase in computational cost. Coupling between particles and grid quantities is achieved through particle deposition and field-weighting techniques borrowed from Particle-In-Cell simulation methods. In this respect, we derive new weight factors in curvilinear coordinates that are more accurate than traditional volume- or area-weighting. A comprehensive suite of numerical benchmarks is presented to assess the accuracy and robustness of the algorithm in Cartesian, cylindrical and spherical coordinates. Particular attention is devoted to the streaming instability which is analyzed in both local and global disk models. The module is part of the PLUTO code for astrophysical gas-dynamics and it is mainly intended for the numerical modeling of protoplanetary disks in which solid and gas interact via aerodynamic drag. △ Less

Submitted 28 August, 2019; originally announced August 2019.

Comments: 22 pages, 13 figures

The Radial Distribution of Dust Particles in the HL Tau Disk from ALMA and VLA Observations

Authors: Carlos Carrasco-González, Anibal Sierra, Mario Flock, Zhaohuan Zhu, Thomas Henning, Claire Chandler, Roberto Galván-Madrid, Enrique Macías, Guillem Anglada, Hendrik Linz, Mayra Osorio, Luis F. Rodríguez, Leonardo Testi, José M. Torrelles, Laura Pérez, Yao Liu

Abstract: Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these… ▽ More Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules. △ Less

Submitted 24 September, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Characterization of Ring Substructures in the Protoplanetary Disk of HD 169142 from Multi-Wavelength ALMA Observations

Authors: Enrique Macias, Catherine Espaillat, Mayra Osorio, Guillem Anglada, Jose M. Torrelles, Carlos Carrasco-Gonzalez, Mario Flock, Hendrik Linz, Gesa H. M. Bertrang, Thomas Henning, Jose F. Gomez, Nuria Calvet, William R. F. Dent

Abstract: We present a detailed multi-wavelength characterization of the multi-ring disk of HD 169142. We report new ALMA observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, te… ▽ More We present a detailed multi-wavelength characterization of the multi-ring disk of HD 169142. We report new ALMA observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is $\sim1$ cm, with slightly flatter power-law size distributions than the ISM-like size distribution ($p\sim3.5$) found in the gaps. In particular, the inner ring ($\sim26$ au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet-disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multi-wavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes. △ Less

Submitted 16 July, 2019; originally announced July 2019.

Comments: 19 pages, 8 figures, accepted for publication in ApJ

Intrinsic polarisation of elongated porous dust grains

Authors: Florian Kirchschlager, Gesa H. -M. Bertrang, Mario Flock

Abstract: ALMA observations revealed recently polarised radiation of several protoplanetary disks in the (sub-)millimetre wavelength range. Besides self-scattering of large particles, thermal emission by elongated grains is a potential source for the detected polarisation signal. We calculate the wavelength dependent absorption and intrinsic polarisation of spheroidally shaped, micrometre and sub-millimetre… ▽ More ALMA observations revealed recently polarised radiation of several protoplanetary disks in the (sub-)millimetre wavelength range. Besides self-scattering of large particles, thermal emission by elongated grains is a potential source for the detected polarisation signal. We calculate the wavelength dependent absorption and intrinsic polarisation of spheroidally shaped, micrometre and sub-millimetre sized dust grains using the discrete dipole approximation. In particular, we analyse the impact of dust grain porosity which appears to be present in disks when small grains coagulate to form larger aggregates. For the first time our results show that (a) the intrinsic polarisation decreases for increasing grain porosity and (b) the polarisation orientation flips by 90 degree for certain ratios of wavelength to grain size. We present a new method to constrain grain porosity and the grain size in protoplanetary disks using multi-wavelength polarisation observations in the far-infrared to millimetre wavelengths. Finally, we find that moderate grain porosities ($\mathcal{P}\lesssim0.7$) potentially explain the observed polarisation fraction in the system HD 142527 while highly porous grains ($\mathcal{P}>0.7$) fail unless the grain's axis ratio is extraordinarily large. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: 10 pages, 10 figures

An analytical model of radial dust trapping in protoplanetary disks

Authors: Anibal Sierra, Susana Lizano, Enrique Macías, Carlos Carrasco-González, Mayra Osorio, Mario Flock

Abstract: We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytic… ▽ More We study dust concentration in axisymmetric gas rings in protoplanetary disks. Given the gas surface density, we derived an analytical total dust surface density by taking into account the differential concentration of all the grain sizes. This model allows us to predict the local dust-to-gas mass ratio and the slope of the particle size distribution, as a function of radius. We test this analytical model comparing it with a 3D magneto-hydrodynamical simulation of dust evolution in an accretion disk. The model is also applied to the disk around HD 169142. By fitting the disk continuum observations simultaneously at $λ= 0.87$, 1.3, 3.0 mm, we obtain a global dust-to-gas mass ratio $ε_{\rm global} = 1.05 \times 10^{-2}$ and a viscosity coefficient $α= 1.35 \times 10^{-2}$. This model can be easily implemented in numerical simulations of accretion disks. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Highly structured disk around the planet host PDS 70 revealed by high-angular resolution observations with ALMA

Authors: M. Keppler, R. Teague, J. Bae, M. Benisty, T. Henning, R. van Boekel, E. Chapillon, P. Pinilla, J. P. Williams, G. H. -M. Bertrang, S. Facchini, M. Flock, Ch. Ginski, A. Juhasz, H. Klahr, Y. Liu, A. Mueller, L. M. Perez, A. Pohl, G. Rosotti, M. Samland, D. Semenov

Abstract: PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the… ▽ More PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology. △ Less

Submitted 17 April, 2019; v1 submitted 20 February, 2019; originally announced February 2019.

Comments: Accepted for publication in A&A

Journal ref: A&A 625, A118 (2019)

Spatial segregation of dust grains in transition disks. SPHERE observations of 2MASS J16083070-3828268 and RXJ1852.3-3700

Authors: M. Villenave, M. Benisty, W. R. F. Dent, F. Menard, A. Garufi, C. Ginski, P. Pinilla, C. Pinte, J. P. Williams, J. de Boer, J. -I. Morino, M. Fukagawa, C. Dominik, M. Flock, T. Henning, A. Juhasz, M. Keppler, G. Muro-Arena, J. Olofsson, L. M. Perez, G. van der Plas, A. Zurlo, M. Carle, P. Feautrier, A. Pavlov , et al. (5 additional authors not shown)

Abstract: Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered li… ▽ More Context. The mechanisms governing the opening of cavities in transition disks are not fully understood. Several processes have been proposed but their occurrence rate is still unknown. Aims. We present spatially resolved observations of two transition disks and aim at constraining their vertical and radial structure using multiwavelength observations. Methods. We have obtained near-IR scattered light observations with VLT/SPHERE of the transition disks J1608 and J1852. We complement our datasets with ALMA observations and with unresolved photometric observations covering a wide range of wavelengths. We performed radiative transfer modeling to analyze the morphology of the disks and compare the results with a sample of 20 other transition disks observed with both SPHERE and ALMA. Results. The scattered light image of J1608 reveals a very inclined disk, with two bright lobes and a large cavity. J1852 shows an inner ring extending beyond the coronagraphic radius up to 15au, a gap and a second ring at 42au. Our radiative transfer model of J1608 indicates that the millimeter-sized grains are less extended vertically and radially than the micron-sized grains, indicating advanced settling and radial drift. We find good agreement with the observations of J1852 with a similar model, but due to the low inclination of the system, the model remains partly degenerate. The analysis of 22 transition disks shows that, in general, the cavities observed in scattered light are smaller than the ones detected at millimeter wavelengths. Conclusions. The analysis of a sample of transition disks indicates that the small grains can flow inward of the region where millimeter grains are trapped. While 15 out of the 22 cavities in our sample could be explained by a planet of less than 13 Jupiter masses, the others either require the presence of a more massive companion or of several low-mass planets. △ Less

Submitted 12 February, 2019; originally announced February 2019.

Comments: Accepted in Astronomy and Astrophysics

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

Imaging the Distribution of Solids in Planet-forming Disks undergoing Hydrodynamical Instabilities with the Next Generation Very Large Array

Authors: L. Ricci, M. Flock, D. Blanco, W. Lyra

Abstract: We present simulations of the capabilities of the Next Generation Very Large Array to image at high angular resolution substructures in the dust emission of protoplanetary disks. The main goal of this study is to investigate the kinds of substructures that are expected by state-of-the-art 3D simulations of disks and that an instrument like the ngVLA, with its current design, can detect. The disk s… ▽ More We present simulations of the capabilities of the Next Generation Very Large Array to image at high angular resolution substructures in the dust emission of protoplanetary disks. The main goal of this study is to investigate the kinds of substructures that are expected by state-of-the-art 3D simulations of disks and that an instrument like the ngVLA, with its current design, can detect. The disk simulations adopted in this investigation consist of global 3D radiation-hydrodynamics models with embedded particles, the latter representing dust grains. Our work shows that the ngVLA can detect and spatially resolve, down to sub-astronomical unit scales in disks in nearby star forming regions, the dust continuum emission at 3mm from azimuthal asymmetric structures, as well as from weak rings and gaps produced in these models as a consequence of the vertical shear instability (VSI). This hydrodynamical instability has been proposed to generate turbulence in regions of weak coupling between the disk gas and magnetic field, as well as to form vortices which may be preferred locations of planetesimal formation. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 8 pages, 2 figures, to appear in the ngVLA memo series

Dust-Pileup at the Dead-Zone Inner Edge and Implications for the Disk Shadow

Authors: Takahiro Ueda, Mario Flock, Satoshi Okuzumi

Abstract: We perform simulations of the dust and gas disk evolution to investigate the observational features of a dust-pileup at the dead-zone inner edge. We show that the total mass of accumulated dust particles is sensitive to the turbulence strength in the dead zone, $α_{\rm dead}$, because of the combined effect of turbulence-induced particle fragmentation (which suppresses particle radial drift) and t… ▽ More We perform simulations of the dust and gas disk evolution to investigate the observational features of a dust-pileup at the dead-zone inner edge. We show that the total mass of accumulated dust particles is sensitive to the turbulence strength in the dead zone, $α_{\rm dead}$, because of the combined effect of turbulence-induced particle fragmentation (which suppresses particle radial drift) and turbulent diffusion. For a typical critical fragmentation velocity of silicate dust particles of $1~{\rm m~s^{-1}}$, the stress to pressure ratio $α_{\rm dead}$ needs to be lower than $3 \times 10^{-4}$ for dust trapping to operate. The obtained dust distribution is postprocessed using the radiative transfer code RADMC-3D to simulate infrared scattered-light images of the inner part of protoplanetary disks with a dust pileup. We find that a dust pileup at the dead-zone inner edge, if present, casts a shadow extending out to $\sim 10~{\rm au}$. In the shadowed region the temperature significantly drops, which in some cases yields even multiple water snow lines. We also find that even without a dust pileup at the dead-zone inner edge, the disk surface can become thermally unstable, and the excited waves can naturally produce shadows and ring-like structures in observed images. This mechanism might account for the ring-like structures seen in the scattered-light images of some disks, such as the TW Hya disk. △ Less

Submitted 23 November, 2018; originally announced November 2018.

Comments: 16 pages, 13 figures, accepted for publication in ApJ

Kinematic links and the co-evolution of MHD winds, jets, and inner disks from a high-resolution optical [OI] survey

Authors: Andrea Banzatti, Ilaria Pascucci, Suzan Edwards, Min Fang, Uma Gorti, Mario Flock

Abstract: We present a survey of optical [OI] emission at 6300 Å toward 65 TTauri stars at the spectral resolution of $\sim7$ km/s. Past work identified a highly blueshifted velocity component (HVC) tracing microjets, and a less blueshifted low-velocity-component (LVC) attributed to winds. We focus here on the LVC kinematics to investigate links between winds, jets, accretion, and disk dispersal. We track t… ▽ More We present a survey of optical [OI] emission at 6300 Å toward 65 TTauri stars at the spectral resolution of $\sim7$ km/s. Past work identified a highly blueshifted velocity component (HVC) tracing microjets, and a less blueshifted low-velocity-component (LVC) attributed to winds. We focus here on the LVC kinematics to investigate links between winds, jets, accretion, and disk dispersal. We track the behavior of four types of LVC components: a broad and narrow component ("BC" and "NC") in LVCs that are decomposed into two Gaussians, which typically have an HVC, and the single-Gaussian LVC profiles separated into those that have an HVC ("SCJ") and those that do not ("SC"). LVC centroid velocities and line widths correlate with HVC equivalent width and accretion luminosity, suggesting that LVC/winds and HVC/jets are kinematically linked and connected to accretion. The deprojected HVC velocity correlates with accretion luminosity, showing that faster jets come with higher accretion. BC and NC kinematics correlate and their blueshifts are maximum at $\sim 35^{\circ}$, suggesting a conical wind geometry with this semi-opening angle. Only SCs include $n_{13-31}$ up to $\sim3$ and their properties correlate with this infrared index, showing that [OI] emission recedes to larger radii as the inner dust is depleted, tracing less dense/hot gas and a decrease in wind velocity. All together, these findings support a scenario where optically thick, accreting inner disks launch radially-extended MHD disk winds that feed jets, and where inner disk winds recede to larger radii and jets disappear in concert with dust depletion. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: Accepted for publication on ApJ

A New Look at T Tauri Star Forbidden Lines: MHD Driven Winds from the Inner Disk

Authors: Min Fang, Ilaria Pascucci, Suzan Edwards, Uma Gorti, Andrea Banzatti, Mario Flock, Patrick Hartigan, Gregory J. Herczeg, Andrea K. Dupree

Abstract: Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($Δv\sim$6\kms) analysis of [S II] $λ$4068, [O I] $λ$5577, and [O I] $λ$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity c… ▽ More Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($Δv\sim$6\kms) analysis of [S II] $λ$4068, [O I] $λ$5577, and [O I] $λ$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity component (HVC) associated with jets, and a low-velocity narrow (LVC-NC) and broad (LVC-BC) components. We confirm previous findings that many LVCs are blueshifted by more than 1.5 kms$^{-1}$ thus most likely trace a slow disk wind. We further show that the profiles of individual components are similar in the three lines. We find that most LVC-BC and NC line ratios are explained by thermally excited gas with temperatures between 5,000$-$10,000 K and electron densities $\sim10^{7}-10^{8}$ cm$^{-3}$. The HVC ratios are better reproduced by shock models with a pre-shock H number density of $\sim10^{6}-10^{7}$ cm$^{-3}$. Using these physical properties, we estimate $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ for the LVC and $\dot{M}_{\rm jet}/\dot{M}_{\rm acc}$ for the HVC. In agreement with previous work, the mass carried out in jets is modest compared to the accretion rate. With the likely assumption that the NC wind height is larger than the BC, the LVC-BC $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ is found to be higher than the LVC-NC. These results suggest that most of the mass loss occurs close to the central star, within a few au, through an MHD driven wind. Depending on the wind height, MHD winds might play a major role in the evolution of the disk mass. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 45 pages, 23 figures, and 7 tables, accepted by ApJ