Search

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

We present a study of the abundance of calcium in the innermost disk of 70 T Tauri stars in the star-forming regions of Chamaeleon I, Lupus and Orion OB1b. We use calcium as a proxy for the refractory material that reaches the inner disk. We used magnetospheric accretion models to analyze the Ca II emission lines and estimate abundances in the accretion flows of the stars, which feed from the inner disks. We find Ca depletion in disks of all three star-forming regions, with 57% of the sample having [Ca/H] < -0.30 relative to the solar abundance. All disks with cavities and/or substructures show depletion, consistent with trapping of refractories in pressure bumps. Significant Ca depletion ([Ca/H] < -0.30) is also measured in 60% of full disks, although some of those disks may have hidden substructures or cavities. We find no correlation between Ca abundance and stellar or disk parameters except for the mass accretion rate onto the star. This could suggest that the inner and outer disks are decoupled, and that the mass accretion rate is related to a mass reservoir in the inner disk, while refractory depletion reflects phenomena in the outer disk related to the presence of structure and forming planets. Our results of refractory depletion and timescales for depletion are qualitatively consistent with expectations of dust growth and radial drift including partitioning of elements and constitute direct evidence that radial drift of solids locked in pebbles takes place in disks., Comment: Accepted for publication in The Astrophysical Journal

Growing evidence suggests that protoplanetary discs may be influenced by late stage infall from the interstellar medium (ISM). It remains unclear the degree to which infall shapes disc populations at ages $\gtrsim 1$~Myr. We explore possible spatial correlations between stellar accretion rates in the Lupus star forming region, which would support the hypothesis that infall can regulate stellar accretion. We consider both the `clustered' stars towards the center of Lupus 3, and the `distributed' stars that are more sparsely distributed across the Lupus complex. We take the observed accretion rates in the literature and explore spatial correlations. In particular, we test whether the clustered stars exhibit a radial gradient in normalised accretion rates, and whether the distributed stars have spatially correlated accretion rates. We find statistically significant correlations for both the clustered and distributed samples. The clustered sample exhibits higher accretion rates in the central region, consistent with the expected Bondi-Hoyle-Lyttleton accretion rate. Stars that are spatially closer among the distributed population also exhibit more similar accretion rates. These results cannot be explained by the stellar mass distribution for either sample. Age gradients are disfavoured, though not discounted, because normalised disc dust masses are not spatially correlated across the region. Spatially correlated stellar accretion rates within the Lupus star forming region argue in favour of an environmental influence on stellar accretion, possibly combined with internal processes in the inner disc. Refined age measurements and searches for evidence of infalling material are potential ways to further test this finding., Comment: 10 pages, 6 figures, accepted for publication in A&A

Classical T Tauri Stars (CTTSs) are highly variable stars that possess gas- and dust-rich disks from which planets form. Much of their variability is driven by mass accretion from the surrounding disk, a process that is still not entirely understood. A multi-epoch optical spectral monitoring campaign of four CTTSs (TW Hya, RU Lup, BP Tau, and GM Aur) was conducted along with contemporaneous HST UV spectra and ground-based photometry in an effort to determine accretion characteristics and gauge variability in this sample. Using an accretion flow model, we find that the magnetospheric truncation radius varies between 2.5-5 R* across all of our observations. There is also significant variability in all emission lines studied, particularly Halpha, Hbeta, and Hgamma. Using previously established relationships between line luminosity and accretion, we find that, on average, most lines reproduce accretion rates consistent with accretion shock modeling of HST spectra to within 0.5 dex. Looking at individual contemporaneous observations, however, these relationships are less accurate, suggesting that variability trends differ from the trends of the population and that these empirical relationships should be used with caution in studies of variability., Comment: 39 pages, 10 figures, 12 tables

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

We have performed a comprehensive study of the Orion Nebula Cluster (ONC) combining the photometric data obtained by the two \textit{HST} Treasury programs that targeted this region. To consistently analyze the rich dataset obtained in a wide variety of filters, we adopted a Bayesian approach to fit the Spectral Energy Distribution of the sources, deriving mass, age, extinction, distance, and accretion for each source in the region. The three dimensional study of mass distribution for bona-fide cluster members shows that mass segregation in the ONC extends to sub-solar masses, while the age distribution strongly supports the idea that star formation in the ONC is best described by a major episode of star formation that happened $\sim 1$ Myr ago. For masses $\gtrsim 0.1$ \Msun, our derived empirical initial mass function (IMF) is in good agreement with a Chabrier system IMF. Both the accretion luminosity (\Lacc) and mass accretion rates (\dMacc) are best described by broken power-law relations. This suggests that for the majority of young circumstellar disks in this cluster the excess emission may be dominated by X-ray-driven photoevaporation by the central star rather than external photoevaporation. If this is the case, the slopes of the power-law relations may be largely determined by the initial conditions set at the onset of the star formation process, which may be quite similar between regions that eventually form clusters of different sizes.

