Showing 1–50 of 123 results for author: Rosotti, G

Seeing the unseen: a method to detect unresolved rings in protoplanetary disks

Authors: Chiara E. Scardoni, Richard A. Booth, Cathie J. Clarke, Giovanni P. Rosotti, Alvaro Ribas

Abstract: While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in pr… ▽ More While high resolution ALMA observations reveal a wealth of substructure in protoplanetary discs, they remain incapable of resolving the types of small scale dust structures predicted, for example, by numerical simulations of the streaming instability. In this Letter, we propose a method to find evidence for unresolved, optically thick dusty rings in protoplanetary disks. We demonstrate that, in presence of unresolved rings, the brightness of an inclined disc exhibits a distinctive emission peak at the minor axis. Furthermore, the azimuthal brightness depends on both the geometry of the rings and the dust optical properties; we can therefore use the azimuthal brightness variations to both detect unresolved rings and probe their properties. By analyzing the azimuthal brightness in the test-case of ring-like substructures formed by streaming instability, we show that the resulting peak is likely detectable by ALMA for typical disc parameters. Moreover, we present an analytic model that not only qualitatively but also quantitatively reproduces the peak found in the simulations, validating its applicability to infer the presence of unresolved rings in observations and characterize their optical properties and shape. This will contribute to the identification of disk regions where streaming instability (and thus planet formation) is occurring. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted by ApJ, 10 pages, 6 figures

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Planet-driven spirals in protoplanetary discs: limitations of the semi-analytical theory for observations

Authors: D. Fasano, A. J. Winter, M. Benisty, G. Rosotti, A. Ruzza, G. Lodato, C. Toci, T. Hilder, A. Izquierdo, D. Price

Abstract: Detecting protoplanets during their formation stage is an important but elusive goal of modern astronomy. Kinematic detections via the spiral wakes in the gaseous disc are a promising avenue to achieve this goal. We aim to test the applicability to observations in the low and intermediate planet mass regimes of a commonly used semi-analytical model for planet induced spiral waves. In contrast with… ▽ More Detecting protoplanets during their formation stage is an important but elusive goal of modern astronomy. Kinematic detections via the spiral wakes in the gaseous disc are a promising avenue to achieve this goal. We aim to test the applicability to observations in the low and intermediate planet mass regimes of a commonly used semi-analytical model for planet induced spiral waves. In contrast with previous works which proposed to use the semi-analytical model to interpret observations, in this study we analyse for the first time both the structure of the velocity and density perturbations. We run a set of FARGO3D hydrodynamic simulations and compare them with the output of the semi-analytic model in the code wakeflow, which is obtained by solving Burgers' equation using the simulations as an initial condition. We find that the velocity field derived from the analytic theory is discontinuous at the interface between the linear and nonlinear regions. After 0.2 r$_p$ from the planet, the behaviour of the velocity field closely follows that of the density perturbations. In the low mass limit, the analytical model is in qualitative agreement with the simulations, although it underestimates the azimuthal width and the amplitude of the perturbations, predicting a stronger decay but a slower azimuthal advance of the shock fronts. In the intermediate regime, the discrepancy increases, resulting in a different pitch angle between the spirals of the simulations and the analytic model. The implementation of a fitting procedure based on the minimisation of intensity residuals is bound to fail due to the deviation in pitch angle between the analytic model and the simulations. In order to apply this model to observations, it needs to be revisited accounting also for higher planet masses. △ Less

Submitted 19 June, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: 10 pages, 7 figures. Accepted for publication in A&A May 27, 2024

Population Synthesis Models Indicate a Need for Early and Ubiquitous Disk Substructures

Authors: Luca Delussu, Tilman Birnstiel, Anna Miotello, Paola Pinilla, Giovanni Rosotti, Sean M. Andrews

Abstract: Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust… ▽ More Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslantα\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Correlations among complex organic molecules around protostars: Effects of physical structure

Authors: P. Nazari, B. Tabone, G. P. Rosotti, E. F. van Dishoeck

Abstract: Column density ratios of complex organic molecules are generally constant across protostellar systems with some low-level scatter. However, the scatter in formamide (NH$_2$CHO) to methanol (CH$_3$OH) column density ratio is one of the highest. This larger scatter is sometimes interpreted as evidence of gas-phase formation of NH$_2$CHO. In this work we propose an alternative interpretation in which… ▽ More Column density ratios of complex organic molecules are generally constant across protostellar systems with some low-level scatter. However, the scatter in formamide (NH$_2$CHO) to methanol (CH$_3$OH) column density ratio is one of the highest. This larger scatter is sometimes interpreted as evidence of gas-phase formation of NH$_2$CHO. In this work we propose an alternative interpretation in which this scatter is produced by differences in the snowline locations related to differences in binding energies of these species and the small-scale structure of the envelope and the disk system. We also include CH$_3$CN in our work as a control molecule which has a similar binding energy to CH$_3$OH. We use radiative transfer models to calculate the emission from these species in protostellar systems with and without disks. The abundances of these species are parameterized in our models. Then we fit the calculated emission lines to find the column densities as done in real observations. We find a correction factor of ~10 to be multiplied by gas-phase $N_{NH_2CHO}/N_{CH_3OH}$ to give the true abundance ratio of these two species in the ices. We find that models with different physical parameters produce a scatter in $N_{NH_2CHO}/N_{CH_3OH}$, comparable with that of observations. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ is larger than that of $N_{CH_3CN}/N_{CH_3OH}$ in models consistent with the observations. We show that the scatter in $N_{NH_2CHO}/N_{CH_3OH}$ will be lower if we correct for the difference in sublimation temperatures of these two species in observations of ~40 protostellar systems with ALMA. The scatter in $N_{NH_2CHO}/N_{CH_3OH}$ can be partially explained by the difference in their binding energies. We conclude that gas-phase chemistry routes for NH$_2$CHO are not necessary to explain the larger scatter of $N_{NH_2CHO}/N_{CH_3OH}$ compared with other ratios. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A

Kaleidoscope of irradiated disks: MUSE observations of proplyds in the Orion Nebula Cluster. I. Sample presentation and ionization front sizes

Authors: Mari-Liis Aru, Karina Mauco, Carlo F. Manara, Thomas J. Haworth, Stefano Facchini, Anna F. McLeod, Anna Miotello, Monika G. Petr-Gotzens, Massimo Robberto, Giovanni P. Rosotti, Silvia Vicente, Andrew Winter, Megan Ansdell

Abstract: In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology… ▽ More In the Orion Nebula Cluster (ONC), protoplanetary disks exhibit ionized gas clouds in the form of a striking teardrop shape as massive stars irradiate the disk material. We present the first spatially and spectrally resolved observations of 12 proplyds, using Integral Field Spectroscopy observations performed with the MUSE instrument in Narrow Field Mode (NFM) on the VLT. We present the morphology of the proplyds in seven emission lines and measure the radius of the ionization front (I-front) of the targets in four tracers, covering transitions of different ionization states for the same element. We also derive stellar masses for the targets. The measurements follow a consistent trend of increasing I-front radius for a decreasing strength of the far-UV radiation as expected from photoevaporation models. By analyzing the ratios of the I-front radii as measured in the emission lines of Ha, [OI] 6300, [OII] 7330, and [OIII] 5007, we observe the ionization stratification, that is, the most ionized part of the flow being the furthest from the disk (and closest to the UV source). The ratios of I-front radii scale in the same way for all proplyds in our sample regardless of the incident radiation. We show that the stratification can help constrain the densities near the I-front by using a 1D photoionization model. We derive the upper limits of photoevaporative mass-loss rates by assuming ionization equilibrium, and estimate values decreasing towards lower impinging radiation. We do not find a correlation between Mloss and stellar mass. The highest mass-loss rate is for the proplyd 244-440. These values of Mloss, combined with estimates of the disk mass with ALMA, confirm previous estimates of the short lifetime of these proplyds. This work demonstrates the potential of this dataset and offers a new set of observables to be used to test current and future models of external photoevaporation. △ Less

Submitted 9 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Submitted to A&A on December 18, 2023. Accepted on April 2, 2024

Journal ref: A&A 687, A93 (2024)

DBNets: A publicly available deep learning tool to measure the masses of young planets in dusty protoplanetary discs

Authors: Alessandro Ruzza, Giuseppe Lodato, Giovanni Pietro Rosotti

