Search

There is still much uncertainty around the mechanism that rules the accretion of proto-planetary disks. In the last years, Magnetohydrondynamic (MHD) wind-driven accretion has been proposed as a valid alternative to the more conventional viscous accretion. In particular, winds have been shown to reproduce the observed correlation between the mass of the disk Md and the mass accretion rate onto the central star Macc, but this has been done only for specific conditions. It is not clear whether this implies fine tuning or if it is a general result. We investigate under which conditions the observed correlation between the mass of the disk Md and the mass accretion rate onto the central star Macc can be obtained. We find that, in the absence of a correlation between the initial mass M0 and the initial accretion timescale tacc,0, the slope of the Md-Macc correlation depends on the value of the spread of the initial conditions of masses and lifetimes of disks. Then, we clarify the conditions under which a disk population can be fitted with a single power-law. Moreover, we derive an analytical expression for the spread of log(Md/Macc) valid when the spread of tacc is taken to be constant. In the presence of a correlation between M0 and tacc,0, we derive an analytical expression for the slope of the Md-Macc correlation in the initial conditions of disks and at late times. We conclude that MHD winds can predict the observed values of the slope and the spread of the Md-Macc correlation under a broad range of initial conditions. This is a fundamental expansion of previous works on the MHD paradigm, exploring the establishment of this fundamental correlation beyond specific initial conditions., Comment: Accepted for publication in A&A. 13 pages, 16 figures

The evolution of protoplanetary disks in regions with massive OB stars is influenced by externally driven winds that deplete the outer parts of disks. These winds have previously been studied via forbidden oxygen emission lines, which also arise in isolated disks in low-mass star forming-regions (SFRs) with weak external UV fields in photoevaporative or magnetic (internal) disk winds. It is crucial to determine how to disentangle external winds from internal ones. Here, we report a proxy for unambiguously identifying externally driven winds with a forbidden line of neutral atomic carbon, [C i] 8727 A. We compare for the first time the spatial location of the emission in the [O i] 5577 A, [O i] 6300 A, and [C i] 8727 A lines traced by VLT/MUSE-NFM, with the ALMA Band 7 continuum disk emission in a sample of 12 proplyds in the Orion Nebula Cluster (ONC). We confirm that the [O i] 5577 A emission is co-spatial with the disk emission, whereas the [O i] 6300 A is emitted both on the disk surface and on the ionization front of the proplyds. We show for the first time that the [C i] 8727 A line is also co-spatial with the disk surface in proplyds, as seen in the MUSE and ALMA data comparison. To verify whether the [C i] 8727 A line is detected in regions where external photoevaporation is not expected, we examine VLT/X-Shooter spectra for young stars in low-mass SFRs. Although the [O i] lines are well detected in all these targets, there is <<10% detection rate in the case of the [C i] 8727 A line. This number increases substantially to a ~40% detection rate in sigma-Orionis, a region with intermediate UV radiation. The spatial location of the [C i] 8727 A line emission and the lack of its detection in isolated disks in low-mass SFRs strongly suggest that this line is a tell-tale tracer of externally driven photoevaporative winds, which agrees with recent excitation models., Comment: Accepted to A&A on October 25, 2024