High-angular resolution imaging by ALMA has revealed the near-universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M-dwarfs in the galaxy. Here we present high-resolution (~50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ~10 au for 2M0434, 2M0436 and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared SEDs to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet-disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ~10 au cavities under current observations and theories. Currently, disks around mid-to-late M stars in Taurus show a higher detection frequency of cavities than earlier type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass., Comment: 25 pages, 6 figures. Accepted for publication in ApJ

This paper aims at revisit the physical and dynamical properties of the warm atomic gas in the inner disk region of classical T Tauri stars (CTTs) and relate them to the properties of the outer dusty disk. We used the high resolution (R=115,000) spectra of 36 CTTs observed as part of the GHOsT project and analysed the profile and luminosity of the brightest optical forbidden lines, namely [OI]630 and 557nm, [SII]406 and 673nm, and [NII]658nm. We find that in about 40% of sources the so-called narrow low-velocity component (NLVC) display a peak velocity compatible with the stellar velocity. In these sources, that typically show lower mass accretion rates and the absence of a jet, the [OI]630nm profiles are well fitted by a simple Keplerian disk model, indicating that the emission from the disk is dominant with respect to the wind contribution. For transitional disks (TD), no correlation is found between $R_{kep}$, derived from the line HWHM, and the size of the dust cavity. We also see an anti-correlation between the [OI] 557/630 nm ratio and $R_{kep}$, which suggests that the [OI] emitting region expands as the gas cools and becomes less dense. We confirmed previous findings on the density and temperature ranges implied by the line ratios, and additionally constrained the ionisation fraction in the NLVC to be < 0.1. We however discuss the limits of applying this diagnostic to winds that are not spatially resolved. For the outflow component, we estimated the mass-loss for both the disk winds and jets and compared the results with X-ray photoevaporative models. We conclude that without better knowledge of the wind geometry, and given the limitation of the diagnostics, the mass-loss in the wind traced by the LVC cannot be constrained better than a factor of 100, with a mass-loss/mass-accretion ratio spanning between ~ 0.01 and more than 1., Comment: Accepted for publication in Astronomy and Astrophysics

Protostellar jets and outflows are signposts of active star formation. In H II regions, molecular tracers like CO only reveal embedded portions of the outflow. Outside the natal cloud, outflows are dissociated, ionized, and eventually completely ablated, leaving behind only the high-density jet core. Before this process is complete, there should be a phase where the outflow is partially molecular and partially ionized. In this paper, we capture the HH 900 outflow while this process is in action. New observations from the ERIS/SPIFFIER near-IR integral field unit (IFU) spectrograph using the K-middle filter ($\lambda$=2.06-2.34 $\mu$m) reveal H$_2$ emission from the dissociating outflow and Br-$\gamma$ tracing its ionized skin. Both lines trace the wide-angle outflow morphology but H$_2$ only extends $\sim$5000 au into the H II region while Br-$\gamma$ extends the full length of the outflow ($\sim$12,650 au), indicating rapid dissociation of the molecules. H$_2$ has higher velocities further from the driving source, consistent with a jet-driven outflow. Diagnostic line ratios indicate that photoexcitation, not just shocks, contributes to the excitation in the outflow. We argue that HH 900 is the first clear example of an evaporating molecular outflow and predict that a large column of neutral material that may be detectable with ALMA accompanies the dissociating molecules. Results from this study will help guide the interpretation of near-IR images of externally irradiated jets and outflows such as those obtained with the James Webb Space Telescope (JWST) in high-mass star-forming regions where these conditions may be common., Comment: MNRAS, accepted

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

Accretion plays a central role in the physics that governs the evolution and dispersal of protoplanetary disks. The primary goal of this paper is to analyze the stability over time of the mass accretion rate onto TW Hya, the nearest accreting solar-mass young star. We measure veiling across the optical spectrum in 1169 archival high-resolution spectra of TW Hya, obtained from 1998--2022. The veiling is then converted to accretion rate using 26 flux-calibrated spectra that cover the Balmer jump. The accretion rate measured from the excess continuum has an average of $2.51\times10^{-9}$~M$_\odot$~yr$^{-1}$ and a Gaussian distribution with a FWHM of 0.22 dex. This accretion rate may be underestimated by a factor of up to 1.5 because of uncertainty in the bolometric correction and another factor of 1.7 because of excluding the fraction of accretion energy that escapes in lines, especially Ly$\alpha$. The accretion luminosities are well correlated with He line luminosities but poorly correlated with H$\alpha$ and H$\beta$ luminosity. The accretion rate is always flickering over hours but on longer timescales has been stable over 25 years. This level of variability is consistent with previous measurements for most, but not all, accreting young stars., Comment: Accepted by ApJ. 31 pages

We present VLT/MUSE Narrow Field Mode (NFM) observations of a pair of disc-bearing young stellar objects towards the Orion Bar: 203-504 and 203-506. Both of these discs are subject to external photoevaporation, where winds are launched from their outer regions due to environmental irradiation. Intriguingly, despite having projected separation from one another of only 1.65{\arcsec} (660au at 400pc), 203-504 has a classic teardrop shaped ``proplyd'' morphology pointing towards $\theta^2$Ori A (indicating irradiation by the EUV of that star, rather than $\theta^1$ Ori C) but 203-506 has no ionisation front, indicating it is not irradiated by stellar EUV at all. However, 203-506 does show [CI] 8727{\AA} and [OI] 6300{\AA} in emission, indicating irradiation by stellar FUV. This explicitly demonstrates the importance of FUV irradiation in driving mass loss from discs. We conclude that shielding of 203-506 from EUV is most likely due to its position on the observers side of an ionized layer lying in the foreground of the Huygens Region. We demonstrate that the outflow HH 519, previously thought to be emanating from 203-504 is actually an irradiated cloud edge and identify a new compact outflow from that object approximately along our line of sight with a velocity $\sim130$\,km\,s$^{-1}$., Comment: 14 pages, 8 figures. Accepted for publication in MNRAS

PDS 70 is so far the only young disc where multiple planets have been detected by direct imaging. The disc has a large cavity when seen at sub-mm and NIR wavelengths, which hosts two massive planets. This makes PDS 70 the ideal target to study the physical conditions in a strongly depleted inner disc shaped by two giant planets, and in particular to test whether disc winds can play a significant role in its evolution. Using X-Shooter and HARPS spectra, we detected for the first time the wind-tracing [O I] 6300AA line, and confirm the low-moderate value of mass-accretion rate in the literature. The [O I] line luminosity is high with respect to the accretion luminosity when compared to a large sample of discs with cavities in nearby star-forming regions. The FWHM and blue-shifted peak of the [O I] line suggest an emission in a region very close to the star, favouring a magnetically driven wind as the origin. We also detect wind emission and high variability in the He I 10830AA line, which is unusual for low-accretors. We discuss that, although the cavity of PDS 70 was clearly carved out by the giant planets, the substantial inner disc wind could also have had a significant contribution to clearing the inner-disc., Comment: 13 pages, 7 figures, accepted for publication in ApJ

The Atacama Large Millimeter/submillimeter Array (ALMA) revolutionised our understanding of protoplanetary discs. However, the available data have not given conclusive answers yet on the underlying disc evolution mechanisms (viscosity or MHD winds). Improving upon the current results, mostly based on the analysis of disc sizes, is difficult because larger, deeper and higher angular resolution surveys would be required, which could be prohibitive even for ALMA. In this Letter, we introduce an alternative method to study disc evolution based on $^{12}$CO fluxes. In fact, fluxes can be readily collected using less time-consuming, lower resolution observations, while tracing the same disc physico-chemical processes as sizes: assuming that $^{12}$CO is optically thick, fluxes scale with the disc surface area. We developed a semi-analytical model to compute $^{12}$CO fluxes and benchmarked it against the results of DALI thermochemical models, recovering an agreement within a factor of three. As a proof of concept we compared our models with Lupus and Upper Sco data, taking advantage of the increased samples, by a factor 1.3 (Lupus) and 3.6 (Upper Sco), when studying fluxes instead of sizes. Models and data agree well only if CO depletion is considered. However, the uncertainties on the initial conditions limited our interpretation of the observations. Our new method can be used to design future ad hoc observational strategies to collect better data and give conclusive answers on disc evolution., Comment: 10 pages, 9 figures; accepted for publication in A&A Letters. Code and data available at https://github.com/fzagaria/COpops.git

Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $\sigma$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3\sigma$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5\sigma$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST., Comment: 15 pages, 5 figures. Accepted for publication in ApJ

Absorption features in stellar atmospheres are often used to calibrate photocentric velocities for kinematic analysis of further spectral lines. The Li feature at $\sim$ 6708 {\AA} is commonly used, especially in the case of young stellar objects for which it is one of the strongest absorption lines. However, this is a complex line comprising two isotope fine-structure doublets. We empirically measure the wavelength of this Li feature in a sample of young stars from the PENELLOPE/VLT programme (using X-Shooter, UVES and ESPRESSO data) as well as HARPS data. For 51 targets, we fit 314 individual spectra using the STAR-MELT package, resulting in 241 accurately fitted Li features, given the automated goodness-of-fit threshold. We find the mean air wavelength to be 6707.856 {\AA}, with a standard error of 0.002 {\AA} (0.09 km/s) and a weighted standard deviation of 0.026 {\AA} (1.16 km/s). The observed spread in measured positions spans 0.145 {\AA}, or 6.5 km/s, which is up to a factor of six higher than typically reported velocity errors for high-resolution studies. We also find a correlation between the effective temperature of the star and the wavelength of the central absorption. We discuss how exclusively using this Li feature as a reference for photocentric velocity in young stars could potentially be introducing a systematic positive offset in wavelength to measurements of further spectral lines. If outflow tracing forbidden lines, such as [O i] 6300 {\AA}, are actually more blueshifted than previously thought, this then favours a disk wind as the origin for such emission in young stars., 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

Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation. Aims. In order to systematically study late infall of gas onto disks, we identify candidate Class II sources in close vicinity to a reflection nebula (RN) that may be undergoing this process. Methods. First we targeted Class II sources with known kilo-au scale gas structures - possibly due to late infall of material - and we searched for RNe in their vicinity in optical and near-infrared images. Second, we compiled a catalogue of Class II sources associated with RNe and looked for the large-scale CO structures in archival ALMA data. Using the catalogues of protostars and RNe, we also estimated the probability of Class II sources interacting with surrounding material. Results. All of the sources with large-scale gas structures also exhibit some reflection nebulosity in their vicinity. Similarly, at least five Class II objects associated with a prominent RNe, and for which adequate ALMA observations are available, were found to have spirals or stream-like structures which may be due to late infall. We report the first detection of these structures around S CrA. Conclusions. Our results suggest that a non-negligible fraction of Class II disks in nearby star-forming regions may be associated with RNe and could therefore be undergoing late accretion of gas. Surveys of RNe and kilo-au scale gas structures around Class II sources will allow us to better understand the frequency and impact of late-infall phenomena., Comment: Accepted for publication in Astronomy & Astrophysics

T Tauri stars produce broad Lyman-alpha emission lines that contribute $\sim$88% of the total UV flux incident on the inner circumstellar disks. Lyman-alpha photons are generated at the accretion shocks and in the protostellar chromospheres and must travel through accretion flows, winds and jets, the protoplanetary disks, and the interstellar medium before reaching the observer. This trajectory produces asymmetric, double-peaked features that carry kinematic and opacity signatures of the disk environments. To understand the link between the evolution of Lyman-alpha emission lines and the disks themselves, we model HST-COS spectra from targets included in Data Release 3 of the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) program. We find that resonant scattering in a simple spherical expanding shell is able to reproduce the high velocity emission line wings, providing estimates of the average velocities within the bulk intervening H I. The model velocities are significantly correlated with the K band veiling, indicating a turnover from Lyman-alpha profiles absorbed by outflowing winds to emission lines suppressed by accretion flows as the hot inner disk is depleted. Just 30% of targets in our sample have profiles with red-shifted absorption from accretion flows, many of which have resolved dust gaps. At this stage, Lyman-alpha photons may no longer intersect with disk winds along the path to the observer. Our results point to a significant evolution of Lyman-alpha irradiation within the gas disks over time, which may lead to chemical differences that are observable with ALMA and JWST., Comment: accepted to ApJ

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.

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.