Abstract: Current methods to characterize embedded planets in protoplanetary disc observations are severely limited either in their ability to fully account for the observed complex physics or in their computational and time costs. To address this shortcoming, we developed DBNets: a deep learning tool, based on convolutional neural networks, that analyses substructures observed in the dust continuum emissio… ▽ More Current methods to characterize embedded planets in protoplanetary disc observations are severely limited either in their ability to fully account for the observed complex physics or in their computational and time costs. To address this shortcoming, we developed DBNets: a deep learning tool, based on convolutional neural networks, that analyses substructures observed in the dust continuum emission of protoplanetary discs to quickly infer the mass of allegedly embedded planets. We focussed on developing a method to reliably quantify not only the planet mass, but also the associated uncertainty introduced by our modelling and adopted techniques. Our tests gave promising results achieving an 87% reduction of the log Mp mean squared error with respect to an analytical formula fitted on the same data (DBNets metrics: lmse 0.016, r2-score 97%). With the goal of providing the final user of DBNets with all the tools needed to interpret their measurements and decide on their significance, we extensively tested our tool on out-of-distribution data. We found that DBNets can identify inputs strongly outside its training scope returning an uncertainty above a specific threshold and we thus provided a rejection criterion that helps determine the significance of the results obtained. Additionally, we outlined some limitations of our tool: it can be reliably applied only on discs observed with inclinations below approximately 60°, in the optically thin regime, with a resolution 8 times better than the gap radial location and with a signal-to-noise ratio higher than approximately ten. Finally, we applied DBNets to 33 actual observations of protoplanetary discs measuring the mass of 48 proposed planets and comparing our results with the available literature. We confirmed that most of the observed gaps imply planets in the sub-Jupiter regime. DBNets is publicly available at dbnets.fisica.unimi.it. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Rotation curves in protoplanetary disks with thermal stratification

Authors: Paola Martire, Cristiano Longarini, Giuseppe Lodato, Giovanni P. Rosotti, Andrew Winter, Stefano Facchini, Caitlyn Hardiman, Myriam Benisty, Jochen Stadler, Andrés F. Izquierdo, Leonardo Testi

Abstract: In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal strat… ▽ More In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal stratification on the azimuthal velocity of protoplanetary disks. High resolution gas observations are showing velocity differences between CO isotopologues, which cannot be adequately explained with vertically isothermal models. The aim of this work is to determine whether a stratified model can explain this discrepancy. We analytically solve the hydrostatic equilibrium for a stratified disk and we derive the azimuthal velocity. We test the model with SPH numerical simulations and then we use it to fit for star mass, disk mass and scale radius of the sources in the MAPS sample. In particular, we use 12CO and 13CO datacubes. △ Less

Submitted 11 March, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, accepted for publication on A&A

Journal ref: A&A 686, A9 (2024)

Accuracy of ALMA estimates of young disk radii and masses. Predicted observations from numerical simulations

Authors: Ngo-Duy Tung, Leonardo Testi, Ugo Lebreuilly, Patrick Hennebelle, Anaëlle Maury, Ralf S. Klessen, Luca Cacciapuoti, Matthias González, Giovanni Rosotti, Sergio Molinari

Abstract: Protoplanetary disks, which are the natural consequence of the gravitational collapse of the dense molecular cloud cores, host the formation of the known planetary systems in our universe. Substantial efforts have been dedicated to investigating the properties of these disks in the more mature Class II stage, either via numerical simulations of disk evolution from a limited range of initial condit… ▽ More Protoplanetary disks, which are the natural consequence of the gravitational collapse of the dense molecular cloud cores, host the formation of the known planetary systems in our universe. Substantial efforts have been dedicated to investigating the properties of these disks in the more mature Class II stage, either via numerical simulations of disk evolution from a limited range of initial conditions or observations of their dust continuum and line emission from specific molecular tracers. The results coming from these two standpoints have been used to draw comparisons. However, few studies have investigated the main limitations at work when measuring the embedded Class 0/I disk properties from observations, especially in a statistical fashion. In this study, we provide a first attempt to compare the accuracy of some critical disk parameters in Class 0/I systems, as derived on real ALMA observational data, with the corresponding physical parameters that can be directly defined by theoreticians and modellers in numerical simulations. The approach we follow here is to provide full post-processing of the numerical simulations and apply it to the synthetic observations the same techniques used by observers to derive the physical parameters. We performed 3D Monte Carlo radiative transfer and mock interferometric observations of the disk populations formed in a magnetohydrodynamic (MHD) simulation model of disk formation through the collapse of massive clumps with the tools \textsc{Radmc-3d} and \textsc{Casa}, respectively, to obtain their synthetic observations. With these observations, we re-employed the techniques commonly used in disk modelling from their continuum emissions to infer the properties that would most likely be obtained with real interferometers. We then demonstrated how these properties may vary with respect to the gas kinematics analyses and dust continuum modelling. △ Less

Submitted 29 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics, 32 pages, 28 figures

Observing planetesimal formation under streaming instability in the rings of HD 163296

Authors: Francesco Zagaria, Cathie J. Clarke, Richard A. Booth, Stefano Facchini, Giovanni P. Rosotti

Abstract: We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a… ▽ More We introduce a new technique to determine the gas turbulence and surface density in bright disc rings, under the assumption that dust growth is limited by turbulent fragmentation at the ring centre. We benchmark this prescription in HD 163296, showing that our measurements are consistent with available turbulence upper limits and agree with independent estimates of the gas surface density within a factor of two. We combine our results with literature measurements of the dust surface density and grain size to determine the dust-to-gas ratio and Stokes number in the 67 au and 100 au rings. Our estimates suggest that particle clumping is taking place under the effect of streaming instability (SI) in the 100 au ring. Even though in the presence of external isotropic turbulence this process might be hindered, we provide evidence that turbulence is non-isotropic in both rings and likely originating from mechanisms (such as ambipolar diffusion) that could ease particle clumping under SI. Finally, we determine the mass accretion rate under the assumption that the disc is in steady state and turbulence regulates angular momentum transport. Our results are in tension with spectroscopic measurements and suggest that other mechanisms might be responsible for accretion, in qualitative agreement with the detection of a magneto-centrifugal wind in this system. Applying our method to larger samples can be used to statistically assess if SI is a viable mechanism to form planetesimals in bright rings. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 13 pages, 4 figures; accepted for publication on ApJL

Synthetic populations of protoplanetary disks. Impact of magnetic fields and radiative transfer

Authors: U. Lebreuilly, P. Hennebelle, T. Colman, A. Maury, N. -D. Tung, L. Testi, R. Klessen, S. Molinari, B. Commerçon, M. González, E. Pacetti, A. Somigliana, G. Rosotti

Abstract: Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to… ▽ More Protostellar disks are the product of angular momentum conservation during the protostellar collapse. Understanding their formation is crucial because they are the birthplace of planets and because their formation is tightly related to star formation. Unfortunately, the initial properties of Class 0 disks and their evolution are still poorly constrained observationally and theoretically. We aim to better understand the mechanisms that set the statistics of disk properties as well as to study their formation in massive protostellar clumps. We also want to provide the community with synthetic disk populations to better interpret young disk observations. We use the ramses code to model star and disk formation in massive protostellar clumps with MHD including the effect of ambipolar diffusion and RT including the stellar radiative feedback. Those simulations, resolved up to the astronomical unit scale, allow to investigate the formation of disk populations. Magnetic fields play a crucial role in disk formation. A weaker initial field leads to larger and massive disks and weakens the stellar radiative feedback by increasing fragmentation. We find that ambipolar diffusion impacts disk and star formation and leads to very different disk magnetic properties. The stellar radiative feedback also have a strong influence, increasing the temperature and reducing fragmentation. Comparing our disk populations with observations reveals that our models with a mass-to-flux ratio of 10 seems to better reproduce observed disk sizes. This also sheds light on a tension between models and observations for the disk masses. The clump properties and physical modeling impact disk populations significantly. The tension between observations and models for disk mass estimates is critical to solve with synthetic observations in future years, in particular for our comprehension of planet formation. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 20 pages, 15 figures, accepted for publication in Astronomy & Astrophysics

Dual-Band Observations of the Asymmetric Ring around CIDA 9A: Dead or Alive?

Authors: Daniel Harsono, Feng Long, Paola Pinilla, Alessia A. Rota, Carlo F. Manara, Gregory J. Herczeg, Doug Johnstone, Giovanni Rosotti, Giuseppe Lodato, Francois Menard, Marco Tazzari, Yangfan Shi

Abstract: While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data.… ▽ More While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

The population of young low-mass stars in Trumpler 14

Authors: Dominika Itrich, Leonardo Testi, Giacomo Beccari, Carlo F. Manara, Megan Reiter, Thomas Preibisch, Anna F. McLeod, Giovanni Rosotti, Ralf Klessen, Sergio Molinari, Patrick Hennebelle

Abstract: Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from out… ▽ More Massive star-forming regions are thought to be the most common birth environments in the Galaxy and the only birth places of very massive stars. Their presence in the stellar cluster alters the conditions within the cluster impacting at the same time the evolution of other cluster members. In principle, copious amounts of ultraviolet radiation produced by massive stars can remove material from outer parts of the protoplanetary disks around low- and intermediate-mass stars in the process of external photoevaporation, effectively reducing the planet-formation capabilities of those disks. Here, we present deep VLT/MUSE observations of low-mass stars in Trumpler 14, one of the most massive, young, and compact clusters in the Carina Nebula Complex. We provide spectral and stellar properties of 717 sources and based on the distribution of stellar ages derive the cluster age of $\sim$1~Myr. The majority of the stars in our sample have masses $\leqslant$1~$M_\odot$, what makes our spectroscopic catalogue the most deep to date in term of masses, and proves that detailed investigations of low-mass stars are possible in the massive but distant regions. Spectroscopic studies of low-mass members of the whole Carina Nebula Complex are missing. Our work provides an important step forward towards filling this gap and set the stage for follow-up investigation of accretion properties in Trumpler 14. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A, 27 pages, 28 figures