The mid-infrared water vapor emission spectrum provides a novel way to characterize the delivery of icy pebbles towards the innermost ($<5$ au) regions of planet-forming disks. Recently, JWST MIRI-MRS showed that compact disks exhibit an excess of low-energy water vapor emission relative to extended multi-gapped disks, suggesting that icy pebble drift is more efficient in the former. We carry out detailed emission line modeling to retrieve the excitation conditions of rotational water vapor emission in a sample of four compact and three extended disks within the JDISC Survey. We present two-temperature H$_2$O slab model retrievals and, for the first time, constrain the spatial distribution of water vapor by fitting parametric radial temperature and column density profiles. Such models statistically outperform the two-temperature slab fits. We find a correlation between the observable hot water vapor mass and stellar mass accretion rate, as well as an anti-correlation between cold water vapor mass and sub-mm dust disk radius, confirming previously reported water line flux trends. We find that the mid-IR spectrum traces H$_2$O with temperatures down to 180-300 K, but the coldest 150-170 K gas remains undetected. Furthermore the H$_2$O temperature profiles are generally steeper and cooler than the expected `super-heated' dust temperature in passive irradiated disks. The column density profiles are used to estimate icy pebble mass fluxes, which suggest that compact and extended disks may produce markedly distinct inner-disk exoplanet populations if local feeding mechanisms dominate their assembly., Comment: Accepted for publication in ApJ

(Abridged) Observational surveys of entire star-forming regions have provided evidence of power-law correlations between the disc properties and the stellar mass, especially the disc mass (${M_d \propto M_*}^{\lambda_m}$) and the accretion rate ($\dot M \propto {M_*}^{\lambda_{acc}}$). Whether the secular disc evolution affects said correlations is still debated: while the purely viscous scenario has been probed, other mechanisms could impact differently. We study the evolution of the slopes $\lambda_m$ and $\lambda_{acc}$ in the wind-driven and hybrid case and compare it to the viscous prediction, using a combination of analytical calculations and numerical simulations (performed with the 1D population synthesis code Diskpop, that we also present and release). Assuming $M_d(0) \propto {M_*}^{\lambda_{m, 0}}$ and $\dot M(0) \propto {M_*}^{\lambda_{acc, 0}}$ as initial conditions, we find that viscous and hybrid accretion preserve the shape of the correlations and evolve their slope; on the other hand, MHD winds change the shape of the correlations, bending them according to the scaling of the accretion timescale with the stellar mass. We also show how a spread in the initial conditions conceals this behaviour. We then analyse the impact of disc dispersal, and find that the currently available sample sizes ($\sim 30$ discs at 5 Myr) introduce stochastic oscillations in the slopes evolution, which dominate over the physical signatures. Increasing the sample size could mitigate this issue: $\sim 140$ discs at 5 Myr, corresponding to the complete Upper Sco sample, would give small enough error bars to use the evolution of the slopes as a proxy for the driving mechanism of disc evolution. Finally, we discuss how the observational claim of steepening slopes necessarily leads to an initially steeper $M_d - M_*$ correlation with respect to $\dot M - M_*$., Comment: 14 pages, 5 figures and 3 tables (plus Appendix). Accepted for publication in A&A. Includes Diskpop and popcorn software release

Detecting planet signatures in protoplanetary disks is fundamental to understanding how and where planets form. In this work, we report dust and gas observational hints of planet formation in the disk around 2MASS-J16120668-301027, as part of the ALMA Large Program "AGE-PRO: ALMA survey of Gas Evolution in Protoplanetary disks". The disk was imaged with the Atacama Large Millimeter/submillimeter Array (ALMA) at Band 6 (1.3 mm) in dust continuum emission and four molecular lines: $^{12}$CO(J=2-1), $^{13}$CO(J=2-1), C$^{18}$O(J=2-1), and H$_2$CO(J=3$_{(3,0)}$-2$_{(2,0)}$). Resolved observations of the dust continuum emission (angular resolution of $\sim 150$ mas, 20 au) show a ring-like structure with a peak at $0.57 ^{\prime \prime}$ (75 au), a deep gap with a minimum at 0.24$^{\prime \prime}$ (31 au), an inner disk, a bridge connecting the inner disk and the outer ring, along with a spiral arm structure, and a tentative detection (to $3\sigma$) of a compact emission at the center of the disk gap, with an estimated dust mass of $\sim 2.7-12.9$ Lunar masses. We also detected a kinematic kink (not coincident with any dust substructure) through several $^{12}$CO channel maps (angular resolution $\sim$ 200 mas, 30 au), located at a radius of $\sim 0.875^{\prime \prime}$ (115.6 au). After modeling the $^{12}$CO velocity rotation around the protostar, we identified a tentative rotating-like structure at the kink location with a geometry similar to that of the disk. We discuss potential explanations for the dust and gas substructures observed in the disk, and their potential connection to signatures of planet formation., Comment: 24 pages, 15 figures

Whether the angular momentum of protoplanetary discs is redistributed by viscosity or extracted by magnetised winds is a long-standing question. Demographic indicators, such as gas disc sizes and stellar accretion rates, have been proposed as ways of distinguishing between these two mechanisms. In this paper, we implement one-dimensional gas simulations to study the evolution of "hybrid" protoplanetary discs simultaneously driven by viscosity and magnetised winds, with dead zones present. We explore how the variations of disc properties, including initial disc sizes, dead zone sizes and angular momentum transport efficiency, affect stellar accretion rates, disc surface density profiles, disc sizes, disc lifetimes, and cumulative mass loss by different processes. Our models show that the expansion of the gas disc size can be sustained when the majority of angular momentum is removed by the magnetised wind for individual protoplanetary discs. However, when we can only observe discs via demographic screenshots, the variation of disc sizes with time is possibly diminished by the disc "personalities", by which we mean the variations of initial disc properties among different discs. Our "hybrid" models re-assess association of the two demographic indicators with mechanisms responsible for angular momentum transport and suggest additional diagnostics are required to assist the differentiation., Comment: 23 pages, 17 figures. Accepted for publication in MNRAS

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., Comment: Accepted by ApJ, 10 pages, 6 figures

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.

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\alpha\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.

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., Comment: Submitted to A&A on December 18, 2023. Accepted on April 2, 2024

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{\deg}, 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.

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., Comment: 16 pages, 12 figures, accepted for publication on A&A

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., Comment: Accepted for publication in Astronomy & Astrophysics, 32 pages, 28 figures

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., Comment: 13 pages, 4 figures; accepted for publication on ApJL

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., Comment: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

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., Comment: Accepted for publication in A&A, 27 pages, 28 figures

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., Comment: 11 pages, 6 figures, 1 table. Published in ApJL

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., Comment: Accepted for publication in A&A; 13 pages, 13 figures, 1 table

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., Comment: 23 pages, 12 figures, accepted in ApJ

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., Comment: Accepted for publication in ApJ Letters

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 ($\alpha$ $\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., Comment: 24 Pages, 15 Figures, Accepted to ApJ

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., Comment: 11 pages, 6 figures, accepted by MNRAS. Welcome any comments and suggestions!

[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., Comment: Accepted for publication in A&A

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., Comment: 19 pages, 6 figures, accpeted for publication in A&A

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., Comment: 24 pages, 7 figures. Accepted for publication on New Astronomy Reviews; small changes following community feedback in v2

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., Comment: Accepted for publication in EPJ Plus Focus Point on Environmental and Multiplicity Effects on Planet Formation G. Lodato and C.F. Manara (Guest editors)

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.

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., Comment: 5 pages, 3 figures

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., Comment: 16 pages, 9 figures, accepted for publication in A&A

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.

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 $\Lambda$ 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 $\Lambda$ translate into a turbulence level of $\alpha_{\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 $\alpha_{\rm turb}\,{\sim}1.2\,{\times}\,10^{-2}$. Due to the degeneracy between $\Lambda$ and the depth of dust surface density drops, the turbulence strength in the D86 gap region is not constrained., Comment: 13 pages, 7 figures, accepted for publication in Science China Physics, Mechanics & Astronomy

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., Comment: Accepted to ApJL, 14 pages, 10 figures

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$., Comment: Accepted for publication in MNRAS - 28 pages, 17 figures

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 $\alpha_{\rm dw}$, which is an analogue to the turbulent viscosity $\alpha_{\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 $\alpha_{\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 $\alpha_{\rm v}$ and $\alpha_{\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., Comment: 14 pages, 9 figures, accepted to MNRAS

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., Comment: Accepted to MNRAS

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}}^{\lambda_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{\lambda_{\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, $\lambda_{\mathrm{m}}$ = $\lambda_{\mathrm{acc}}$, and that this limit is uniquely determined by the initial slopes $\lambda_{\mathrm{m}, 0}$ and $\lambda_{\mathrm{acc}, 0}$. We recover the increasing trend claimed from the observations when the difference in the initial values, $\delta_0 = \lambda_{\mathrm{m}, 0} - \lambda_{\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 $\lambda_{\mathrm{m}, 0} \in [1.2, 2.1]$, $\lambda_{\mathrm{acc}, 0} \in [0.7, 1.5]$., Comment: Accepted for publication in MNRAS. 14 pages, 8 figures

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 $\alpha\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., Comment: 19 pages, 10 figures; accepted for publication in MNRAS

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 $\mu$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., Comment: 6 pages, 4 figures, accepted for publication in MNRAS Letters

High spatial resolution CO observations of mid-inclination (30-75{\deg}) 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}}$., Comment: 31 pages, 17 figures, accepted for publication in ApJ. Image cubes available at https://zenodo.org/record/6410045

Resolved dust continuum and CO line ALMA imaging, and in some cases detection of H$\alpha$ 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., Comment: 17 pages, 10 figures, Accepted by MNRAS

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 $\alpha$ 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., Comment: Accepted for publication in ApJ

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., Comment: 16 pages. Accepted to MNRAS

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., Comment: 13 pages, 6 figures; accepted for publication in MNRAS

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 $\alpha$ 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., Comment: Replaced with the published version in A&A. 19+11 pages, 11+10 figures

It is still unclear whether the evolution of protoplanetary disks, a key ingredient in the theory of planet formation, is driven by viscous turbulence or magnetic disk winds. As viscously evolving disks expand outward over time, the evolution of disk sizes is a discriminant test for studying disk evolution. However, it is unclear how the observed disk size changes over time if disk evolution is driven by magnetic disk winds. Combining the thermochemical code DALI with the analytical wind-driven disk evolution model presented in Tabone et al. (2021a), we study the time evolution of the observed gas outer radius as measured from CO rotational emission ($R_{\rm CO, 90\%}$). The evolution of $R_{\rm CO, 90\%}$ is driven by the evolution of the disk mass, as the physical radius stays constant over time. For a constant $\alpha_{\rm DW}$, an extension of the $\alpha-$Shakura-Sunyaev parameter to wind-driven accretion, $R_{\rm CO, 90\%}$ decreases linearly with time. Its initial size is set by the disk mass and the characteristic radius $R_c$, but only $R_c$ affects the evolution of $R_{\rm CO, 90\%}$, with a larger $R_c$ resulting in a steeper decrease of $R_{\rm CO, 90\%}$. For a time-dependent $\alpha_{\rm DW}$ $R_{\rm CO, 90\%}$ stays approximately constant during most of the disk lifetime until $R_{\rm CO, 90\%}$ rapidly shrinks as the disk dissipates. The constant $\alpha_{\rm DW}$-models are able to reproduce the observed gas disk sizes in the $\sim1-3$ Lupus and $\sim5-11$ Myr old Upper Sco star-forming regions. However, they likely overpredict the gas disk size of younger $(\lessapprox0.7\ \mathrm{Myr})$ disks., Comment: 22 pages, 16 figures, accepted in ApJ

The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solutions of disc evolution and a disc population synthesis approach. We show that MHD disc winds alone can account for both disc dispersal and accretion properties. The decline of disc fraction over time is reproduced by assuming that the initial accretion timescale (a generalization of the viscous timescale) varies from disc to disc and that the decline of the magnetic field strength is slower than that of the gas. The correlation between accretion rate and disc mass, and the dispersion of the data around the mean trend as observed in Lupus is then naturally reproduced. The model also accounts for the rapidity of the disc dispersal. This paves the way for planet formation models in the paradigm of wind-driven accretion., Comment: Accepted for publication in MNRAS Letters

The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscous evolution for the MHD wind case. The transport of angular momentum and mass induced by the wind is parameterized by an $\alpha$-like parameter and by the magnetic lever arm parameter $\lambda$. Extensions of the paradigmatic Lynden-Bell and Pringle similarity solutions to the wind case are presented. We show that wind-driven accretion leads to a steeper decrease in the disc mass and accretion rate than in viscous models due to the absence of disc spreading. If the decline of the magnetic field strength is slower than that of the gas surface density, the disc is dispersed after a finite time. The evolution of the disc in the $\dot{M}_*-M_D$ plane is sensitive to the wind and turbulence parameters. A disc population evolving under the action of winds can exhibit a correlation between $\dot{M}_*$ and $M_D$ depending on the initial conditions. The simplified framework proposed in this work opens to a new avenue to test the effectiveness of wind-driven accretion from the observed disc demographics and constitutes an important step to include wind-driven accretion in planet population synthesis models., Comment: Accepted for publication in MNRAS

The inspiral of supermassive black-hole binaries in gas-rich environment is driven by the presence of an accretion disc and viscous interactions tend to align the spin of the black holes with the orbital angular momentum of the disc. Recent work introduced a new iterative approach to describe the alignment process and the resulting non-linear evolution of the surrounding warped accretion disc. Their model predicted that black-hole spins reach either full alignment or a critical obliquity where solutions to the warp equations cease to exist. In this paper, we show that this critical region corresponds to the disc breaking phenomenon, where the disc is disrupted into two or more discrete sections. We use 3D hydrodynamical simulations to (i) recover the predictions of the semi-analytic model and (ii) unveil a richer phenomenology where the disc exhibits either unsuccessful, single and multiple breaks. We additionally identify hydrodynamic effects such as spiral arms that are able to stabilise the disc against breaking beyond criticality. Our results show that when disc breaking occurs, the ability of black holes and disc to align is compromised and in some cases even prevented as the binary inspirals., Comment: Accepted for publication in MNRAS. 15 pages, 11 figures. Movies of the simulations in this paper can be found at https://youtube.com/playlist?list=PLVCh7-wtOgbJ_NOMDJ3JYOgEAAr0LsJNy

We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at $R=94\pm6$ au, $\phi=50\pm3^\circ$ (P94), is near the centre of one of the gaps in dust continuum emission, and is consistent with a planet mass of 1 $\rm M_{\rm Jup}$. The other planet, located at $R=261\pm4$ au, $\phi=57\pm1^\circ$ (P261), is in the region where a velocity kink was previously observed in $^{12}$CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disc, and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles we detect gas gaps at $R=38$ au, $R=88$ au, and $R=136$ au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at $R=88$ au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the DISCMINER is capable of finding localised, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary discs., Comment: 22 pages, 14 figures, 1 table. Accepted for publication in ApJ

Stars in globular clusters formed and evolved in the most extreme environment: high density and low metallicity. If the formation of stars and planets are at all sensitive to environmental conditions, this should therefore be evident in globular clusters. Observations have indicated that hot Jupiters are at least an order of magnitude less prevalent in the central region of the globular cluster 47 Tucanae than in the field. In this work, we explore the claims in the literature for additional consequences for the low mass stellar initial mass function. Tidal capture, the mechanism that produces X-ray binaries in globular clusters, applies also to brown dwarfs (BDs). This process produces tight stellar-BD binaries that would be detectable by transit surveys. Applying a Monte Carlo dynamical evolution model, we compute the overall BD capture rates. We find that the number of captures is lower than previous estimates. Capture efficiency increases steeply with stellar mass, which means that mass segregation reduces capture efficiency as BDs and low mass stars occupy the same regions. The result of this effect is that the current constraints on the short period companion fraction remains marginally consistent with initially equal numbers of BDs and stars. However, our findings suggest that expanding the sample in 47 Tuc or surveying other globular clusters for close sub-stellar companions can yield constraints on the sub-stellar initial mass function in these environments. We estimate the capture rates in other globular clusters and suggest that 47 Tuc remains a promising target for future transit surveys., Comment: 15 pages, 13 figures, accepted for publication in MNRAS