Light curves of young star systems show photometric variability due to different kinematic, and physical processes. One of the main contributors to the photometric variability is the changing mass accretion rate, which regulates the interplay between the forming young star and the protoplanetary disk. We collected high-resolution spectroscopy in eight different epochs, as well as ground-based and space-borne multi-epoch optical and infrared photometry of WX Cha, an M0 binary system, with an almost edge-on disk (i = 87degrees) in the Chamaeleon I star-forming region. Spectroscopic observations cover 72 days, the ground-based optical monitoring covers 42 days while space-borne TESS photometry extends for 56 days. The multi-wavelength light curves exhibit quasi-periodic variability of 0.35 - 0.53 mag in the near-infrared, and of 1.3 mag in g band. We studied the variability of selected emission lines that trace the accretion, computed the accretion luminosity and the mass accretion rate using empirical relations and obtained values of the accretion luminosity between 1.6 and 3.2 Lsun and mass accretion rate between 3.31x10{-7} Msun/yr and 7.76x10^{-7} Msun/yr. Our results show that WX Cha is accreting at a rate larger than what is typical for T Tauri stars in the same star-forming region with the same stellar parameters. We theorize that this is due to the higher disk mass of WX Cha than what is usual for stars with similar stellar mass, and to the binary nature of the system. Daily changes in the accretion luminosity and in the extinction can explain the photometric variability.