Journal ref: A&A 685, A100 (2024)

The time evolution of $M_{\mathrm{d}}/\dot M$ in protoplanetary discs as a way to disentangle between viscosity and MHD winds

Authors: Alice Somigliana, Leonardo Testi, Giovanni Rosotti, Claudia Toci, Giuseppe Lodato, Benoît Tabone, Carlo Manara, Marco Tazzari

Abstract: As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - m… ▽ More As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - mainly due to instrumental limitations. In this work, we present a viable method to distinguish between the viscous and MHD framework based on the different evolution of the distribution in the disk mass ($M_{\mathrm{d}}$) - accretion rate ($\dot M$) plane of a disk population. With a synergy of analytical calculations and 1D numerical simulations, performed with the population synthesis code \texttt{Diskpop}, we find that both mechanisms predict the spread of the observed ratio $M_{\mathrm{d}}/\dot M$ in a disk population to decrease over time; however, this effect is much less pronounced in MHD-dominated populations as compared to purely viscous populations. Furthermore, we demonstrate that this difference is detectable with the current observational facilities: we show that convolving the intrinsic spread with the observational uncertainties does not affect our result, as the observed spread in the MHD case remains significantly larger than in the viscous scenario. While the most recent data available show a better agreement with the wind model, ongoing and future efforts to obtain direct gas mass measurements with ALMA and ngVLA will cause a reassessment of this comparison in the near future. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 11 pages, 6 figures, 1 table. Published in ApJL

Journal ref: ApJL 954 (2023) L13

Constraining turbulence in protoplanetary discs using the gap contrast: an application to the DSHARP sample

Authors: E. Pizzati, G. P. Rosotti, B. Tabone

Abstract: Constraining the strength of gas turbulence in protoplanetary discs is an open problem that has relevant implications for the physics of gas accretion and planet formation. In this work, we gauge the amount of turbulence in 6 of the discs observed in the DSHARP programme by indirectly measuring the vertical distribution of their dust component. We employ the differences in the gap contrasts observ… ▽ More Constraining the strength of gas turbulence in protoplanetary discs is an open problem that has relevant implications for the physics of gas accretion and planet formation. In this work, we gauge the amount of turbulence in 6 of the discs observed in the DSHARP programme by indirectly measuring the vertical distribution of their dust component. We employ the differences in the gap contrasts observed along the major and the minor axes due to projection effects, and build a radiative transfer model to reproduce these features for different values of the dust scale heights. We find that (a) the scale heights that yield a better agreement with data are generally low ($\lesssim 4$ AU at a radial distance of $100$ AU), and in almost all cases we are only able to place upper limits on their exact values; these conclusions imply (assuming an average Stokes number of $\approx10^{-2}$) low turbulence levels of $α_{\rm SS}\lesssim10^{-3}-10^{-4}$; (b) for the 9 other systems we considered out of the DSHARP sample, our method yields no significant constraints on the disc vertical structure; we conclude that this is because these discs have either a low inclination or gaps that are not deep enough. Based on our analysis we provide an empirical criterion to assess whether a given disc is suitable to measure the vertical scale height. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: Accepted for publication in MNRAS. 13 pages + appendix, 12 figures

Radio multiwavelength analysis of the compact disk CX Tau: Presence of strong free-free variability or anomalous microwave emission

Authors: Pietro Curone, Leonardo Testi, Enrique Macias, Marco Tazzari, Stefano Facchini, Jonathan P. Williams, Cathie J. Clarke, Antonella Natta, Giovanni Rosotti, Claudia Toci, Giuseppe Lodato

Abstract: Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosp… ▽ More Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosphere. This work is focused on CX Tau, a ${\sim}0.4\,M_\odot$ star with an extended gas emission and a compact and apparently structureless dust disk, with an average millimeter flux when compared to Class II sources in Taurus. We present observations from the Karl G. Jansky Very Large Array (VLA) across four bands (between 9.0 mm and 6.0 cm) and combine them with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). This multiwavelength approach allows us to separate the dust continuum from other emissions. After isolating the dust thermal emission, we derived an upper limit of the dust disk extent at 1.3 cm which is consistent with theoretical predictions of a radial drift-dominated disk. Centimeter data show a peculiar behavior: deep observations at 6.0 cm did not detect the source, while at 1.3 cm the flux density is anomalously higher than adjacent bands. Intraband spectral indices suggest a dominant contribution from free-free emission, whereas gyro-synchrotron emission is excluded. To explain these observations, we propose a strong variability among the free-free emission with timescales shorter than a month. Another possible interpretation is the presence of anomalous microwave emission from spinning dust grains. △ Less

Submitted 22 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: Accepted for publication in A&A; 13 pages, 13 figures, 1 table

Journal ref: A&A 677, A118 (2023)

How large is a disk -- what do protoplanetary disk gas sizes really mean?

Authors: Leon Trapman, Giovanni Rosotti, Ke Zhang, Benoit Tabone

Abstract: It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instability or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successf… ▽ More It remains unclear what mechanism is driving the evolution of protoplanetary disks. Direct detection of the main candidates, either turbulence driven by magnetorotational instability or magnetohydrodynamical disk winds, has proven difficult, leaving the time evolution of the disk size as one of the most promising observables able to differentiate between these two mechanisms. But to do so successfully, we need to understand what the observed gas disk size actually traces. We studied the relation between $R_{\rm CO,\ 90\%}$, the radius that encloses 90% of the $^{12}$CO flux, and $R_c$, the radius that encodes the physical disk size, in order to provide simple prescriptions for conversions between these two sizes. For an extensive grid of thermochemical models we calculate $R_{\rm CO,\ 90\%}$ from synthetic observations and relate properties measured at this radius, such as the gas column density, to bulk disk properties, such as $R_c$ and the disk mass $M_{\rm disk}$. We found an empirical correlation between the gas column density at $R_{\rm CO,\ 90\%}$ and disk mass: $N_{\rm gas}(R_{\rm CO,\ 90\%}) \approx 3.73\times10^{21}(M_{\rm disk}/\mathrm{M}_{\odot})^{0.34}\ \mathrm{cm}^{-2}$. Using this correlation we derive an analytical prescription of $R_{\rm CO,\ 90\%}$ that only depends on $R_c$ and $M_{\rm disk}$. We derive $R_c$ for disks in Lupus, Upper Sco, Taurus and DSHARP, finding that disks in the older Upper Sco region are significantly smaller ($\langle R_c \rangle$ = 4.8 au) than disks in the younger Lupus and Taurus regions ($\langle R_c \rangle$ = 19.8 and 20.9 au, respectively). This temporal decrease in $R_c$ goes against predictions of both viscous and wind-driven evolution, but could be a sign of significant external photoevaporation having truncated disks in Upper Sco. △ Less

Submitted 14 July, 2023; originally announced July 2023.

Comments: 23 pages, 12 figures, accepted in ApJ

JWST reveals excess cool water near the snowline in compact disks, consistent with pebble drift

Authors: Andrea Banzatti, Klaus M. Pontoppidan, John Carr, Evan Jellison, Ilaria Pascucci, Joan Najita, Carlos E. Munoz-Romero, Karin I. Oberg, Anusha Kalyaan, Paola Pinilla, Sebastiaan Krijt, Feng Long, Michiel Lambrechts, Giovanni Rosotti, Gregory J. Herczeg, Colette Salyk, Ke Zhang, Edwin Bergin, Nick Ballering, Michael R. Meyer, Simon Bruderer, the JDISCS collaboration

Abstract: Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles dri… ▽ More Previous analyses of mid-infrared water spectra from young protoplanetary disks observed with the Spitzer-IRS found an anti-correlation between water luminosity and the millimeter dust disk radius observed with ALMA. This trend was suggested to be evidence for a fundamental process of inner disk water enrichment, used to explain properties of the Solar System 40 years ago, in which icy pebbles drift inward from the outer disk and sublimate after crossing the snowline. Previous analyses of IRS water spectra, however, were uncertain due to the low spectral resolution that blended lines together. We present new JWST-MIRI spectra of four disks, two compact and two large with multiple radial gaps, selected to test the scenario that water vapor inside the snowline is regulated by pebble drift. The higher spectral resolving power of MIRI-MRS now yields water spectra that separate individual lines, tracing upper level energies from 900 K to 10,000 K. These spectra clearly reveal excess emission in the low-energy lines in compact disks, compared to the large disks, demonstrating an enhanced cool component with $T \approx$ 170-400 K and equivalent emitting radius $R_{\rm{eq}}\approx$ 1-10 au. We interpret the cool water emission as ice sublimation and vapor diffusion near the snowline, suggesting that there is indeed a higher inwards mass flux of icy pebbles in compact disks. Observation of this process opens up multiple exciting prospects to study planet formation chemistry in inner disks with JWST. △ Less

Submitted 3 September, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ Letters

The Effect of Dust Evolution and Traps on Inner Disk Water Enrichment

Authors: Anusha Kalyaan, Paola Pinilla, Sebastiaan Krijt, Andrea Banzatti, Giovanni Rosotti, Gijs D. Mulders, Michiel Lambrechts, Feng Long, Gregory J. Herczeg

Abstract: Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance wit… ▽ More Substructures in protoplanetary disks can act as dust traps that shape the radial distribution of pebbles. By blocking the passage of pebbles, the presence of gaps in disks may have a profound effect on pebble delivery into the inner disk, crucial for the formation of inner planets via pebble accretion. This process can also affect the delivery of volatiles (such as H$_2$O) and their abundance within the water snow line region (within a few au). In this study, we aim to understand what effect the presence of gaps in the outer gas disk may have on water vapor enrichment in the inner disk. Building on previous work, we employ a volatile-inclusive multi-Myr disk evolution model that considers an evolving ice-bearing drifting dust population, sensitive to dust-traps, which loses its icy content to sublimation upon reaching the snow line. We find that vapor abundance in the inner disk is strongly affected by fragmentation velocity (v$_{\rm f}$) and turbulence, which control how intense vapor enrichment from pebble delivery is, if present, and how long it may last. Generally, for disks with low to moderate turbulence ($α$ $\le$ 1 $\times$ 10$^{-3}$) and for a range of v$_{\rm f}$, radial location, and gap depth (especially that of the innermost gaps), can significantly alter enrichment. Shallow inner gaps may continuously leak material from beyond it, despite the presence of additional deep outer gaps. We finally find that the for realistic v$_{\rm f}$ ($\le$ 10 m s$^{-1}$), presence of gaps is more important than planetesimal formation beyond the snow line in regulating pebble and volatile delivery into the inner disk. △ Less

Submitted 4 July, 2023; originally announced July 2023.

Comments: 24 Pages, 15 Figures, Accepted to ApJ

The distribution of accretion rates as a diagnostic of protoplanetary disc evolution

Authors: R. Alexander, G. Rosotti, P. J. Armitage, G. J. Herczeg, C. F. Manara, B. Tabone

Abstract: We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent ("viscous") transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetised disc winds. We find that a homogeneous sample of $\gtrsim$300 observed accretion rates would be sufficient to d… ▽ More We show that the distribution of observed accretion rates is a powerful diagnostic of protoplanetary disc physics. Accretion due to turbulent ("viscous") transport of angular momentum results in a fundamentally different distribution of accretion rates than accretion driven by magnetised disc winds. We find that a homogeneous sample of $\gtrsim$300 observed accretion rates would be sufficient to distinguish between these two mechanisms of disc accretion at high confidence, even for pessimistic assumptions. Current samples of T Tauri star accretion rates are not this large, and also suffer from significant inhomogeneity, so both viscous and wind-driven models are broadly consistent with the existing observations. If accretion is viscous, the observed accretion rates require low rates of disc photoevaporation ($\lesssim$$10^{-9}$M$_{\odot}$yr$^{-1}$). Uniform, homogeneous surveys of stellar accretion rates can therefore provide a clear answer to the long-standing question of how protoplanetary discs accrete. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

Distinguishing Magnetized Disc Winds from Turbulent Viscosity through Substructure Morphology in Planet-forming Discs

Authors: Yinhao Wu, Yi-Xian Chen, Haochang Jiang, Ruobing Dong, Enrique Macías, Min-Kai Lin, Giovanni P. Rosotti, Vardan Elbakyan

Abstract: The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morpholog… ▽ More The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morphology of planet-forming discs in the ALMA continuum. We model the gap-opening process of planets in discs with both viscous evolution and wind-driven accretion by 2D multi-fluid hydrodynamical simulations. Our results show that gap-opening planets in wind-driven accreting discs generate characteristic substructures that differ from those in purely viscous discs. Specifically, we demonstrate that discs, where wind-driven accretion dominates the production of substructures, exhibit significant asymmetries. Based on the diverse outputs of mock images in the ALMA continuum, we roughly divide the planet-induced features into four regimes (moderate-viscosity dominated, moderate-wind dominated, strong-wind dominated, inviscid). The classification of these regimes sets up a potential method to constrain the strength of magnetized disc wind and viscosity based on the observed gap and ring morphology. We discuss the asymmetry feature in our mock images and its potential manifestation in ALMA observations. △ Less

Submitted 27 May, 2023; v1 submitted 14 May, 2023; originally announced May 2023.

Comments: 11 pages, 6 figures, accepted by MNRAS. Welcome any comments and suggestions!

The Disc Miner II: Revealing Gas substructures and Kinematic signatures from Planet-disc interaction through line profile analysis

Authors: Andres F. Izquierdo, Leonardo Testi, Stefano Facchini, Giovanni P. Rosotti, Ewine van Dishoeck, Lisa Wölfer, Teresa Paneque-Carreño

Abstract: [Abridged] The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the DISCMINER modelling fr… ▽ More [Abridged] The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the DISCMINER modelling framework, in three bright CO lines, 12CO, 13CO, and C18O $J=2-1$, to provide a comprehensive summary of the kinematic and column density substructures that planets might be actively sculpting in these discs. We find convincing evidence for the presence of four giant planets located at wide orbits in three of the discs in the sample: two around HD 163296, one in MWC 480, and one in AS 209. One of the planet candidates in HD 163296, P94, previously associated with velocity signatures detected in lower velocity resolution 12CO data, is confirmed and linked to localised velocity and line width perturbations in 13CO and C18O too. We highlight that line widths are also powerful tracers of planet-forming sites as they are sensitive to turbulent motions triggered by planet-disc interactions. In MWC 480, we identified non-axisymmetric line width enhancements around the radial separation of candidate planet-driven buoyancy spirals, which we used to narrow the location of the possible planet to an orbital radius of $R=245$ au and $\rm{PA}=193^\circ$. In the disc of AS 209, we found excess 12CO line widths centred at $R=210$ au, $\rm{PA}=151^\circ$, spanning around the immediate vicinity of a circumplanetary disc candidate proposed previously, which further supports its presence. Our simultaneous analysis of multiple tracers and observables aims to lay the groundwork for robust studies of molecular line properties focused on the search for young planets in discs. △ Less

Submitted 15 May, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A

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

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Diverse outcomes of binary-disk interactions

Authors: Yapeng Zhang, Christian Ginski, Jane Huang, Alice Zurlo, Hervé Beust, Jaehan Bae, Myriam Benisty, Antonio Garufi, Michiel R. Hogerheijde, Rob G. van Holstein, Matthew Kenworthy, Maud Langlois, Carlo F. Manara, Paola Pinilla, Christian Rab, Álvaro Ribas, Giovanni P. Rosotti, Jonathan Williams

Abstract: Circumstellar disks do not evolve in isolation, as about half of solar-type stars were born in binary or multiple systems. Resolving disks in binary systems provides the opportunity to examine the influence of stellar companions on the outcomes of planet formation. We aim to investigate and compare disks in stellar multiple systems with near-infrared scattered-light imaging as part of the Disk Evo… ▽ More Circumstellar disks do not evolve in isolation, as about half of solar-type stars were born in binary or multiple systems. Resolving disks in binary systems provides the opportunity to examine the influence of stellar companions on the outcomes of planet formation. We aim to investigate and compare disks in stellar multiple systems with near-infrared scattered-light imaging as part of the Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS) program. We used polarimetric differential imaging with SPHERE/IRDIS at the VLT to search for scattered light from the circumstellar disks in three multiple systems, CHX 22, S CrA, and HP Cha. We performed astrometric and orbit analyses for the stellar companions using archival HST, VLT/NACO, and SPHERE data. Combined with the age and orbital constraints, the observed disk structures provide insights into the evolutionary history and the impact of the stellar companions. The small grains in CHX 22 form a tail-like structure surrounding the close binary, which likely results from a close encounter and capture of a cloudlet. S CrA shows intricate structures (tentative ringed and spiral features) in the circumprimary disk as a possible consequence of perturbations by companions. The circumsecondary disk is truncated and connected to the primary disk via a streamer, suggesting tidal interactions. In HP Cha, the primary disk is less disturbed and features a tenuous streamer, through which the material flows towards the companions. The comparison of the three systems spans a wide range of binary separation (50 - 500 au) and illustrates the decreasing influence on disk structures with the distance of companions. This agrees with the statistical analysis of exoplanet population in binaries, that planet formation is likely obstructed around close binary systems, while it is not suppressed in wide binaries. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 19 pages, 6 figures, accpeted for publication in A&A