The coevolution of T Tauri stars and their surrounding protoplanetary disks dictates the timescales of planet formation. In this paper, we present magnetospheric accretion and inner disk wall model fits to NUV-NIR spectra of nine classical T Tauri stars in Orion OB1b as part of the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey. Using NUV-optical spectra from the Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program and optical-NIR spectra from the PENELLOPE VLT Large Programme, we find that the accretion rates of these targets are relatively high for the region's intermediate age of 5.0 Myr; rates range from $0.5-17.2 \times 10^{-8}$ M$_{\odot}$/yr, with a median value of $1.2\times 10^{-8}$ M$_{\odot}$/yr. The NIR excesses can be fit with 1200-1800 K inner disk walls located at 0.05-0.10 AU from the host stars. We discuss the significance of the choice in extinction law, as the measured accretion rate depends strongly on the adopted extinction value. This analysis will be extended to the complete sample of T Tauri stars being observed through ULLYSES to characterize accretion and inner disks in star-forming regions of different ages and stellar populations., Comment: 22 pages, 3 figures, 8 tables. Accepted for publication in AJ

Classical T Tauri stars are low-mass objects, which are still accreting material from the surrounding circumstellar disk. The accretion process is essential in the formation of Sun-like stars and in setting the properties of the disk at the time when planet formation occurs. We constructed a complex dataset in order to examine the accretion process of VW Cha, a classical T Tauri multiple system with the aim of studying the physical origin of the photometric and spectroscopic variability of the system. The TESS Space Telescope observed VW Cha between 2019 April 22 and June 19, and we complemented these data with contemporaneous ground-based $I_CJHK$ band photometric measurements. In addition, we obtained high-resolution optical spectra with the VLT/ESPRESSO and the 2.2\,m/FEROS instruments. Analyzing these data, we found that the TESS light curve shows photometric variations on timescales from minutes to weeks with a peak-to-peak amplitude of $\sim$0.8 mag. The near-infrared light curves follow the shape of the optical measurements, however, the peak-to-peak amplitudes are slightly increasing towards the longer wavelengths. We took spectra in both fainter and brighter photometric states of the system, allowing us to examine the origin of a photometric brightening event. Our results show that this brightening event can be explained by increased accretion. In addition, our spectroscopic data also suggest that the primary component of VW Cha is a spectroscopic binary, as it was proposed in earlier works., Comment: Accepted for publication in ApJ. 18 pages, 10 figures

Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged), Comment: Submitted to A&A. Currently in revision. 28 pages, 13 figures

While protoplanetary disks are often treated as isolated systems in planet formation models, observations increasingly suggest that vigorous interactions between Class II disks and their environments are not rare. DO Tau is a T Tauri star that has previously been hypothesized to have undergone a close encounter with the HV Tau system. As part of the DESTINYS ESO Large Programme, we present new VLT/SPHERE polarimetric observations of DO Tau and combine them with archival HST scattered light images and ALMA observations of CO isotopologues and CS to map a network of complex structures. The SPHERE and ALMA observations show that the circumstellar disk is connected to arms extending out to several hundred au. HST and ALMA also reveal stream-like structures northeast of DO Tau, some of which are at least several thousand au long. These streams appear not to be gravitationally bound to DO Tau, and comparisons with previous Herschel far-IR observations suggest that the streams are part of a bridge-like structure connecting DO Tau and HV Tau. We also detect a fainter redshifted counterpart to a previously known blueshifted CO outflow. While some of DO Tau's complex structures could be attributed to a recent disk-disk encounter, they might be explained alternatively by interactions with remnant material from the star formation process. These panchromatic observations of DO Tau highlight the need to contextualize the evolution of Class II disks by examining processes occurring over a wide range of size scales., Comment: 42 pages, 20 figures, accepted by ApJ, reduced data available at https://zenodo.org/record/6408903 (typos from v1 fixed)

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