Journal ref: A&A 672, A145 (2023)

Empirical constraints on turbulence in proto-planetary discs

Authors: Giovanni P. Rosotti

Abstract: Proto-planetary discs, the birth environment of planets, are an example of a structure commonly found in astrophysics, accretion discs. Identifying the mechanism responsible for accretion is a long-standing problem, dating back several decades. The common picture is that accretion is a consequence of turbulence, with several instabilities proposed for its origin. While traditionally this field use… ▽ More Proto-planetary discs, the birth environment of planets, are an example of a structure commonly found in astrophysics, accretion discs. Identifying the mechanism responsible for accretion is a long-standing problem, dating back several decades. The common picture is that accretion is a consequence of turbulence, with several instabilities proposed for its origin. While traditionally this field used to be a purely theoretical endeavour, the landscape is now changing thanks mainly to new observational facilities such as the ALMA radio interferometer. Thanks to large improvements in spatial and spectral resolution and sensitivity (which have enabled the study of disc substructure, kinematics and surveys of large disc populations), multiple techniques have been devised to observationally measure the amount of turbulence in discs. This review summarises these techniques, ranging from attempts at direct detection of turbulence from line broadening, to more indirect approaches that rely on properties of the dust or consider the evolution of global disc properties (such as masses, radii and accretion rates) for large samples, and what their findings are. Multiple lines of evidence suggest that discs are in fact not as turbulent as thought one decade ago. On the other hand, direct detection of turbulence in some discs and the finite radial extent of dust substructures and in some cases the finite vertical extent strongly indicate that turbulence must be present at some level in proto-planetary discs. It is still an open question whether this amount of turbulence is enough to power accretion or if this is instead driven by other mechanisms, such as MHD winds. △ Less

Submitted 3 April, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 24 pages, 7 figures. Accepted for publication on New Astronomy Reviews; small changes following community feedback in v2

Dust dynamics in planet-forming discs in binary systems

Authors: Francesco Zagaria, Giovanni P. Rosotti, Richard D. Alexander, Cathie J. Clarke

Abstract: In multiple stellar systems interactions among the companion stars and their discs affect planet formation. In the circumstellar case tidal truncation makes protoplanetary discs smaller, fainter and less long-lived than those evolving in isolation, thereby reducing the amount of material (gas and dust) available to assemble planetary embryos. On the contrary, in the circumbinary case the reduced a… ▽ More In multiple stellar systems interactions among the companion stars and their discs affect planet formation. In the circumstellar case tidal truncation makes protoplanetary discs smaller, fainter and less long-lived than those evolving in isolation, thereby reducing the amount of material (gas and dust) available to assemble planetary embryos. On the contrary, in the circumbinary case the reduced accretion can increase the disc lifetime, with beneficial effects on planet formation. In this chapter we review the main observational results on discs in multiple stellar systems and discuss their possible explanations, focusing on recent numerical simulations, mainly dealing with dust dynamics and disc evolution. Finally, some open issues and future research directions are examined. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: Accepted for publication in EPJ Plus Focus Point on Environmental and Multiplicity Effects on Planet Formation G. Lodato and C.F. Manara (Guest editors)

The Mass Accretion Rate and Stellar Properties in Class I Protostars

Authors: Eleonora Fiorellino, Lukasz Tychoniec, Fernando Cruz-Saenz de Miera, Simone Antoniucci, Agnes Kospal, Carlo F. Manara, Brunella Nisini, Giovanni Rosotti

Abstract: Stars collect most of their mass during the protostellar stage, yet the accretion luminosity and stellar parameters, which are needed to compute the mass accretion rate, are poorly constrained for the youngest sources. The aim of this work is to fill this gap, computing the stellar properties and the accretion rates for a large sample of Class I protostars located in nearby (< 500 pc) star-forming… ▽ More Stars collect most of their mass during the protostellar stage, yet the accretion luminosity and stellar parameters, which are needed to compute the mass accretion rate, are poorly constrained for the youngest sources. The aim of this work is to fill this gap, computing the stellar properties and the accretion rates for a large sample of Class I protostars located in nearby (< 500 pc) star-forming regions and analysing their interplay. We used a self-consistent method to provide accretion and stellar parameters using SED modeling and veiling information from near-IR observations, when possible. We calculated accretion and stellar properties for the first time for 50 young stars. We focused our analysis on the 39 confirmed protostars, finding that their mass accretion rate varies between about 10^(-8) and about 10^(-4) Msun/yr in a stellar mass range between about 0.1 and 3 Msun. We find systematically larger mass accretion rates for our Class I sample than in Class II objects. Although the mass accretion rate we found is high, it still suggests that either stars collect most of its mass before Class I stage, or eruptive accretion is needed during the overall protostellar phase. Indeed, our results suggest that for a large number of protostars the disk can be unstable, which can result in accretion bursts and disk fragmentation in the past or in the future. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Importance of source structure on complex organics emission III. Effect of disks around massive protostars

Authors: P. Nazari, B. Tabone, G. P. Rosotti

Abstract: Complex organic molecules are only detected toward a fraction of high-mass protostars. The goal of this work is to investigate whether high-mass disks can explain the lack of methanol emission from some massive protostellar systems. We consider an envelope-only and an envelope-plus-disk model and use RADMC-3D to calculate the methanol emission. High and low millimeter (mm) opacity dust are conside… ▽ More Complex organic molecules are only detected toward a fraction of high-mass protostars. The goal of this work is to investigate whether high-mass disks can explain the lack of methanol emission from some massive protostellar systems. We consider an envelope-only and an envelope-plus-disk model and use RADMC-3D to calculate the methanol emission. High and low millimeter (mm) opacity dust are considered for both models separately and the methanol abundance is parameterized. Viscous heating is included due to the high accretion rates of these objects in the disk. In contrast with low-mass protostars, the presence of a disk does not significantly affect the temperature structure and methanol emission. The shadowing effect of the disk is not as important for high-mass objects and the disk mid-plane is hot because of viscous heating, which is effective due to the high accretion rates. Consistent with observations of infrared absorption lines toward high-mass protostars, we find a vertical temperature inversion, i.e. higher temperatures in the disk mid-plane than the disk surface, at radii < 50au for the models with $L=10^4$ L$_{\odot}$ and large mm opacity dust as long as the envelope mass is >550 M$_{\odot}$. The large observed scatter in methanol emission from massive protostars can be mostly explained toward lower luminosity objects with the envelope-plus-disk models including low and high mm opacity dust. The methanol emission variation toward sources with high luminosities cannot be explained by models with or without a disk. However, the $L/M$ of these objects suggest that they could be associated with hypercompact/ultracompact HII regions. Therefore, the low methanol emission toward the high-luminosity sources can be explained by them hosting an HII region where methanol is absent. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 25 pages, 24 figures, Accepted for publication in A&A

Journal ref: A&A 671, A107 (2023)

An analytic solution to measure the gas size in protoplanetary discs in the viscous self-similar scenario

Authors: Claudia Toci, Giuseppe Lodato, Francesco Gerardo Livio, Giovanni Rosotti, Leon Trapman

Abstract: In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical soluti… ▽ More In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical solution for the evolution of the disc radii in viscously evolving protoplanetary discs using $^{12}$CO as a tracer, under the assumption that the $^{12}$CO radius is the radius where the surface density of the disc is equal to the threshold for CO photo-dissociation. We discuss the properties of the solution and the limits of its applicability as a simple numerical prescription to evaluate the observed disc radii of populations of discs. Our results suggest that, in addition to photo-dissociation, also freeze out plays an important role in setting the disc size. We find an effective reduction of the CO abundance by about two orders of magnitude at the location of CO photo-dissociation, which however should not be interpreted as the bulk abundance of CO in the disc. The use of our analytical solution will allow to compute disc sizes for large quantities of models without using expensive computational resources such as radiative transfer calculations. △ Less

Submitted 10 November, 2022; v1 submitted 26 October, 2022; originally announced October 2022.

Comments: 5 pages, 3 figures

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)

The relation between the Mass Accretion Rate and the Disk Mass in Class I Protostars

Authors: Eleonora Fiorellino, Lukasz Tychoniec, Carlo F. Manara, Giovanni Rosotti, Simone Antoniucci, Fernando Cruz-Saenz de Miera, Agnes Kospal, Brunella Nisini

Abstract: The evidence of a relation between the mass accretion rate and the disk mass is established for young, Class II pre-main sequence stars. This observational result opened an avenue to test theoretical models and constrain the initial conditions of the disk formation, fundamental in the understanding of the emergence of planetary systems. However, it is becoming clear that the planet formation start… ▽ More The evidence of a relation between the mass accretion rate and the disk mass is established for young, Class II pre-main sequence stars. This observational result opened an avenue to test theoretical models and constrain the initial conditions of the disk formation, fundamental in the understanding of the emergence of planetary systems. However, it is becoming clear that the planet formation starts even before the Class II stage, in disks around Class 0 and I protostars. We show for the first time evidence for a correlation between the mass accretion rate and the disk mass for a large sample of Class I young stars located in nearby (< 500 pc) star-forming regions. We fit our sample, finding that the Class I objects relation has a slope flatter than Class II stars, and have higher mass accretion rates and disk masses. The results are put in context of the disk evolution models. △ Less