We introduce the STAR-MELT Python package that we developed to facilitate the analysis of time-resolved emission line spectroscopy of young stellar objects. STAR-MELT automatically extracts, identifies and fits emission lines. We summarise our analysis methods that utilises the time domain of high-resolution stellar spectra to investigate variability in the line profiles and corresponding emitting regions. This allows us to probe the innermost disc and accretion structures of YSOs. Local temperatures and densities can be determined using Boltzmann statistics, the Saha equation, and the Sobolev large velocity gradient approximation. STAR-MELT allows for new results to be obtained from archival data, as well as facilitating timely analysis of new data as it is obtained. We present the results of applying STAR-MELT to three YSOs, using spectra from UVES, XSHOOTER, FEROS, HARPS, and ESPaDOnS. We demonstrate what can be achieved for data with disparate time sampling, for stars with different inclinations and variability types. For EX Lupi, we confirm the presence of a localised and stable stellar-surface hot spot associated with the footprint of the accretion column. For GQ Lupi A, we find that the maximum infall rate from an accretion column is correlated with lines produced in the lowest temperatures. For CVSO109 we investigate the rapid temporal variability of a redshifted emission wing, indicative of rotating and infalling material in the inner disc. Our results show that STAR-MELT is a useful tool for such analysis, as well as other applications for emission lines., Comment: Accepted for publication in MNRAS

A large fraction of observed protoplanetary disks in nearby Star-Forming Regions (SFRs) are fainter than expected in CO isotopologue emission. Disks not detected in 13CO line emission are also faint and often unresolved in the continuum emission at an angular resolution of around 0.2 arcseconds. Focusing on the Lupus SFR, the aim of this work is to investigate whether the population of CO-faint disks comprises radially extended and low mass disks - as commonly assumed so far - or if it is of intrinsically radially compact disks, an interpretation that we propose in this paper. The latter scenario was already proposed for individual sources or small samples of disks, while this work targets a large population of disks in a single SFR for which statistical arguments can be made. A new grid of physical-chemical models of compact disks has been run with DALI in order to cover a region of the parameter space that had not been explored before with this code. Such models have been compared with 12CO and 13CO ALMA observations of faint disks in Lupus. Disks that are not detected in 13CO emission and with faint or undetected 12CO emission are consistent with compact disk models. For radially compact disk, the emission of CO isotopologues is mostly optically thick and it scales with the surface area: i.e., it is fainter for smaller objects. The fraction of compact disks is potentially between roughly 50% and 60% of the entire Lupus sample. Deeper observations of 12CO and 13CO at a moderate angular resolution will allow us to distinguish whether faint disks are intrinsically compact, or if they are extended but faint, without the need of resolving them. If the fainter end of the disk population observed by ALMA disk surveys is consistent with such objects being very compact, this will either create a tension with viscous spreading or require MHD winds or external processes to truncate the disks., Comment: 10 pages, 3 main figures, 3 appendix figures. Accepted to A&A

We analysed archival molecular line data of pre-main sequence (PMS) stars in the Lupus and Taurus star-forming regions obtained with ALMA surveys with an integration time of a few minutes per source. We stacked the data of $^{13}$CO and C$^{18}$O (J = 2-1 & 3-2) and CN (N = 3-2, J = 7/2-5/2) lines to enhance the signal-to-noise ratios, and measured the stellar masses of 45 out of 67 PMS stars from the Keplerian rotation in their circumstellar disks. The measured dynamical stellar masses were compared to the stellar masses estimated from the spectroscopic measurements with seven different stellar evolutionary models. We found that the magnetic model of Feiden (2016) provides the best estimate of the stellar masses in the mass range of $0.6~M_{\odot}\leq M_{\star} \leq 1.3~M_{\odot}$ with a deviation of $<$0.7$\sigma$ from the dynamical masses, while all the other models underestimate the stellar masses in this mass range by 20% to 40%. In the mass range of $<0.6~M_{\odot}$, the stellar masses estimated with the magnetic model of Feiden (2016) have a larger deviation ($>2\sigma$) from the dynamical masses, and other, non-magnetic stellar evolutionary models of Siess et al. (2000), Baraffe et al. (2015) and Feiden (2016) show better agreements with the dynamical masses with the deviations of 1.4$\sigma$ to 1.6$\sigma$. Our results show the mass dependence of the accuracy of these stellar evolutionary models., Comment: 31 pages, 10 figures, Accepted by ApJ

We aim to estimate if structures, such as cavities, rings, and gaps, are common in disks around VLMS and to test models of structure formation in these disks. We also aim to compare the radial extent of the gas and dust emission in disks around VLMS, which can give us insight about radial drift. We studied six disks around VLMS in the Taurus star-forming region using ALMA Band 7 ($\sim 340\,$GHz) at a resolution of $\sim0.1''$. The targets were selected because of their high disk dust content in their stellar mass regime. Our observations resolve the disk dust continuum in all disks. In addition, we detect the $^{12}$CO ($J=3-2$) emission line in all targets and $^{13}$CO ($J=3-2$) in five of the six sources. The angular resolution allows the detection of dust substructures in three out of the six disks, which we studied by using UV-modeling. Central cavities are observed in the disks around stars MHO\,6 (M5.0) and CIDA\,1 (M4.5), while we have a tentative detection of a multi-ringed disk around J0433. Single planets of masses $0.1\sim0.4\,M_{\rm{Jup}}$ would be required. The other three disks with no observed structures are the most compact and faintest in our sample. The emission of $^{12}$CO and $^{13}$CO is more extended than the dust continuum emission in all disks of our sample. When using the $^{12}$CO emission to determine the gas disk extension $R_{\rm{gas}}$, the ratio of $R_{\rm{gas}}/R_{\rm{dust}}$ in our sample varies from 2.3 to 6.0, which is consistent with models of radial drift being very efficient around VLMS in the absence of substructures. Our observations do not exclude giant planet formation on the substructures observed. A comparison of the size and luminosity of VLMS disks with their counterparts around higher mass stars shows that they follow a similar relation.