Submitted 9 September, 2022; originally announced September 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)

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

Authors: Brodie J. Norfolk, Christophe Pinte, Josh Calcino, Iain Hammond, Nienke van der Marel, Daniel J. Price, Sarah T. Maddison, Valentin Christiaens, Jean-Francois Gonzalez, Dori Blakely, Giovanni Rosotti, Christian Ginski

Abstract: Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner sp… ▽ More Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL, 14 pages, 10 figures

PENELLOPE III. The peculiar accretion variability of XX Cha and its impact on the observed spread of accretion rates

Authors: R. A. B. Claes, C. F. Manara, R. Garcia-Lopez, A. Natta, M. Fang, Z. P. Fockter, P. Ábrahám, J. M. Alcalá, J. Campbell-White, A. Caratti o Garatti, E. Covino, D. Fedele, A. Frasca, J. F. Gameiro, G. J. Herczeg, Á. Kóspál, M. G. Petr-Gotzens, G. Rosotti, L. Venuti, G. Zsidi

Abstract: The processes regulating protoplanetary disk evolution are constrained by studying how mass accretion rates scale with stellar and disk properties. The spread in these relations can be used as a constraint to the models of disk evolution, but only if the impact of accretion variability is correctly accounted for. While the effect of variability might be substantial in the embedded phases of star f… ▽ More The processes regulating protoplanetary disk evolution are constrained by studying how mass accretion rates scale with stellar and disk properties. The spread in these relations can be used as a constraint to the models of disk evolution, but only if the impact of accretion variability is correctly accounted for. While the effect of variability might be substantial in the embedded phases of star formation, it is often considered limited at later stages. Here we report on the observed large variation in the accretion rate for one target, XX Cha, and we discuss the impact on population studies of classical T Tauri stars. The mass accretion rate determined by fitting the UV-to-near-infrared spectrum in recent X-Shooter observations is compared with the one measured with the same instrument 11 years before. XX Cha displays an accretion variability of almost 2 dex between 2010 and 2021. Although the timescales on which this variability happens are uncertain, XX Cha displays an extreme accretion variability for a classical T Tauri star. If such behavior is common among classical T Tauri stars, possibly on longer timescales than previously probed, it could be relevant for discussing the disk evolution models constrained by the observed spread in accretion rates. Finally, we remark that previous studies of accretion variability based on spectral lines may have underestimated the variability of some targets. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A letters

Journal ref: A&A 664, L7 (2022)

Distribution of solids in the rings of the HD 163296 disk: a multiwavelength study

Authors: G. Guidi, A. Isella, L. Testi, C. J. Chandler, H. B. Liu, H. M. Schmid, G. Rosotti, C. Meng, J. Jennings, J. P. Williams, J. M. Carpenter, I. de Gregorio-Monsalvo, H. Li, S. F. Liu, S. Ortolani, S. P. Quanz, L. Ricci, M. Tazzari

Abstract: In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic… ▽ More In this paper we analyze new observations from ALMA and VLA, at a high angular resolution corresponding to 5 - 8 au, of the protoplanetary disk around HD 163296 to determine the dust spatial distribution and grain properties. We fit the spectral energy distribution as a function of the radius at five wavelengths from 0.9 to 9\,mm, using a simple power law and a physical model based on an analytic description of radiative transfer that includes isothermal scattering. We considered eight dust populations and compared the models' performance using Bayesian evidence. Our analysis shows that the moderately high optical depth ($τ$>1) at $λ\leq$ 1.3 mm in the dust rings artificially lower the millimeter spectral index, which should therefore not be considered as a reliable direct proxy of the dust properties and especially the grain size. We find that the outer disk is composed of small grains on the order of 200 $μ$m with no significant difference between rings at 66 and 100 au and the adjacent gaps, while in the innermost 30 au, larger grains ($\geq$mm) could be present. We show that the assumptions on the dust composition have a strong impact on the derived surface densities and grain size. In particular, increasing the porosity of the grains to 80\% results in a total dust mass about five times higher with respect to grains with 25\% porosity. Finally, we find that the derived opacities as a function of frequency deviate from a simple power law and that grains with a lower porosity seem to better reproduce the observations of HD163296. While we do not find evidence of differential trapping in the rings of HD163296, our overall results are consistent with the postulated presence of giant planets affecting the dust temperature structure and surface density, and possibly originating a second-generation dust population of small grains. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 27 pages, 28 figures, Accepted for publication in Astronomy and Astrophysics (A&A)

Journal ref: A&A 664, A137 (2022)

Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments

Authors: Andrew J. Winter, Cathie J. Clarke, Giovanni Rosotti, Mirek Giersz

Abstract: Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularis… ▽ More Short period, massive planets, known as hot Jupiters (HJs), have been discovered around $\sim 1$ percent of local field stars. The inward migration necessary to produce HJs may be `low eccentricity', due to torques in the primordial disc, or `high eccentricity' (HEM). The latter involves exciting high orbital eccentricity, allowing sufficiently close passages with the host star to raise circularising tides in the planet. We present an analytic framework for quantifying the role of dynamical encounters in high density environments during HEM. We show that encounters can enhance or suppress HEM, depending on the local stellar density and the initial semi-major axis $a_0$. For moderate densities, external perturbations can excite large eccentricities that allow a planet to circularise over the stellar lifetime. At extremely high densities, these perturbations can instead result in tidal disruption of the planet, thus yielding no HJ. This may explain the apparent excess of HJs in M67 compared with their local field star abundance versus their apparent deficit in 47 Tuc. Applying our analytic framework, we demonstrate that for an initial massive planet population similar to the field, the expected HJ occurrence rate in 47 Tuc is $f_\mathrm{HJ}=2.2\times 10^{-3}$, which remains consistent with present constraints. Future large (sample sizes $\gtrsim 10^5$) or sensitive transit surveys of stars in globular clusters are required to refute the hypothesis that the initial planet population is similar to the solar neighbourhood average. Non-detection in such surveys would have broad consequences for planet formation theory, implying planet formation rates in globular clusters must be suppressed across a wide range of $a_0$. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS - 28 pages, 17 figures

Gap opening by planets in discs with magnetised winds

Authors: Vardan Elbakyan, Yinhao Wu, Sergei Nayakshin, Giovanni Rosotti

Abstract: Planets open deep gaps in protoplanetary discs when their mass exceeds a gap opening mass, $M_{\rm gap}$. We use one- and two-dimensional simulations to study planet gap opening in discs with angular momentum transport powered by MHD disc winds. We parameterise the efficiency of the MHD disc wind angular momentum transport through a dimensionless parameter $α_{\rm dw}$, which is an analogue to the… ▽ More Planets open deep gaps in protoplanetary discs when their mass exceeds a gap opening mass, $M_{\rm gap}$. We use one- and two-dimensional simulations to study planet gap opening in discs with angular momentum transport powered by MHD disc winds. We parameterise the efficiency of the MHD disc wind angular momentum transport through a dimensionless parameter $α_{\rm dw}$, which is an analogue to the turbulent viscosity $α_{\rm v}$. We find that magnetised winds are much less efficient in counteracting planet tidal torques than turbulence is. For discs with astrophysically realistic values of $α_{\rm dw}$, $M_{\rm gap}$ is always determined by the residual disc turbulence, and is a factor of a few to ten smaller than usually obtained for viscous discs. We introduce a gap opening criterion applicable for any values of $α_{\rm v}$ and $α_{\rm dw}$ that may be useful for planet formation population synthesis. We show that in discs powered by magnetised winds, growing planets detach from the disc at planet masses below $\sim 0.1M_{\rm Jup}$ inside 10 AU. This promotes formation of super-Earth planets rather than gas giants in this region, in particular precluding formation of hot jupiters in situ. On larger scales, ALMA gap opening planet candidates may be less massive than currently believed. Future high-resolution observations with instruments such as the extended ALMA, ngVLA, and SKA are likely to show abundant narrow annular features at $R < 10$ AU due to ubiquitous super-Earth planets. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Comments: 14 pages, 9 figures, accepted to MNRAS

Super-resolution trends in the ALMA Taurus survey: Structured inner discs and compact discs

Authors: Jeff Jennings, Marco Tazzari, Cathie J. Clarke, Richard A. Booth, Giovanni P. Rosotti