We analyse spatially resolved ALMA observations at 0.9, 1.3, and 3.1 mm for the 26 brightest protoplanetary discs in the Lupus star-forming region. We characterise the discs multi-wavelength brightness profiles by fitting the interferometric visibilities in a homogeneous way, obtaining effective disc sizes at the three wavelengths, spectral index profiles and optical depth estimates. We report three fundamental discoveries: first, the millimeter continuum size - luminosity relation already observed at 0.9 mm is also present at 1.3 mm with an identical slope, and at 3.1 mm with a steeper slope, confirming that emission at longer wavelengths becomes increasingly optically thin. Second, when observed at 3.1 mm the discs appear to be only 9% smaller than when observed at 0.9 mm, in tension with models of dust evolution which predict a starker difference. Third, by forward modelling the sample of measurements with a simple parametric disc model, we find that the presence of large grains ($a_\mathrm{max}>1 $mm) throughout the discs is the most favoured explanation for all discs as it reproduces simultaneously their spectral indices, optical depth, luminosity, and radial extent in the 0.9-1.3 mm wavelength range. We also find that the observations can be alternatively interpreted with the discs being dominated by optically thick, unresolved, substructures made of mm-sized grains with a high scattering albedo., Comment: MNRAS Accepted. 20 pages, 12 figures, 5 tables. Detailed fit results in Appendix E (27 pages, 78 figures), available as supplementary material. Machine-readable version of Table 1, Table 2, and Table 3 available at https://zenodo.org/record/4756381

We present the first ALMA survey of protoplanetary discs at 3 mm, targeting 36 young stellar objects in the Lupus star-forming region with deep observations (sensitivity 20-50 microJy/beam) at ~0.35" resolution (~50 au). Building on previous ALMA surveys at 0.89 and 1.3 mm that observed the complete sample of Class II discs in Lupus at a comparable resolution, we aim to assess the level of grain growth in the relatively young Lupus region. We measure 3 mm integrated fluxes, from which we derive disc-averaged 1-3 mm spectral indices. We find that the mean spectral index of the observed Lupus discs is $\alpha_\mathrm{1-3 mm}=2.23\pm0.06$, in all cases $\alpha_\mathrm{1-3 mm}<3.0$, with a tendency for larger spectral indices in the brightest discs and in transition discs. Furthermore, we find that the distribution of spectral indices in Lupus discs is statistically indistinguishable from that of the Taurus and Ophiuchus star-forming regions. Assuming the emission is optically thin, the low values $\alpha_\mathrm{1-3 mm}\leq 2.5$ measured for most discs can be interpreted with the presence of grains larger than 1 mm. The observations of the faint discs in the sample can be explained without invoking the presence of large grains, namely through a mixture of optically thin and optically thick emission from small grains. However, the bright (and typically large) discs do inescapably require the presence of millimeter-sized grains in order to have realistic masses. Based on a disc mass argument, our results challenge previous claims that the presence of optically thick sub-structures may be a universal explanation for the empirical millimeter size-luminosity correlation observed at 0.89 mm., Comment: MNRAS Accepted; 12 pages, 5 figures, 2 tables. Machine-readable tables available at https://zenodo.org/record/4756282

Protoplanetary disk surveys by the Atacama Large Millimeter/sub-millimeter Array (ALMA) are now probing a range of environmental conditions, from low-mass star-forming regions like Lupus to massive OB clusters like $\sigma$ Orionis. Here we conduct an ALMA survey of protoplanetary disks in $\lambda$ Orionis, a ~5 Myr old OB cluster in Orion, with dust mass sensitivities comparable to the surveys of nearby regions (~0.4 $M_\oplus$). We assess how massive OB stars impact planet formation, in particular from the supernova that may have occurred ~1 Myr ago in the core of $\lambda$ Orionis; studying these effects is important as most planetary systems, including our Solar System, are likely born in cluster environments. We find that the effects of massive stars, in the form of pre-supernova feedback and/or a supernova itself, do not appear to significantly reduce the available planet-forming material otherwise expected at the evolved age of $\lambda$ Orionis. We also compare a lingering massive "outlier" disk in $\lambda$ Orionis to similar systems in other evolved regions, hypothesizing that these outliers host companions in their inner disks that suppress disk dispersal to extend the lifetimes of their outer primordial disks. We conclude with numerous avenues for future work, highlighting how $\lambda$ Orionis still has much to teach us about perhaps one of the most common types of planet-forming environments in the Galaxy., Comment: 21 pages, 8 figures, 3 tables. Accepted to AAS Journals

In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $\kappa_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $\kappa_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks., Comment: 16 pages, 10 figures, accepted for publication in A&A

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with., Comment: Accepted for publication in A&A; 19 pages, 13 figures. Key results shown in Fig. 8 & 10

We report the discovery of a low-mass stellar companion around the young Herbig Be star MWC 297. We performed multi-epoch high-contrast imaging in the near infrared (NIR) with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. The companion is found at projected separation of 244.7$\pm$13.2 au and a position angle of 176.4$\pm$0.1 deg. The large separation supports formation via gravitational instability. From the spectrum, we estimate a mass of 0.1-0.5 M$_{\odot}$, the range conveying uncertainties in the extinction of the companion and in evolutionary models at young ages. The orbit coincides with a gap in the dust disk inferred from the Spectral Energy Distribution (SED). The young age ($\lesssim$ 1 Myr) and mass ratio with the central star ($\sim 0.01$) makes the companion comparable to PDS 70~b, suggesting a relation between formation scenarios and disk dynamics., Comment: 4 figures

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths., Comment: ApJ accepted, 25 pages, 11 figures

The high-quality OmegaCAM photometry of the 3x3 deg around the Orion Nebula Cluster (ONC) in r, and i filters by Beccari et al.(2017) revealed three well-separated pre-main sequences in the color-magnitude diagram (CMD). The objects belonging to the individual sequences are concentrated towards the center of the ONC. The authors concluded that there are two competitive scenarios: a population of unresolved binaries and triples with an exotic mass ratio distribution, or three stellar populations with different ages. We use Gaia DR2 in combination with the photometric OmegaCAM catalog to test and confirm the presence of the putative three stellar populations. We also study multiple stellar systems in the ONC for the first time using Gaia DR2. We confirm that the second and third sequence members are more centrally concentrated towards the center of the ONC. In addition we find an indication that the parallax and proper motion distributions are different among the members of the stellar sequences. The age difference among stellar populations is estimated to be 1-2 Myr. We use Gaia measurements to identify and remove as many unresolved multiple system candidates as possible. Nevertheless we are still able to recover two well-separated sequences with evidence for the third one, supporting the existence of the three stellar populations. We were able to identify a substantial number of wide binary objects (separation between 1000-3000 au). This challenges previously inferred values that suggested no wide binary stars exist in the ONC. Our inferred wide-binary fraction is approx 5%. We confirm the three populations correspond to three separated episodes of star formation. Based on this result, we conclude that star formation is not happening in a single burst in this region. (abridged), Comment: Astronomy and Astrophysics (A&A) accepted. 12 pages, 9 figures + appendix. New version with language corrections and new ID values in Tab.A.1

We describe EPIC 205718330 and EPIC 235240266, two systems identified in the K2 data whose light curves contain episodic drops in brightness with shapes and durations similar to those of the young "dipper" stars, yet shallower by ~1-2 orders of magnitude. These "little dippers" have diverse profile shapes with durations of ~0.5-1.0 days and depths of ~0.1-1.0% in flux; however, unlike most of the young dipper stars, these do not exhibit any detectable infrared excess indicative of protoplanetary disks, and our ground-based follow-up spectra lack any signatures of youth while indicating these objects as kinematically old. After ruling out instrumental and/or data processing artifacts as sources of the dimming events, we investigate possible astrophysical mechanisms based on the light curve and stellar properties. We argue that the little dippers are consistent with transits of star-grazing exocomets, and speculate that they are signposts of massive non-transiting exoplanets driving the close-approach orbits., Comment: 14 pages, 9 Figures; accepted to MNRAS

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet., Comment: Accepted for publication in A&A, 11 pages, 8 figures

Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases., Comment: 23 pages, 11 figures, ApJ accepted

ALMA surveys of nearby star-forming regions have shown that the dust mass in the disk is correlated with the stellar mass, but with a large scatter. This scatter could indicate either different evolutionary paths of disks or different initial conditions within a single cluster. We present ALMA Cycle 3 follow-up observations for 14 Class II disks that were low S/N detections or non-detections in our Cycle 2 survey of the $\sim 2$ Myr-old Chamaeleon I star-forming region. With 5 times better sensitivity, we detect millimeter dust continuum emission from six more sources and increase the detection rate to 94\% (51/54) for Chamaeleon I disks around stars earlier than M3. The stellar-disk mass scaling relation reported in \citet{pascucci2016} is confirmed with these updated measurements. Faint outliers in the $F_{mm}$--$M_*$ plane include three non-detections (CHXR71, CHXR30A, and T54) with dust mass upper limits of 0.2 M$_\oplus$ and three very faint disks (CHXR20, ISO91, and T51) with dust masses $\sim 0.5$ M$_\oplus$. By investigating the SED morphology, accretion property and stellar multiplicity, we suggest for the three millimeter non-detections that tidal interaction by a close companion ($<$100 AU) and internal photoevaporation may play a role in hastening the overall disk evolution. The presence of a disk around only the secondary star in a binary system may explain the observed stellar SEDs and low disk masses for some systems., Comment: ApJ accepted

Context. Transition disks (TDs) are circumstellar disks with inner regions highly depleted in dust. TDs are observed in a small fraction of disk-bearing objects at ages of 1-10 Myr. They are important laboratories to study evolutionary effects in disks, from photoevaporation to planet-disk interactions. Aims. We report the discovery of a large inner dust-empty region in the disk around the very low mass star CIDA 1 (M$_{\star} \sim 0.1-0.2$ M$_{\odot}$). Methods. We used ALMA continuum observations at 887$\mu$m, which provide a spatial resolution of $0."21\times0."12$ ($\sim$15$\times$8 au in radius at 140 pc). Results. The data show a dusty ring with a clear cavity of radius $\sim$20 au, the typical characteristic of a TD. The emission in the ring is well described by a narrow Gaussian profile. The dust mass in the disk is $\sim$17 M$_{\oplus}$. CIDA 1 is one of the lowest mass stars with a clearly detected millimeter cavity. When compared to objects of similar stellar mass, it has a relatively massive dusty disk (less than $\sim5$% of Taurus Class II disks in Taurus have a ratio of $M_{\rm{disk}}/M_{\star}$ larger than CIDA 1) and a very high mass accretion rate (CIDA 1 is a disk with one of the lowest values of $M_{\rm{disk}}/\dot M$ ever observed). In light of these unusual parameters, we discuss a number of possible mechanisms that can be responsible for the formation of the dust cavity (e.g., photoevaporation, dead zones, embedded planets, close binary). We find that an embedded planet of a Saturn mass or a close binary are the most likely possibilities., Comment: Accepted for publication in A&A, minor correction in one equation