Abstract: The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank')… ▽ More The 1.33 mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (>~55 au), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of nonparametric visibility fits for data at the survey's 0.12" resolution. Then we use the nonparametric model in Frankenstein ('frank') to identify new substructure in three compact and seven extended sources. Among the new features we identify three trends: a higher occurrence rate of substructure in the survey's compact discs than previously seen, underresolved (potentially azimuthally asymmetric) substructure in the innermost disc of extended sources, and a 'shoulder' on the trailing edge of a ring in discs with strong depletion at small radii. Noting the shoulder morphology is present in multiple discs observed at higher resolution, we postulate it is tracing a common physical mechanism. We further demonstrate how a super-resolution frank brightness profile is useful in motivating an accurate parametric model, using the highly structured source DL Tau in which frank finds two new rings. Finally we show that sparse (u, v) plane sampling may be masking the presence of substructure in several additional compact survey sources. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: Accepted to MNRAS

Distributions of gas and small and large grains in the LkH$α\,330$ disk trace a young planetary system

Authors: P. Pinilla, M. Benisty, N. T. Kurtovic, J. Bae, R. Dong, Z. Zhu, S. Andrews, J. Carpenter, C. Ginski, J. Huang, A. Isella, L. Pérez, L. Ricci, G. Rosotti, M. Villenave, D. Wilner

Abstract: [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''… ▽ More [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$α\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction. △ Less

Submitted 25 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A, minor changes after language edition

Journal ref: A&A 665, A128 (2022)

Inward and outward migration of massive planets: moving towards a stalling radius

Authors: C. E. Scardoni, C. J. Clarke, G. P. Rosotti, R. A. Booth, R. D. Alexander, E. Ragusa

Abstract: Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simu… ▽ More Recent studies on the planet-dominated regime of Type II migration showed that, contrary to the conventional wisdom, massive planets can migrate outwards. Using `fixed-planet' simulations these studies found a correlation between the sign of the torques acting on the planet and the parameter $K'$ (which describes the depth of the gap carved by the planet in the disc). We perform `live-planet' simulations exploring a range of $K'$ and disc mass values to test and extend these results. The excitation of planet eccentricity in live-planet simulations breaks the direct dependence of migration rate (rate of change of semi-major axis) on the torques imposed, an effect that `fixed-planet' simulations cannot treat. By disentangling the contribution to the torque due to the semi-major axis evolution from that due to the eccentricity evolution, we recover the relation between the magnitude and sign of migration and $K'$ and argue that this relation may be better expressed in terms of the related gap depth parameter $K$. We present a toy model in which the sign of planetary migration changes at a limiting value of $K$, through which we explore planets' migration in viscously evolving discs. The existence of the torque reversal shapes the planetary system's architecture by accumulating planets either at the stalling radius or in a band around it (defined by the interplay between the planet migration and the disc evolution). In either case, planets pile up in the area $1-10$ au, disfavouring hot Jupiter formation through Type II migration in the planet-dominated regime. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 14 figures, accepted for publication in MNRAS

On the time evolution of the $M_{\rm d} - M_\star$ and $\dot M - M_\star$ correlations for protoplanetary discs: the viscous timescale increases with stellar mass

Authors: Alice Somigliana, Claudia Toci, Giovanni Rosotti, Giuseppe Lodato, Marco Tazzari, Carlo Manara, Leonardo Testi, Federico Lepri

Abstract: Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their… ▽ More Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their subsequent evolution is still uncertain. We conduct a theoretical analysis of the impact of viscous evolution on power-law initial conditions for a population of protoplanetary discs. We find that, for evolved populations, viscous evolution enforces the two correlations to have the same slope, $λ_{\mathrm{m}}$ = $λ_{\mathrm{acc}}$, and that this limit is uniquely determined by the initial slopes $λ_{\mathrm{m}, 0}$ and $λ_{\mathrm{acc}, 0}$. We recover the increasing trend claimed from the observations when the difference in the initial values, $δ_0 = λ_{\mathrm{m}, 0} - λ_{\mathrm{acc}, 0}$, is larger than $1/2$; moreover, we find that this increasing trend is a consequence of a positive correlation between the viscous timescale and the stellar mass. We also present the results of disc population synthesis numerical simulations, that allow us to introduce a spread and analyse the effect of sampling, which show a good agreement with our analytical predictions. Finally, we perform a preliminary comparison of our numerical results with observational data, which allows us to constrain the parameter space of the initial conditions to $λ_{\mathrm{m}, 0} \in [1.2, 2.1]$, $λ_{\mathrm{acc}, 0} \in [0.7, 1.5]$. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS. 14 pages, 8 figures

Modelling the secular evolution of proto-planetary disc dust sizes -- A comparison between the viscous and magnetic wind case

Authors: Francesco Zagaria, Giovanni P. Rosotti, Cathie J. Clarke, Benoît Tabone

Abstract: For many years proto-planetary discs have been thought to evolve viscously: angular momentum redistribution leads to accretion and outward disc spreading. Recently, the hypothesis that accretion is due, instead, to angular momentum removal by magnetic winds gained new popularity: no disc spreading is expected in this case. In this paper, we run several one-dimensional gas and dust simulations to m… ▽ More For many years proto-planetary discs have been thought to evolve viscously: angular momentum redistribution leads to accretion and outward disc spreading. Recently, the hypothesis that accretion is due, instead, to angular momentum removal by magnetic winds gained new popularity: no disc spreading is expected in this case. In this paper, we run several one-dimensional gas and dust simulations to make predictions on the time-evolution of disc sizes \textit{in the dust} and to assess whether they can be used to understand how discs evolve. We show that viscous and magnetic wind models have very different dust disc radii. In particular, MHD wind models are compact and their sizes either remain constant or decrease with time. On the contrary, discs become larger with time in the viscous case (when $α\gtrsim10^{-3}$). Although current observations lack enough sensitivity to discriminate between these two scenarios, higher-sensitivity surveys could be fruitful to this goal on a $1\,{\rm to}\,10\,{\rm Myr}$ age range. When compared with the available ALMA Band~7 data, both viscous and magnetic wind models are compatible with the observationally-inferred dust radii in Lupus, Chamaeleon~I and Upper Sco. Furthermore, in the drift-dominated regime, the size-luminosity correlation is reproduced in Lupus, both in Band~7 and 3, while in Upper Sco a different slope than in the data is predicted. Sub-structures (potentially undetected) can explain several outliers with large observed sizes. Higher-angular-resolution observations will be helpful to test our predictions in the case of more compact discs, expected in both frameworks, particularly at the age of Upper Sco. △ Less

Submitted 22 May, 2022; originally announced May 2022.

Comments: 19 pages, 10 figures; accepted for publication in MNRAS

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

Authors: Charles J. Law, Sage Crystian, Richard Teague, Karin I. Öberg, Evan A. Rich, Sean M. Andrews, Jaehan Bae, Kevin Flaherty, Viviana V. Guzmán, Jane Huang, John D. Ilee, Joel H. Kastner, Ryan A. Loomis, Feng Long, Laura M. Pérez, Sebastián Pérez, Chunhua Qi, Giovanni P. Rosotti, Dary Ruíz-Rodríguez, Takashi Tsukagoshi, David J. Wilner

Abstract: High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure o… ▽ More High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

Unveiling the outer dust disc of TW Hya with deep ALMA observations

Authors: John D. Ilee, Catherine Walsh, Jeff Jennings, Richard A. Booth, Giovanni P. Rosotti, Richard Teague, Takashi Tsukagoshi, Hideko Nomura

Abstract: The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALM… ▽ More The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

ALMA constraints on assembly of Core Accretion planets

Authors: Sergei Nayakshin, Vardan Elbakyan, Giovanni Rosotti

Abstract: Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$α$ emission, hint that young massive planets are abundant at wide separations in protoplanetary discs. Here, we show how these observations can probe the runaway phase of planetary growth in the Core Accretion theory. Planets in this phase have the right range of masses to account for the predominantly moderate cont… ▽ More Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$α$ emission, hint that young massive planets are abundant at wide separations in protoplanetary discs. Here, we show how these observations can probe the runaway phase of planetary growth in the Core Accretion theory. Planets in this phase have the right range of masses to account for the predominantly moderate contrast gaps and rings seen in ALMA observations. However, we find that these planets gain mass and migrate inward very rapidly. As a result, the phase when they could produce gaps with properties similar to those observed is very short, i.e., $t_{\rm gap} \lesssim 0.1$~Myr, independently of the disc viscosity parameter. This would require many tens to hundreds of gas giant planets to be born per ALMA system, violating the available mass budget of solids in realistic discs. This also predicts preponderance of discs with very wide gaps or complete inner disc holes, which is not observed. We show that suppression of both planet accretion and migration by a factor of at least ten is a possible solution to these serious problems. Future population synthesis models of planet formation should aim to address both exoplanetary data of older discless planetary systems and ALMA discs with embedded planets in one framework. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 17 pages, 10 figures, Accepted by MNRAS

Gas Disk Sizes from CO Line Observations: A Test of Angular Momentum Evolution

Authors: Feng Long, Sean M. Andrews, Giovanni Rosotti, Daniel Harsono, Paola Pinilla, David J. Wilner, Karin I Öberg, Richard Teague, Leon Trapman, Benoît Tabone

Abstract: The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emis… ▽ More The size of a disk encodes important information about its evolution. Combining new Submillimeter Array (SMA) observations with archival Atacama Large Millimeter Array (ALMA) data, we analyze mm continuum and CO emission line sizes for a sample of 44 protoplanetary disks around stars with masses of 0.15--2\,$M_{\odot}$ in several nearby star-forming regions. Sizes measured from $^{12}$CO line emission span from 50 to 1000\,au. This range could be explained by viscous evolution models with different $α$ values (mostly of $10^{-4}-10^{-3}$) and/or a spread of initial conditions. The CO sizes for most disks are also consistent with MHD wind models that directly remove disk angular momentum, but very large initial disk sizes would be required to account for the very extended CO disks in the sample. As no CO size evolution is observed across stellar ages of 0.5--20\,Myr in this sample, determining the dominant mechanism of disk evolution will require a more complete sample for both younger and more evolved systems. We find that the CO emission is universally more extended than the continuum emission by an average factor of $2.9\pm1.2$. The ratio of the CO to continuum sizes does not show any trend with stellar mass, mm continuum luminosity, or the properties of substructures. The GO Tau disk has the most extended CO emission in this sample, with an extreme CO to continuum size ratio of 7.6. Seven additional disks in the sample show high size ratios ($\gtrsim4$) that we interpret as clear signs of substantial radial drift. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Importance of source structure on complex organics emission II. Can disks explain lack of methanol emission from some low-mass protostars?

Authors: P. Nazari, B. Tabone, G. P. Rosotti, M. L. van Gelder, R. Meshaka, E. F. van Dishoeck

Abstract: Some protostellar systems show little or no millimetre line emission of complex organics. This can be interpreted as a low abundance of these molecules, alternatively they could be present in the system but are not seen in the gas. The goal is to investigate the latter hypothesis for methanol. We will attempt to answer the question: Does the presence of a disk and optically thick dust reduce metha… ▽ More Some protostellar systems show little or no millimetre line emission of complex organics. This can be interpreted as a low abundance of these molecules, alternatively they could be present in the system but are not seen in the gas. The goal is to investigate the latter hypothesis for methanol. We will attempt to answer the question: Does the presence of a disk and optically thick dust reduce methanol emission even if methanol is abundant in the ices and gas? Using the radiative transfer code RADMC-3D, methanol emission lines from an envelope-only model and an envelope-plus-disk model are calculated and compared with each other and the observations. Methanol emission from the envelope-only model is always stronger than from the envelope-plus-disk model by at least a factor 2 as long as the disk radius is larger than 30 au (for L=8 L$_{\odot}$). In most cases, this is due to lower temperatures and, hence, the smaller amount of warm methanol inside the snow surface of the envelope-plus-disk model. The intensities drop by more than an order of magnitude for models including high mm opacity dust grains and disk radii of at least 50 au (for L=8 L$_{\odot}$) due to continuum over-subtraction. The line intensities from the envelope-only models overproduce the observations of protostars with lower methanol emission even with large dust optical depth effects. The envelope-plus-disk models can explain the bulk of the observations. However, they can only reproduce the observations of sources with high luminosities and low methanol emission when dust optical depth effects especially continuum over-subtraction in the disk becomes effective. Therefore, both the effects of disk and dust optical depth should be considered to explain the observations. In conclusion: Absence of methanol emission does not imply absence of methanol molecules in either gas or ice. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: Accepted to A&A

Journal ref: A&A 663, A58 (2022)

Characteristics of small protoplanetary disc warps in kinematic observations

Authors: Alison K. Young, Richard Alexander, Giovanni Rosotti, Christophe Pinte

Abstract: Many circumstellar discs appear to have misaligned central regions that give rise to shadows seen in scattered light observations. Small warps ($<20^\circ$ misalignment) are probably more common but are also more difficult to detect than the large misalignments studied previously. We present the characteristics of CO emission that may be used to identify a small disc warp, found from synthetic… ▽ More Many circumstellar discs appear to have misaligned central regions that give rise to shadows seen in scattered light observations. Small warps ($<20^\circ$ misalignment) are probably more common but are also more difficult to detect than the large misalignments studied previously. We present the characteristics of CO emission that may be used to identify a small disc warp, found from synthetic $^{13}$CO maps of a model misaligned circumbinary disc. The spectra are not symmetrical, so fitting a Keplerian model is not appropriate and can hide a warp or lead to spurious features such as spirals appearing in the residuals. We quantify the observed warp structure by fitting sinusoids to concentric annuli of the disc. From this we can trace the radial variation of the peak velocity and of the azimuth of the peak velocity, i.e., the twist. At near face-on inclinations, these radial profiles reveal the warp structure. The twist remains detectable at moderate inclinations (${i_{\rm outer~disc}\lesssim 35^{\circ}}$) in the absence of radial flows but the measured inclination must be accurate to $\lesssim 5^{\circ}$ to allow detection of the radial variation. The observed twist does not provide a direct measure of the warp structure because of its dependence on optical depth. The warp causes broad asymmetries in the channel maps that span several channels and that are distinct from localised features caused by embedded planets and gravitational instability. We suspect that kinematic evidence of warps may have been missed and we suggest a few examples where the data may be revisited. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 16 pages. Accepted to MNRAS

Demographics of young stars and their protoplanetary disks: lessons learned on disk evolution and its connection to planet formation

Authors: C. F. Manara, M. Ansdell, G. P. Rosotti, A. M. Hughes, P. J. Armitage, G. Lodato, J. P. Williams

Abstract: Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys… ▽ More Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys of multiple star-forming regions cover both spectroscopy of young stars and mm interferometry of their protoplanetary disks. This provides an unprecedented statistical sample of stellar masses and mass accretion rates, as well as disk masses and radii, for almost 1000 young stellar objects within 300 pc from us, while also sampling different evolutionary stages, ages, and environments. At the same time, surveys of debris disks are revealing the bulk properties of this class of more evolved objects. This chapter reviews the statistics of these measured global star and disk properties and discusses their constraints on theoretical models describing global disk evolution. Our comparisons of observations to theoretical model predictions extends beyond the traditional viscous evolution framework to include analytical descriptions of magnetic wind effects. Finally, we discuss how recent observational results can provide a framework for models of planet population synthesis and planet formation. △ Less

Submitted 1 March, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Final accepted version after interaction with the referees and community feedback. 21 pages (24 with references), 8 figures. Data table available at http://ppvii.org/chapter/15/

Stellar multiplicity affects the correlation between proto-planetary disc masses and accretion rates: binaries explain high-accretors in Upper Sco

Authors: Francesco Zagaria, Cathie J. Clarke, Giovanni P. Rosotti, Carlo F. Manara

Abstract: In recent years a correlation between mass accretion rates onto new-born stars and their proto-planetary disc masses was detected in nearby young star-forming regions. Although such a correlation can be interpreted as due to viscous-diffusion processes in the disc, highly-accreting sources with low disc masses in more evolved regions remain puzzling. In this paper, we hypothesise that the presence… ▽ More In recent years a correlation between mass accretion rates onto new-born stars and their proto-planetary disc masses was detected in nearby young star-forming regions. Although such a correlation can be interpreted as due to viscous-diffusion processes in the disc, highly-accreting sources with low disc masses in more evolved regions remain puzzling. In this paper, we hypothesise that the presence of a stellar companion truncating the disc can explain these outliers. Firstly, we searched the literature for information on stellar multiplicity in Lupus, Chamaeleon~I and Upper Sco, finding that roughly 20 per cent of the discs involved in the correlation are in binaries or higher-order multiple stellar systems. We prove with high statistical significance that at any disc mass these sources have systematically higher accretion rates than those in single-stars, with the bulk of the binary population being clustered around $M_\mathrm{disc}/\dot{M}_\mathrm{acc}\approx0.1\,\mathrm{Myr}$. We then run coupled gas and dust one-dimensional evolutionary models of tidally truncated discs to be compared with the data. We find that these models are able to reproduce well most of the population of observed discs in Lupus and Upper Sco, even though the unknown eccentricity of each binary prevents an object by object comparison. In the latter region, the agreement improves if the grain coagulation efficiency is reduced, as may be expected in discs around close binaries. Finally, we mention that thermal winds and sub-structures can be important in explaining few outlying sources. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 13 pages, 6 figures; accepted for publication in MNRAS

A large population study of protoplanetary disks: Explaining the millimeter size-luminosity relation with or without sub-structure

Authors: Apostolos Zormpas, Tilman Birnstiel, Giovanni P. Rosotti, Sean M. Andrews

Abstract: Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolut… ▽ More Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $α$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure. △ Less

Submitted 23 May, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: Replaced with the published version in A&A. 19+11 pages, 11+10 figures

Journal ref: A&A 661, A66 (2022)