In recent ALMA surveys, the gas distributions and velocity structures of most of the protoplanetary disks can still not be imaged at high S/N due to the short integration time. In this work, we re-analyzed the ALMA 13CO (3-2) and C18O (3-2) data of 88 young stellar objects in Lupus with the velocity-aligned stacking method to enhance S/N and to study the kinematics and disk properties traced by molecular lines. This method aligns spectra at different positions in a disk based on the projected Keplerian velocities at their positions and then stacks them. This method enhances the S/N ratios of molecular-line data and allows us to obtain better detections and to constrain dynamical stellar masses and disk orientations. We obtain 13CO detections in 41 disks and C18O detections in 18 disks with 11 new detections in 13CO and 9 new detections in C18O after applying the method. We estimate the disk orientations and the dynamical stellar masses from the 13CO data. Our estimated dynamical stellar masses correlate with the spectroscopic stellar masses, and in a subsample of 16 sources, where the inclination angles are better constrained, the two masses are in a good agreement within the uncertainties and with a mean difference of 0.15 Msun. With more detections of fainter disks, our results show that high gas masses derived from the 13CO and C18O lines tend to be associated with high dust masses estimated from the continuum emission. Nevertheless, the scatter is large (0.9 dex), implying large uncertainties in deriving the disk gas mass from the line fluxes. We find that with such large uncertainties it is expected that there is no correlation between the disk gas mass and the mass accretion rate with the current data. Deeper observations to detect disks with gas masses <1E-5 Msun in molecular lines are needed to investigate the correlation between the disk gas mass and the mass accretion rate., Comment: Submitted to A&A

Transition disks with large dust cavities around young stars are promising targets for studying planet formation. Previous studies have revealed the presence of gas cavities inside the dust cavities hinting at recently formed, giant planets. However, many of these studies are biased towards the brightest disks in the nearby star forming regions, and it is not possible to derive reliable statistics that can be compared with exoplanet populations. We present the analysis of 11 transition disks with large cavities (>20 AU radius) from a complete disk survey of the Lupus star forming region, using ALMA Band 7 observations at 0.3" (22-30 AU radius) resolution of the 345 GHz continuum, 13CO and C18O 3-2 observations and the Spectral Energy Distribution of each source. Gas and dust surface density profiles are derived using the physical-chemical modeling code DALI. This is the first study of transition disks of large cavities within a complete disk survey within a star forming region. The dust cavity sizes range from 20-90 AU radius and in three cases, a gas cavity is resolved as well. The deep drops in gas density and large dust cavity sizes are consistent with clearing by giant planets. The fraction of transition disks with large cavities in Lupus is ~11%, which is inconsistent with exoplanet population studies of giant planets at wide orbits. Furthermore, we present a hypothesis of an evolutionary path for large massive disks evolving into transition disks with large cavities., Comment: 29 pages, 15 figures, Accepted by ApJ

The early stages of the pre-Main-Sequence (PMS) evolution of young stars are characterized by a substantial interaction between the central star and the surrounding disk. Among the processes that take place, the accretion of material from the disk onto the star is a crucial one, as it traces the disk evolution and, in turn, affects how planets are formed and evolve in disks. Even if accretion has been studied for many years, it is still the focus of a large body of work, as we do not fully understand how it happens., Comment: Perspective article for the Star Formation Newsletter 300. 6 pages, 4 figures. Also available at http://www.ifa.hawaii.edu/~reipurth/newsletter/newsletter300.pdf

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year., Comment: 23 pages, 17 figures, accepted to A&A

The determination of the distance to dark star-forming clouds is a key parameter to derive the properties of the cloud itself, and of its stellar content. This parameter is still loosely constrained even in nearby star-forming regions. We want to determine the distances to the clouds in the Chamaeleon-Musca complex and to explore the connection between these clouds and the large scale cloud structures in the galaxy. We use the newly estimated distances obtained from the parallaxes measured by the Gaia satellite and included in the Tycho-Gaia Astrometric Solution catalog. When known members of a region are included in this catalog we use their parallaxes to infer the distance to the cloud. Otherwise, we analyze the dependence of the color excess on the distance of the stars and look for a turn-on of this excess, which is a proxy of the position of the front-edge of the star-forming cloud. We are able to measure the distance to the three Chamaeleon clouds. The distance to Chamaeleon I is 179 pc, 20 pc further away than previously assumed. The Chamaeleon II cloud is located at the distance of 181 pc, which agrees with previous estimates. We are able to measure for the first time a distance to the Chamaeleon III cloud of 199 pc. Finally, the distance of the Musca cloud is smaller than 603 pc. These estimates do not allow us to distinguish between the possibility that the Chamaeleon clouds are part of a sheet of clouds parallel to the galactic plane, or perpendicular to it. Gaia Data Release 2 will allow us to put more stringent constraints on the distances to these clouds by giving us access to parallax measurements for a larger number of members of these regions., Comment: Accepted for publication on A&A. Abstract shortened for arxiv constraints