To understand the formation of planetary systems, one needs to understand the initial conditions of planet formation, i.e. the young gas-rich planet forming disks. Spatially resolved high-contrast observations are of particular interest, since substructures in disks, linked to planet formation, can be detected and close companions or even planets in formation embedded in the disk can be revealed. In this study we present the first result of the DESTINYS survey (Disk Evolution Study Through Imaging of Nearby Young Stars). DESTINYS is an ESO/SPHERE large program that aims at studying disk evolution in scattered light, mainly focusing on a sample of low-mass stars (<1$M_\odot$) in nearby (~200 pc) star-forming regions. In this particular study we present the observations of the ET Cha (RECX 15) system, a nearby 'old' classical T Tauri star (5-8 Myr, ~100 pc), which is still strongly accreting. We use SPHERE/IRDIS in H-band polarimetric imaging mode to obtain high contrast images of the ET Cha system to search for scattered light from the circumstellar disk as well as thermal emission from close companions. We additionally employ VLT/NACO total intensity archival data taken in 2003. We report here the discovery of a low-mass (sub)stellar companion with SPHERE/IRDIS to ET Cha. We are estimating the mass of this new companion based on photometry. Depending on the system age it is a 5 Myr, 50 $M_{Jup}$ brown dwarf or an 8 Myr, 0.10 $M_\odot$ M-type pre-main-sequence star. We explore possible orbital solutions and discuss the recent dynamic history of the system. Independent of the precise companion mass we find that the presence of the companion likely explains the small size of the disk around ET Cha. The small separation of the binary pair indicates that the disk around the primary component is likely clearing from the outside in, explaining the high accretion rate of the system.

We report on the discovery of periodic dips in the multiband lightcurve of ISO-ChaI 52, a young stellar object in the Chamaeleon I dark cloud. This is one among the peculiar objects that display very low or negligible accretion both in their UV continuum and spectral lines, although they present a remarkable infrared excess emission characteristic of optically-thick circumstellar disks. We have analyzed a VLT/X-Shooter spectrum with the tool ROTFIT to determine the stellar parameters. The latter, along with photometry from our campaign with the REM telescope and from the literature, have allowed us to model the spectral energy distribution and to estimate the size and temperature of the inner and outer disk. From the rotational period of the star/disk system of 3.45 days we estimate a disk inclination of 36$^\circ$. The depth of the dips in different bands has been used to gain information about the occulting material. A single extinction law is not able to fit the observed behavior, while a two-component model of a disk warp composed of a dense region with a gray extinction and an upper layer with an ISM-type extinction provides a better fit of the data.

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.

Determining the mechanisms that drive the evolution of protoplanetary disks is a necessary step to understand how planets form. Here we measured the mass accretion rate for young stellar objects at age >5 Myr, a critical test for the current models of disk evolution. We present the analysis of the spectra of 36 targets in the ~5-10 Myr old Upper Scorpius region for which disk masses were measured with ALMA. We find that the mass accretion rates in this sample of old but still survived disks are similarly high as those of the younger (<3 Myr old) star-forming regions of Lupus and Cha I, when considering the dependence on stellar and disk mass. In particular, several disks show high mass accretion rates >10^-9 Msun/yr while having low disk masses. Furthermore, the median values of the measured mass accretion rates in the disk mass ranges where our sample is complete at a level ~60-80% are compatible in these three regions. At the same time, the spread of mass accretion rates at any given disk mass is still >0.9 dex even at age>5 Myr. These results are in contrast with simple models of viscous evolution, which would predict that the values of the mass accretion rate diminish with time, and a tighter correlation with disk mass at age>5 Myr. Similarly, simple models of internal photoevaporation cannot reproduce the observed mass accretion rates, while external photoevaporation might explain the low disk masses and high accretion rates. A partial possible solution to the discrepancy with the viscous models is that the gas-to-dust ratio of the disks at >5 Myr is significantly different and higher than the canonical 100, as suggested by some dust and gas disk evolution models. The results shown here require the presence of several inter-playing processes, such as detailed dust evolution, external photoevaporation and possibly MHD winds, to explain the secular evolution of protoplanetary disks.

In order to obtain a complete census of the stellar and sub-stellar population, down to a few M$_{Jup}$ in the $\sim1$ Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the $1.4 \mu$m H$_2$O absorption feature and an adjacent line-free continuum region. Out of $4,504$ detected sources, $3,352$ (about $75\%$) appear fainter than m$_{130}=14$ (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass (M$\simeq 0.072 M_\odot$) at $\sim 1$ Myr. Of these, however, only $742$ sources have a negative F130M-139N color index, indicative of the presence of H$_2$O vapor in absorption, and can therefore be classified as bona-fide M and L dwarfs, with effective temperatures T$\lesssim 2850$ K at an assumed $1$ Myr cluster age. On our color-magnitude diagram, this population of sources with H$_2$O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M-F139N color index, and can be traced down to the sensitivity limit of our survey, m$_{130}\simeq 21.5$, corresponding to a $1$ Myr old $\simeq 3 $M$_{Jup}$, planetary mass object under about 2 magnitudes of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of $1, 2$ and $3$ Myr (down to $\sim 1 $M$_{Jup}$) fail to reproduce the observed H$_2$O color index at M$\lesssim 20 $M$_{Jup}$. We perform a Bayesian analysis to determine extinction, mass and effective temperature of each sub-stellar member of our sample, together with its membership probability.

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.

Recent ALMA surveys of protoplanetary disks have shown that for most disks the extent of the gas emission is greater than the extent of the thermal emission of the millimeter-sized dust. Both line optical depth and the combined effect of radially dependent grain growth and radial drift may contribute to this observed effect. For a sample of 10 disks from the Lupus survey we investigate how well dust-based models without radial dust evolution reproduce the observed 12CO outer radius, and determine whether radial dust evolution is required to match the observed gas-dust size difference. We used the thermochemical code DALI to obtain 12CO synthetic emission maps and measure gas and dust outer radii (Rco, Rmm) using the same methods as applied to the observations, which were compared to observations on a source-by-source basis. For 5 disks we find that the observed gas-dust size difference is larger than the gas-dust size difference due to optical depth, indicating that we need both dust evolution and optical depth effects to explain the observed gas-dust size difference. For the other 5 disks the observed gas-dust size difference can be explained using only line optical depth effects. We also identify 6 disks not included in our initial sample but part of a survey of the same star-forming region that show significant 12CO emission beyond 4 x Rmm. These disks, for which no Rco is available, likely have gas-dust size differences greater than 4 and are difficult to explain without substantial dust evolution. Our results suggest that radial drift and grain growth are common features among both bright and fain disks. The effects of radial drift and grain growth can be observed in disks where the dust and gas radii are significantly different, while more detailed models and deeper observations are needed to see this effect in disks with smaller differences.

Very recently, a second companion on wider orbit has been discovered around GQ Lup. This is a low-mass accreting star partially obscured by a disk seen at high inclination. If detected, this disk may be compared to the known disk around the primary. We detected this disk on archive HST and WISE data. The extended spectral energy distribution provided by these data confirms the presence of accretion from Halpha emission and UV excess, and shows an IR excess attributable to a warm disk. In addition, we resolved the disk on the HST images. This is found to be roughly aligned with the disk of the primary. Both of them are roughly aligned with the Lupus I dust filament containing GQ Lup.

The UV radiation field is a critical regulator of gas-phase chemistry in surface layers of disks around young stars. In an effort to understand the relationship between photocatalyzing UV radiation fields and gas emission observed at infrared and sub-mm wavelengths, we present an analysis of new and archival HST, Spitzer, ALMA, IRAM, and SMA data for five targets in the Lupus cloud complex and 14 systems in Taurus-Auriga. The HST spectra were used to measure LyA and FUV continuum fluxes reaching the disk surface, which are responsible for dissociating relevant molecular species (e.g. HCN, N2). Semi-forbidden C II] 2325 and UV-fluorescent H2 emission were also measured to constrain inner disk populations of C+ and vibrationally excited H2. We find a significant positive correlation between 14 micron HCN emission and fluxes from the FUV continuum and C II] 2325, consistent with model predictions requiring N2 photodissociation and carbon ionization to trigger the main CN/HCN formation pathways. We also report significant negative correlations between sub-mm CN emission and both C II] and FUV continuum fluxes, implying that CN is also more readily dissociated in disks with stronger FUV irradiation. No clear relationships are detected between either CN or HCN and LyA or UV-H2 emission. This is attributed to the spatial stratification of the various molecular species, which span several vertical layers and radii across the inner and outer disk. We expect that future observations with JWST will build on this work by enabling more sensitive IR surveys than were possible with Spitzer.

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.

Substellar companions at wide separation around stars hosting planets or brown dwarfs (BDs) yet close enough for their formation in the circumstellar disc are of special interest. In this letter we report the discovery of a wide (projected separation $\sim$16.0arcsec, or 2400 AU, and position angle 114.61$^\circ$) companion of the GQ Lup A-B system, most likely gravitationally bound to it. A VLT/X-Shooter spectrum shows that this star, 2MASS J15491331-3539118, is a bonafide low-mass ($\sim$0.15 M$_\odot$) young stellar object (YSO) with stellar and accretion/ejection properties typical of Lupus YSOs of similar mass, and with kinematics consistent with that of the GQ Lup A-B system. A possible scenario for the formation of the triple system is that GQ Lup A and 2MASS J15491331-3539118 formed by fragmentation of a turbulent core in the Lup I filament, while GQ Lup B, the BD companion of GQ Lup A at 0.7arcsec, formed in situ by the fragmentation of the circumprimary disc. The recent discoveries that stars form along cloud filaments would favour the scenario of turbulent fragmentation for the formation of GQ Lup A and 2MASS J15491331-3539118.

Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on timescales ~Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts also specific disc isochrones in the accretion rate - disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) timescale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data in the accretion rate - disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate.

The Young Suns Exoplanet Survey (YSES) consists of a homogeneous sample of 70 young, solar-mass stars located in the Lower Centaurus-Crux subgroup of the Scorpius-Centaurus association with an average age of $15\pm3\,$Myr. We report the detection of a co-moving companion around the K3IV star TYC 8998-760-1 (2MASSJ13251211-6456207) that is located at a distance of $94.6\pm0.3\,$pc using SPHERE/IRDIS on the VLT. Spectroscopic observations with VLT/X-SHOOTER constrain the mass of the star to $1.00\pm0.02\,M_{\odot}$ and an age of $16.7\pm1.4\,$Myr. The companion TYC 8998-760-1 b is detected at a projected separation of 1.71'', which implies a projected physical separation of $162\,$au. Photometric measurements ranging from $Y$ to $M$ band provide a mass estimate of $14\pm3\,M_\mathrm{jup}$ by comparison to BT-Settl and AMES-dusty isochrones, corresponding to a mass ratio of $q=0.013\pm0.003$ with respect to the primary. We rule out additional companions to TYC 8998-760-1 that are more massive than $12\,M_\mathrm{jup}$ and farther than $12\,$au away from the host. Future polarimetric and spectroscopic observations of this system with ground and space based observatories will facilitate testing of formation and evolution scenarios shaping the architecture of the circumstellar environment around this 'young Sun'.

RX~J1604.3-2130A is a young, dipper-type, variable star in the Upper Scorpius association, suspected to have an inclined inner disk with respect to its face-on outer disk. We study the eclipses to constrain the inner disk properties.We use time-resolved photometry from the Rapid Eye Mount telescope and Kepler2 data to study the multi-wavelength variability, and archival optical and IR data to track accretion, rotation, and changes in disk structure. The observations reveal details of the structure and matter transport through the inner disk. The eclipses show 5d quasi-periodicity, with the phase drifting in time and some periods showing increased/decreased eclipse depth and frequency. Dips are consistent with extinction by slightly processed dust grains in an inclined, irregularly-shaped inner disk locked to the star through two relatively stable accretion structures. The grains are located near the dust sublimation radius ($\sim$0.06 au) at the corotation radius, and can explain the shadows observed in the outer disk. The total mass (gas and dust) required to produce the eclipses and shadows is a few % of a Ceres mass. Such amount of mass is accreted/replenished by accretion in days to weeks, which explains the variability from period to period. Spitzer and WISE variability reveal variations in the dust content in the innermost disk on a few years timescale, which is consistent with small imbalances (compared to the stellar accretion rate) in the matter transport from the outer to the inner disk. A decrease in the accretion rate is observed at the times of less eclipsing variability and low mid-IR fluxes, confirming this picture. The v$sini$=16km/s confirms that the star cannot be aligned with the outer disk, but is likely close to equator-on and to be aligned with the inner disk. This anomalous orientation is a challenge for standard theories of protoplanetary disk formation.

We present new 890 $\mu m$ continuum ALMA observations of 5 brown dwarfs (BDs) with infrared excess in Lupus I and III -- which, in combination with 4 BDs previously observed, allowed us to study the mm properties of the full known BD disk population of one star-forming region. Emission is detected in 5 out of the 9 BD disks. Dust disk mass, brightness profiles and characteristic sizes of the BD population are inferred from continuum flux and modeling of the observations. Only one source is marginally resolved, allowing for the determination of its disk characteristic size. We conduct a demographic comparison between the properties of disks around BDs and stars in Lupus. Due to the small sample size, we cannot confirm or disprove if the disk mass over stellar mass ratio drops for BDs, as suggested for Ophiuchus. Nevertheless, we find that all detected BD disks have an estimated dust mass between 0.2 and 3.2 $M_{\bigoplus}$; these results suggest that the measured solid masses in BD disks can not explain the observed exoplanet population, analogous to earlier findings on disks around more massive stars. Combined with the low estimated accretion rates, and assuming that the mm-continuum emission is a reliable proxy for the total disk mass, we derive ratios of $\dot{M}_{\mathrm{acc}} / M_{\mathrm{disk}}$ significantly lower than in disks around more massive stars. If confirmed with more accurate measurements of disk gas masses, this result could imply a qualitatively different relationship between disk masses and inward gas transport in BD disks.

While planets are commonly discovered around main-sequence stars, the processes leading to their formation are still far from being understood. Current planet population synthesis models, which aim to describe the planet formation process from the protoplanetary disk phase to the time exoplanets are observed, rely on prescriptions for the underlying properties of protoplanetary disks where planets form and evolve. The recent development in measuring disk masses and disk-star interaction properties, i.e., mass accretion rates, in large samples of young stellar objects demand a more careful comparison between the models and the data. We performed an initial critical assessment of the assumptions made by planet synthesis population models by looking at the relation between mass accretion rates and disk masses in the models and in the currently available data. We find that the currently used disk models predict mass accretion rate in line with what is measured, but with a much lower spread of values than observed. This difference is mainly because the models have a smaller spread of viscous timescales than what is needed to reproduce the observations. We also find an overabundance of weakly accreting disks in the models where giant planets have formed with respect to observations of typical disks. We suggest that either fewer giant planets have formed in reality or that the prescription for planet accretion predicts accretion on the planets that is too high. Finally, the comparison of the properties of transition disks with large cavities confirms that in many of these objects the observed accretion rates are higher than those predicted by the models. On the other hand, PDS70, a transition disk with two detected giant planets in the cavity, shows mass accretion rates well in line with model predictions.

We investigate ongoing accretion activity in young stars in the TW Hydrae association (TWA, ~8-10 Myr), an ideal target to probe the final stages of disk accretion down to brown dwarf masses. Our sample comprises eleven TWA members with infrared excess, amounting to 85% of the total TWA population with disks, with spectral types between M0 and M9, and masses between 0.58 and 0.02 Msol. We employed homogeneous spectroscopic data from 300 to 2500 nm, obtained with X-shooter, to derive individual extinction, stellar parameters, and accretion parameters simultaneously. We then examined Balmer lines and forbidden emission lines to probe the physics of the star-disk interaction environment. We detected signatures of ongoing accretion for 70% of our TWA targets. This implies a fraction of accretors of 13-17% across the entire TWA (accounting for potentially accreting members not included in our survey). The spectral emission associated with these stars reveals a more evolved stage of these accretors compared to younger PMS populations: (i) a large fraction (~50%) exhibit nearly symmetric, narrow H_alpha line profiles; (ii) over 80% exhibit Balmer decrements consistent with moderate accretion and optically thin emission; (iii) less than a third exhibit forbidden line emission in [O I] 6300A, indicative of winds and outflows activity. However, the distribution in accretion rates (Mdot) derived for the TWA sample follows closely that of younger regions for Mstar~0.1-0.3 Msun. An overall correlation between Mdot and Mstar is detected, best reproduced by Mdot~Mstar^(2.1+/-0.5). At least in the lowest Mstar regimes, stars that still retain a disk at ages ~8-10 Myr are found to exhibit statistically similar, albeit moderate, accretion levels as those measured around younger objects. This slow Mdot evolution may be associated with longer evolutionary timescales of disks around low-mass stars.

The mass accretion rate is a crucial parameter for the study of the evolution of accretion discs around young low-mass stellar and substellar objects (YSOs). We revisit the case of MYLup, an object where VLT/X-Shooter data suggested a negligible mass accretion rate, and show it to be accreting on a level similar to other ClassII YSOs in Lupus based on HST observations. In our HST-COS and -STIS spectra, we find many emission lines, as well as substantial FUV continuum excess emission, which can be ascribed to active accretion. The total luminosity of the CIV 1549A doublet is 4.1$\times10^{-4}$ Lsun. Using scalings between accretion luminosity, Lacc, and CIV luminosity from the literature, we derive Lacc $\sim$2$\times10^{-1}$ Lsun, which is more than an order of magnitude higher than the upper limit estimated from the X-Shooter observations. We discuss possible reasons for the X-Shooter-HST discrepancy, the most plausible being that the low contrast between the continuum excess emission and the photospheric+chromospheric emission at optical wavelengths in MYLup hampered detection of excess emission. The luminosity of the FUV continuum and CIV lines, strong H$_{2}$ fluorescence, and a "1600 A Bump" place MYLup in the class of accreting objects with gas-rich discs. So far, MYLup is the only peculiar case in which a significant difference between the HST and X-Shooter Macc estimates exists that is not ascribable to variability. The mass accretion rate inferred from the revisited Lacc estimate is Macc $\sim$ 1( $^{+1.5}_{-0.5}$ )$\times10^{-8}$ Msun/yr. This value is consistent with the typical value derived for accreting YSOs of similar mass in Lupus and points to less clearing of the inner disc than indicated by near- and mid-infrared observations. This is confirmed by ALMA data, which show that the gaps and rings seen in the sub-millimetre are relatively shallow.

The impact of stellar multiplicity on the evolution of planet-forming disks is still the subject of debate. Here we present and analyze disk structures around ten multiple stellar systems that were included in an unbiased, high spatial resolution survey performed with ALMA of 32 protoplanetary disks in the Taurus star-forming region. At the unprecedented spatial resolution of ~0.12" we detect and spatially resolve the disks around all primary stars, and those around eight secondary and one tertiary star. The dust radii of disks around multiple stellar systems are smaller than those around single stars in the same stellar mass range and in the same region. The disks in multiple stellar systems also show a steeper decay of the millimeter continuum emission at the outer radius than disks around single stars, suggestive of the impact of tidal truncation on the shape of the disks in multiple systems. However, the observed ratio between the dust disk radii and the observed separation of the stars in the multiple systems is consistent with analytic predictions of the effect of tidal truncation only if the eccentricities of the binaries are rather high (typically >0.5), or if the observed dust radii are a factor of two smaller than the gas radii, as is typical for isolated systems. Similar high-resolution studies targeting the gaseous emission from disks in multiple stellar systems are required to resolve this question.

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.

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)

The large majority of protoplanetary disks have very compact ($\lesssim15\,$AU) continuum emission at mm wavelengths. However, high angular resolution observations that resolve these small disks are still lacking, due to their intrinsically fainter emission compared with large bright disks. In this letter, we present $1.3\,$mm ALMA data of the faint ($\sim10\,$mJy) disk orbiting the TTauri star CX Tau at a resolution of $\sim40\,$mas, $\sim5\,$AU in diameter. The mm-dust disk is compact, with a 68$\%$ enclosing flux radius of 14$\,$AU, and the intensity profile exhibits a sharp drop between 10-20$\,$AU, and a shallow tail between 20-40$\,$AU. No clear signatures of substructure in the dust continuum are observed, down to the same sensitivity level of the DSHARP large program. However, the angular resolution does not allow to detect substructures at a scale of the disk aspect ratio in the inner regions. The radial intensity profile resembles well the inner regions of more extended disks imaged at the same resolution in DSHARP, but with no rings present in the outer disk. No inner cavity is detected, even though the disk has been classified as a transition disk from the spectral energy distribution in the near infrared. The emission of $^{12}$CO is much more extended, with a 68$\%$ enclosing flux radius of 75$\,$AU. The large difference of the mm dust and gas extents ($>5$) strongly points to the occurrence of radial drift, and matches well the predictions of theoretical models.

The dense clusters within the Serpens Molecular Cloud are among the most active regions of nearby star formation. In this paper, we use Gaia DR2 parallaxes and proper motions to statistically measure 1167 kinematic members of Serpens, few of which were previously identified, to evaluate the star formation history of the complex. The optical members of Serpens are concentrated in three distinct groups located at 380 to 480 pc; the densest clusters are still highly obscured by optically-thick dust and have few optical members. The total population of young stars and protostars in Serpens is at least 2000 stars, including past surveys that were most sensitive to protostars and disks, and may be far higher. Distances to dark clouds measured from deficits in star counts are consistent with the distances to the optical star clusters. The Serpens Molecular Cloud is seen in the foreground of the Aquila Rift, dark clouds located at 600 to 700 pc, and behind patchy extinction, here called the Serpens Cirrus, located at ~250 pc. Based on the lack of a distributed population of older stars, the star formation rate throughout the Serpens Molecular Cloud increased by at least a factor of 20 within the past ~5 Myr. The optically bright stars in Serpens Northeast are visible because their natal molecular cloud has been eroded and not because they were flung outwards from a central factory of star formation. The separation between subclusters of 20 to 100 pc and the absence of an older population leads to speculation that an external forcing was needed to trigger the active star formation.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

RW Aur A is a CTTS that has suddenly undergone three major dimming events since 2010. We aim to understand the dimming properties, examine accretion variability, and derive the physical properties of the inner disc traced by the CO ro-vibrational emission at NIR wavelengths (2.3 mic). We compared two epochs of X-Shooter observations, during and after the dimming. We modelled the rarely detected CO bandhead emission in both epochs to examine whether the inner disc properties had changed. The SED was used to derive the extinction properties of the dimmed spectrum and compare the infrared excess between the two epochs. Lines tracing accretion were used to derive the mass accretion rate in both states. The CO originates from a region with physical properties of T=3000 K, N$_{CO}$=1x10$^{21}$ cm$^{-2}$ and vsini=113 km/s. The extinction properties of the dimming layer were derived with the effective optical depth ranging from teff 2.5-1.5 from the UV to the NIR. The inferred mass accretion rate Macc is $1.5x 10^{-8}$ Msun/yr and $\sim 2x 10^{-8}$ Msun/yr after and during the dimming respectively. By fitting the SED, additional emission is observed in the IR during the dimming event from dust grains with temperatures of 500-700K. The physical conditions traced by the CO are similar for both epochs, indicating that the inner gaseous disc properties do not change during the dimming events. The extinction curve is flatter than that of the ISM, and large grains of a few hundred microns are thus required. When we correct for the observed extinction, Macc is constant in the two epochs, suggesting that the accretion is stable and therefore does not cause the dimming. The additional hot emission in the NIR is located at about 0.5 au from the star. The dimming events could be due to a dust-laden wind, a severe puffing-up of the inner rim, or a perturbation caused by the recent star-disc encounter.

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters.

Aims. R Coronae Australis (R CrA) is the brightest star of the Coronet nebula of the Corona Australis (CrA) star forming region. It has very red colors, probably due to dust absorption and it is strongly variable. High contrast instruments allow for an unprecedented direct exploration of the immediate circumstellar environment of this star. Methods. We observed R CrA with the near-IR channels (IFS and IRDIS) of SPHERE at VLT. In this paper, we used four different epochs, three of them from open time observations while one is from the SPHERE guaranteed time. The data were reduced using the DRH pipeline and the SPHERE Data Center. On the reduced data we implemented custom IDL routines with the aim to subtract the speckle halo.We have also obtained pupil-tracking H-band (1.45-1.85 micron) observations with the VLT/SINFONI near-infrared medium-resolution (R~3000) spectrograph. Results. A companion was found at a separation of 0.156" from the star in the first epoch and increasing to 0.18400 in the final one. Furthermore, several extended structures were found around the star, the most noteworthy of which is a very bright jet-like structure North-East from the star. The astrometric measurements of the companion in the four epochs confirm that it is gravitationally bound to the star. The SPHERE photometry and the SINFONI spectrum, once corrected for extinction, point toward an early M spectral type object with a mass between 0.3 and 0.55 M?. The astrometric analyis provides constraints on the orbit paramenters: e~0.4, semi-major axis at 27-28 au, inclination of ~ 70{\deg} and a period larger than 30 years. We were also able to put constraints of few MJup on the mass of possible other companions down to separations of few tens of au.

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.

Low brightness dips have been recently observed in images of protoplanetary disks, and they are believed to be shadows by the inner disk. We present VLT/SPHERE polarimetric differential imaging of the transition disk around the dipper star RX J1604.3-2130. We gathered 11 epochs that cover a large temporal baseline, to search for variability over timescales of years, months, weeks, and days. Our observations unambiguously reveal two dips along an almost face-on narrow ring (with a width of $\sim$20 au), and the location of the peak of this ring is at $\sim$65 au. The ring lies inside the ring-like structure observed with ALMA, which peaks at $\sim$83 au. This segregation can result from particle trapping in pressure bumps, potentially due to planet(s). We find that the dips are variable, both in morphology and in position. The eastern dip, at a position angle (PA) of $\sim$83.7$\pm$13.7$^\circ$, has an amplitude that varies between 40% to 90%, and its angular width varies from 10$^{\circ}$ to 34$^{\circ}$. The western dip, at a PA of $\sim$265.90$\pm$13.0$^\circ$, is more variable, with amplitude and width variations of 31% to 95% and 12$^{\circ}$ to 53$^{\circ}$, respectively. The separation between the dips is 178.3$^{\circ}\pm$14.5$^{\circ}$, corresponding to a large misalignment between the inner and outer disk, supporting the classification of J1604 as an aperiodic dipper. The variability indicates that the innermost regions are highly dynamic, possibly due to a massive companion or to a complex magnetic field topology.

The cyanide radical CN is abundant in protoplanetary disks, with line fluxes often comparable to those of $^{13}$CO. It is known to be sensitive to UV irradiation of the upper disk atmosphere, with models predicting ring-shaped emission. We seek to characterize the CN emission from 94 Class II disks in the Lupus star forming region, compare it to observations in other regions, and interpret our observations with a grid of models. The CN emission morphology is discussed for two primordial disks, Sz 71 and Sz 98, and modeled in more detail. ALMA observed CN N=3-2 in Lupus disks down to sensitivities better than previous surveys. Models constructed with the physico-chemical code DALI are used to study the integrated fluxes of the disks and resolved emission of CN in disks without (dust) substructures. We find that CN N=3-2 is bright, and detected in 38% of sources, but its disk-integrated flux is not strongly correlated to either $^{13}$CO or continuum flux. Compared to pre-ALMA single-dish surveys, no significant difference in the CN flux distributions in Lupus and Taurus-Auriga is found, although $\rho$ Ophiuchus disks may be fainter on average. We find ring-shaped CN emission with peak radii of ~50 AU in two resolved disks. We conclude that a large fraction of sources is faint in CN; only exponential gas surface density cutoffs at $R_{c} \leq 15$ AU can reconcile observations with models. This is the first observational evidence of such a compact gas disk population in Lupus. Absolute intensities and the emission morphology of CN are reproduced by DALI models without the need for any continuum substructure; they are unrelated to the CO snowline location. These observations and the successful modeling of these rings provide a new probe of the structure and conditions in disks, and particularly their incident UV radiation field, if disk size is determined by the data.

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.

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.

When and how planets form in protoplanetary disks is still a topic of discussion. Exoplanet detection surveys and protoplanetary disk surveys are now providing results that allow us to have new insights. We collect the masses of confirmed exoplanets and compare their dependence with stellar mass with the same dependence for protoplanetary disk masses measured in ~1-3 Myr old star-forming regions. The latter are recalculated by us using the new estimates of their distances derived from Gaia DR2 parallaxes. We note that single and multiple exoplanetary systems form two different populations, probably pointing to a different formation mechanism for massive giant planets around very low mass stars. While expecting that the mass in exoplanetary systems is much lower than the measured disk masses, we instead find that exoplanetary systems masses are comparable or higher than the most massive disks. This same result is found also by converting the measured planet masses into heavy-element content (core masses for the giant planets and full masses for the super-Earth systems) and by comparing this value with the disk dust masses. Unless disk dust masses are heavily underestimated, this is a big conundrum. An extremely efficient recycling of dust particles in the disk cannot solve this conundrum. This implies that either the cores of planets have formed very rapidly (<0.1-1 Myr) and large amount of gas is expelled on the same timescales from the disk, or that disks are continuously replenished of fresh planet-forming material from the environment. These hypotheses can be tested by measuring disk masses in even younger targets and by better understanding if and how the disks are replenished by their surroundings.

Extensive surveys of star-forming regions with Spitzer have revealed populations of disk-bearing young stellar objects. These have provided crucial constraints, such as the timescale of dispersal of protoplanetary disks, obtained by carefully combining infrared data with spectroscopic or X-ray data. While observations in various regions agree with the general trend of decreasing disk fraction with age, the Lupus V and VI regions appeared to have been at odds, having an extremely low disk fraction. Here we show, using the recent Gaia data release 2 (DR2), that these extremely low disk fractions are actually due to a very high contamination by background giants. Out of the 83 candidate young stellar objects (YSOs) in these clouds observed by Gaia, only five have distances of 150 pc, similar to YSOs in the other Lupus clouds, and have similar proper motions to other members in this star-forming complex. Of these five targets, four have optically thick (Class II) disks. On the one hand, this result resolves the conundrum of the puzzling low disk fraction in these clouds, while, on the other hand, it further clarifies the need to confirm the Spitzer selected diskless population with other tracers, especially in regions at low galactic latitude like Lupus V and VI. The use of Gaia astrometry is now an independent and reliable way to further assess the membership of candidate YSOs in these, and potentially other, star-forming regions.

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.

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.

We want to study whether the astrometric and photometric accuracies obtained for the Carte du Ciel plates digitized with a commercial digital camera are high enough for scientific exploitation of the plates. We use a digital camera Canon EOS~5Ds, with a 100mm macrolens for digitizing. We analyze six single-exposure plates and four triple-exposure plates from the Helsinki zone of Carte du Ciel (+39 degr < delta < 47 degr). Each plate is digitized using four images, with a significant central area being covered twice for quality control purposes. The astrometric calibration of the digitized images is done with the data from the Gaia TGAS (Tycho-Gaia Astrometric Solution) of the first Gaia data release (Gaia DR1), Tycho-2, HSOY (Hot Stuff for One Year), UCAC5 (USNO CCD Astrograph Catalog), and PMA catalogs. The best astrometric accuracy is obtained with the UCAC5 reference stars. The astrometric accuracy for single-exposure plates is sigma(R.A.)=0.16" and sigma(Dec.)=0.15" expressed as a Gaussian deviation of the astrometric residuals. For triple-exposure plates the astrometric accuracy is sigma(R.A.)=0.12" and sigma(Dec.)=0.13". The 1-sigma uncertainty of photometric calibration is about 0.28 mag and 0.24 mag for single- and triple-exposure plates, respectively. We detect the photographic adjacency (Kostinsky) effect in the triple-exposure plates. We show that accuracies at least of the level of scanning machines can be achieved with a digital camera, without any corrections for possible distortions caused by our instrumental setup. This method can be used to rapidly and inexpensively digitize and calibrate old photographic plates enabling their scientific exploitation.

A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $\alpha_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

In the present study we aim to investigate the circumstellar environment of the spectroscopic binary T Tauri star CS Cha. From unresolved mid- to far-infrared photometry it is predicted that CS Cha hosts a disk with a large cavity. In addition, SED modeling suggests significant dust settling, pointing towards an evolved disk that may show signs of ongoing or completed planet formation. We observed CS Cha with the high contrast imager VLT/SPHERE in polarimetric differential imaging mode to resolve the circumbinary disk in near infrared scattered light. These observations were followed-up by VLT/NACO L-band observations and complemented by archival VLT/NACO K-band and HST/WFPC2 I-band data. We resolve the compact circumbinary disk around CS Cha for the first time in scattered light. We find a smooth, low inclination disk with an outer radius of $\sim$55 au (at 165 pc). We do not detect the inner cavity but find an upper limit for the cavity size of $\sim$15 au. Furthermore, we find a faint co-moving companion with a projected separation of 210 au from the central binary outside of the circumbinary disk. The companion is detected in polarized light and shows an extreme degree of polarization (13.7$\pm$0.4 \% in J-band). The companion's J- and H-band magnitudes are compatible with masses of a few M$_\mathrm{Jup}$. However, K-, L- and I-band data draw this conclusion into question. We explore with radiative transfer modeling whether an unresolved circum-companion disk can be responsible for the high polarization and complex photometry. We find that the set of observations is best explained by a heavily extincted low mass ($\sim 20 \mathrm{M}_\mathrm{Jup}$) brown dwarf or high mass planet with an unresolved disk and dust envelope.

We analyze the dust morphology of 29 transition disks (TDs) observed with ALMA at (sub-) millimeter-emission. We perform the analysis in the visibility plane to characterize the total flux, cavity size, and shape of the ring-like structure. First, we found that the $M_{\rm{dust}}-M_\star$ relation is much flatter for TDs than the observed trends from samples of class II sources in different star forming regions. This relation demonstrates that cavities open in high (dust) mass disks, independent of the stellar mass. The flatness of this relation contradicts the idea that TDs are a more evolved set of disks. Two potential reasons (not mutually exclusive) may explain this flat relation: the emission is optically thick or/and millimeter-sized particles are trapped in a pressure bump. Second, we discuss our results of the cavity size and ring width in the context of different physical processes for cavity formation. Photoevaporation is an unlikely leading mechanism for the origin of the cavity of any of the targets in the sample. Embedded giant planets or dead zones remain as potential explanations. Although both models predict correlations between the cavity size and the ring shape for different stellar and disk properties, we demonstrate that with the current resolution of the observations, it is difficult to obtain these correlations. Future observations with higher angular resolution observations of TDs with ALMA will help to discern between different potential origins of cavities in TDs.

We present ALMA Band 6 observations of a complete sample of protoplanetary disks in the young (1-3 Myr) Lupus star-forming region, covering the 1.33 mm continuum and the 12CO, 13CO, and C18O J=2-1 lines. The spatial resolution is 0.25 arcsec with a medium 3-sigma continuum sensitivity of 0.30 mJy, corresponding to M_dust ~ 0.2 M_earth. We apply "Keplerian masking" to enhance the signal-to-noise ratios of our 12CO zero-moment maps, enabling measurements of gas disk radii for 22 Lupus disks; we find that gas disks are universally larger than mm dust disks by a factor of two on average, likely due to a combination of the optically thick gas emission as well as the growth and inward drift of the dust. Using the gas disk radii, we calculate the dimensionless viscosity parameter, alpha_visc, finding a broad distribution and no correlations with other disk or stellar parameters, suggesting that viscous processes have not yet established quasi-steady states in Lupus disks. By combining our 1.33 mm continuum fluxes with our previous 890 micron continuum observations, we also calculate the mm spectral index, alpha_mm, for 70 Lupus disks; we find an anti-correlation between alpha_mm and mm flux for low-mass disks (M_dust < 5), followed by a flattening as disks approach alpha_mm = 2, which could indicate faster grain growth in higher-mass disks, but may also reflect their larger optically thick components. In sum, this work demonstrates the continuous stream of new insights into disk evolution and planet formation that can be gleaned from unbiased ALMA disk surveys.

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.

Disk warps around classical T Tauri stars (CTTS) can periodically obscure the central star for some viewing geometries. For these so-called AA Tau-like variables, the obscuring material is located in the inner disk and absorption spectroscopy allows one to characterize its dust and gas content. Since the observed emission from CTTS consists of several components (photospheric, accretion, jet, and disk emission), which can all vary with time, it is generally challenging to disentangling disk features from emission variability. Multi-epoch, flux-calibrated, broadband spectra provide us with the necessary information to cleanly separate absorption from emission variability. We applied this method to three epochs of VLT/X-Shooter spectra of the CTTS V354 Mon (CSI Mon-660) located in NGC 2264 and find that: (a) the accretion emission remains virtually unchanged between the three epochs; (b) the broadband flux evolution is best described by disk material obscuring part of the star, and (c) the Na and K gas absorption lines show only a minor increase in equivalent width during phases of high dust extinction. The limits on the absorbing gas column densities indicate a low gas-to-dust ratio in the inner disk, less than a tenth of the ISM value. We speculate that the evolutionary state of V354 Mon, rather old with a low accretion rate, is responsible for the dust excess through an evolution toward a dust dominated disk or through the fragmentation of larger bodies that drifted inward from larger radii in a still gas dominated disk.

We present a study of molecular gas in the inner disk $\left(r < 20 \, \text{AU} \right)$ around RY Lupi, with spectra from HST-COS, HST-STIS, and VLT-CRIRES. We model the radial distribution of flux from hot gas in a surface layer between $r = 0.1-10$ AU, as traced by Ly$\alpha$-pumped H$_2$. The result shows H$_2$ emission originating in a ring centered at $\sim$3 AU that declines within $r < 0.1$ AU, which is consistent with the behavior of disks with dust cavities. An analysis of the H$_2$ line shapes shows that a two-component Gaussian profile $\left(\text{FWHM}_{broad, H_2} = 105 \pm 15 \, \text{km s}^{-1}; \, \text{FWHM}_{narrow, H_2} = 43 \pm 13 \, \text{km s}^{-1} \right)$ is statistically preferred to a single-component Gaussian. We interpret this as tentative evidence for gas emitting from radially separated disk regions $\left( \left \langle r_{broad, H_2} \right \rangle \sim 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, H_2} \right \rangle \sim 3 \pm 2 \, \text{AU} \right)$. The 4.7 $\mu$m $^{12}$CO emission lines are also well fit by two-component profiles $\left( \left \langle r_{broad, CO} \right \rangle = 0.4 \pm 0.1 \, \text{AU}; \, \left \langle r_{narrow, CO} \right \rangle = 15 \pm 2 \, \text{AU} \right)$. We combine these results with 10 $\mu$m observations to form a picture of gapped structure within the mm-imaged dust cavity, providing the first such overview of the inner regions of a young disk. The HST SED of RY Lupi is available online for use in modeling efforts.

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.

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.

We use the first data release of the Gaia mission to explore the three dimensional arrangement and the age ordering of the many stellar groups towards the Orion OB association, aiming at a new classification and characterization of the stellar population. We make use of the parallaxes and proper motions provided in the Tycho Gaia Astrometric Solution (TGAS) sub-set of the Gaia catalogue, and of the combination of Gaia and 2MASS photometry. In TGAS we find evidence for the presence of a young population, at a parallax $\varpi \sim 2.65 \, \mathrm{mas}$, loosely distributed around some known clusters: 25 Ori, $\epsilon$ Ori and $\sigma$ Ori, and NGC 1980 ($\iota$ Ori). The low mass counterpart of this population is visible in the color-magnitude diagrams constructed by combining Gaia and 2MASS photometry. We study the density distribution of the young sources in the sky. We find the same groups as in TGAS, and also some other density enhancements that might be related to the recently discovered Orion X group, the Orion dust ring, and to the $\lambda$ Ori complex. We estimate the ages of this population and we infer the presence of an age gradient going from 25 Ori (13-15 Myr) to the ONC (1-2 Myr). We confirm this age ordering by repeating the Bayesian fit using the Pan-STARRS1 data. The estimated ages towards the NGC 1980 cluster span a broad range of values. This can either be due to the presence of two populations coming from two different episodes of star formation or to a large spread along the line of sight of the same population. Our results form the first step towards using the Gaia data to unravel the complex star formation history of the Orion region in terms of the different star formation episodes, their duration, and their effects on the surrounding interstellar medium.

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.

We have analysed the [OI]6300 A line in a sample of 131 young stars with discs in the Lupus, Chamaeleon and signa Orionis star forming regions, observed with the X-shooter spectrograph at VLT. The stars have mass accretion rates spanning from 10^{-12} to 10^{-7} Mo/yr. The line profile was deconvolved into a low velocity component (LVC, < 40 km/s) and a high velocity component (HVC, > 40 km/s ), originating from slow winds and high velocity jets, respectively. The LVC is by far the most frequent component, with a detection rate of 77%, while only 30% of sources have a HVC. The [OI]6300 luminosity of both the LVC and HVC, when detected, correlates with stellar and accretion parameters of the central sources (i.e. Lstar , Mstar , Lacc , Macc), with similar slopes for the two components. The line luminosity correlates better with the accretion luminosity than with the stellar luminosity or stellar mass. We suggest that accretion is the main drivers for the line excitation and that MHD disc-winds are at the origin of both components. In the sub-sample of Lupus sources observed with ALMA a relationship is found between the HVC peak velocity and the outer disc inclination angle, as expected if the HVC traces jets ejected perpendicularly to the disc plane. Mass loss rates measured from the HVC span from ~ 10^{-13} to ~10^{-7} Mo/yr. The corresponding Mloss/Macc ratio ranges from ~0.01 to ~0.5, with an average value of 0.07. However, considering the upper limits on the HVC, we infer a ratio < 0.03 in more than 40% of sources. We argue that most of these sources might lack the physical conditions needed for an efficient magneto-centrifugal acceleration in the star-disc interaction region. Systematic observations of populations of younger stars, that is, class 0/I, are needed to explore how the frequency and role of jets evolve during the pre-main sequence phase.

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.

We investigate the relation between disk mass and mass accretion rate to constrain the mechanism of angular momentum transport in protoplanetary disks. Dust mass and mass accretion rate in Chamaeleon I are correlated with a slope close to linear, similar to the one recently identified in Lupus. We investigate the effect of stellar mass and find that the intrinsic scatter around the best-fit Mdust-Mstar and Macc-Mstar relations is uncorrelated. Disks with a constant alpha viscosity can fit the observed relations between dust mass, mass accretion rate, and stellar mass, but over-predict the strength of the correlation between disk mass and mass accretion rate when using standard initial conditions. We find two possible solutions. 1) The observed scatter in Mdust and Macc is not primoridal, but arises from additional physical processes or uncertainties in estimating the disk gas mass. Most likely grain growth and radial drift affect the observable dust mass, while variability on large time scales affects the mass accretion rates. 2) The observed scatter is primordial, but disks have not evolved substantially at the age of Lupus and Chamaeleon I due to a low viscosity or a large initial disk radius. More accurate estimates of the disk mass and gas disk sizes in a large sample of protoplanetary disks, either through direct observations of the gas or spatially resolved multi-wavelength observations of the dust with ALMA, are needed to discriminate between both scenarios or to constrain alternative angular momentum transport mechanisms such as MHD disk winds.

In this paper, we compare simple viscous diffusion models for the disc evolution with the results of recent surveys of the properties of young protoplanetary discs. We introduce the useful concept of `disc isochrones' in the accretion rate - disc mass plane and explore a set of Montecarlo realization of disc initial conditions. We find that such simple viscous models can provide a remarkable agreement with the available data in the Lupus star forming region, with the key requirement that the average viscous evolutionary timescale of the discs is comparable to the cluster age. Our models produce naturally a correlation between mass accretion rate and disc mass that is shallower than linear, contrary to previous results and in agreement with observations. We also predict that a linear correlation, with a tighter scatter, should be found for more evolved disc populations. Finally, we find that such viscous models can reproduce the observations in the Lupus region only in the assumption that the efficiency of angular momentum transport is a growing function of radius, thus putting interesting constraints on the nature of the microscopic processes that lead to disc accretion.

Understanding how young stars and their circumstellar disks form and evolve is key to explain how planets form. The evolution of the star and the disk is regulated by different processes, both internal to the system or related to their environment. The former include accretion of material onto the central star, wind emission, and photoevaporation of the disk due to high-energy radiation from the central star. These are best studied spectroscopically, and the distance to the star is a key parameter in all these studies. Here we present new estimates of the distance to a complex of nearby star-forming clouds obtained combining TGAS distances with measurement of extinction on the line of sight. Furthermore, we show how we plan to study the effects of the environment on the evolution of disks with Gaia, using a kinematic modelling code we have developed to model young star-forming regions.

The formation of planets strongly depends on the total amount as well as on the spatial distribution of solids in protoplanetary disks. Thanks to the improvements in resolution and sensitivity provided by ALMA, measurements of the surface density of mm-sized grains are now possible on large samples of disks. Such measurements provide statistical constraints that can be used to inform our understanding of the initial conditions of planet formation. We analyze spatially resolved observations of 36 protoplanetary disks in the Lupus star forming complex from our ALMA survey at 890 micron, aiming to determine physical properties such as the dust surface density, the disk mass and size and to provide a constraint on the temperature profile. We fit the observations directly in the uv-plane using a two-layer disk model that computes the 890 micron emission by solving the energy balance at each disk radius. For 22 out of 36 protoplanetary disks we derive robust estimates of their physical properties. The sample covers stellar masses between ~0.1 and ~2 Solar masses, and we find no trend between the average disk temperatures and the stellar parameters. We find, instead, a correlation between the integrated sub-mm flux (a proxy for the disk mass) and the exponential cut-off radii (a proxy of the disk size) of the Lupus disks. Comparing these results with observations at similar angular resolution of Taurus-Auriga/Ophiuchus disks found in literature and scaling them to the same distance, we observe that the Lupus disks are generally fainter and larger at a high level of statistical significance. Considering the 1-2 Myr age difference between these regions, it is possible to tentatively explain the offset in the disk mass/disk size relation with viscous spreading, however with the current measurements other mechanisms cannot be ruled out.

With the purpose of performing a homogeneous determination of elemental abundances for members of the Lupus T association, we analyzed three chemical elements: lithium, iron, and barium. The aims were: to derive the Li abundance for ~90% of known class II stars in the Lupus I, II, III, IV clouds; to perform chemical tagging of a region where few Fe abundance measurements have been obtained in the past, and no determination of the Ba content has been done up to now. We also investigated possible Ba enhancement, as this element has become increasingly interesting in the last years following the evidence of Ba over-abundance in young clusters, the origin of which is still unknown. Using X-shooter@VLT, we analyzed the spectra of 89 cluster members, both class II and III stars. We measured the strength of the Li line and derived the abundance of this element through equivalent width measurements and curves of growth. For six class II stars we also measured the Fe and Ba abundances using the spectral synthesis and the code MOOG. The veiling contribution was taken into account for all three elements. We find a dispersion in the strength of the Li line at low Teff and identify three targets with severe Li depletion. The nuclear age inferred for these highly Li-depleted stars is around 15 Myr, which exceeds the isochronal one. As in other star-forming regions, no metal-rich members are found in Lupus, giving support to a recent hypothesis that the Fe abundance distribution of most of the nearby young regions could be the result of a common and widespread star formation episode involving the Galactic thin disk. We find that Ba is over-abundant by ~0.7 dex with respect to the Sun. Since current theoretical models cannot reproduce this Ba abundance pattern, we investigated whether this unusually large Ba content might be related to effects due to stellar parameters, stellar activity, and accretion.

The dependence of the mass accretion rate on the stellar properties is a key constraint for star formation and disk evolution studies. Here we present a study of a sample of stars in the Chamaeleon I star forming region carried out using the VLT/X-Shooter spectrograph. The sample is nearly complete down to M~0.1Msun for the young stars still harboring a disk in this region. We derive the stellar and accretion parameters using a self-consistent method to fit the broad-band flux-calibrated medium resolution spectrum. The correlation between the accretion luminosity to the stellar luminosity, and of the mass accretion rate to the stellar mass in the logarithmic plane yields slopes of 1.9 and 2.3, respectively. These slopes and the accretion rates are consistent with previous results in various star forming regions and with different theoretical frameworks. However, we find that a broken power-law fit, with a steeper slope for stellar luminosity smaller than ~0.45 Lsun and for stellar masses smaller than ~ 0.3 Msun, is slightly preferred according to different statistical tests, but the single power-law model is not excluded. The steeper relation for lower mass stars can be interpreted as a faster evolution in the past for accretion in disks around these objects, or as different accretion regimes in different stellar mass ranges. Finally, we find two regions on the mass accretion versus stellar mass plane empty of objects. One at high mass accretion rates and low stellar masses, which is related to the steeper dependence of the two parameters we derived. The second one is just above the observational limits imposed by chromospheric emission. This empty region is located at M~0.3-0.4Msun, typical masses where photoevaporation is known to be effective, and at mass accretion rates ~10^-10 Msun/yr, a value compatible with the one expected for photoevaporation to rapidly dissipate the inner disk.

Studies of the formation and evolution of young stars and their disks rely on the knowledge of the stellar parameters of the young stars. The derivation of these parameters is commonly based on comparison with photospheric template spectra. Furthermore, chromospheric emission in young active stars impacts the measurement of mass accretion rates, a key quantity to study disk evolution. Here we derive stellar properties of low-mass pre-main sequence stars without disks, which represent ideal photospheric templates for studies of young stars. We also use these spectra to constrain the impact of chromospheric emission on the measurements of mass accretion rates. The spectra in reduced, flux-calibrated, and corrected for telluric absorption form are made available to the community. We derive the spectral type for our targets by analyzing the photospheric molecular features present in their VLT/X-Shooter spectra by means of spectral indices and comparison of the relative strength of photospheric absorption features. We also measure effective temperature, gravity, projected rotational velocity, and radial velocity from our spectra by fitting them with synthetic spectra with the ROTFIT tool. The targets have negligible extinction and spectral type from G5 to M8. We perform synthetic photometry on the spectra to derive the typical colors of young stars in different filters. We measure the luminosity of the emission lines present in the spectra and estimate the noise due to chromospheric emission in the measurements of accretion luminosity in accreting stars. We provide a calibration of the photospheric colors of young PMS stars as a function of their spectral type in a set of standard broad-band optical and near-infrared filters. For stars with masses of ~ 1.5Msun and ages of ~1-5 Myr, the chromospheric noise converts to a limit of measurable mass accretion rates of ~ 3x10^-10 Msun/yr.

The mass of a protoplanetary disk limits the formation and future growth of any planet. Masses of protoplanetary disks are usually calculated from measurements of the dust continuum emission by assuming an interstellar gas-to-dust ratio. To investigate the utility of CO as an alternate probe of disk mass, we use ALMA to survey $^{13}$CO and C$^{18}$O J = $3-2$ line emission from a sample of 93 protoplanetary disks around stars and brown dwarfs with masses from 0.03 -- 2 M$_{\odot}$ in the nearby Chamaeleon I star-forming region. We detect $^{13}$CO emission from 17 sources and C$^{18}$O from only one source. Gas masses for disks are then estimated by comparing the CO line luminosities to results from published disk models that include CO freeze-out and isotope-selective photodissociation. Under the assumption of a typical ISM CO-to-H$_2$ ratios of $10^{-4}$, the resulting gas masses are implausibly low, with an average gas mass of $\sim$ 0.05 M$_{Jup}$ as inferred from the average flux of stacked $^{13}$CO lines. The low gas masses and gas-to-dust ratios for Cha I disks are both consistent with similar results from disks in the Lupus star-forming region. The faint CO line emission may instead be explained if disks have much higher gas masses, but freeze-out of CO or complex C-bearing molecules is underestimated in disk models. The conversion of CO flux to CO gas mass also suffers from uncertainties in disk structures, which could affect gas temperatures. CO emission lines will only be a good tracer of the disk mass when models for C and CO depletion are confirmed to be accurate.

As part of the Accretion Discs in H$\alpha$ with OmegaCAM (ADHOC) survey, we imaged in r, i and H-alpha a region of 12x8 square degrees around the Orion Nebula Cluster. Thanks to the high-quality photometry obtained, we discovered three well-separated pre-main sequences in the color-magnitude diagram. The populations are all concentrated towards the cluster's center. Although several explanations can be invoked to explain these sequences we are left with two competitive, but intriguing, scenarios: a population of unresolved binaries with an exotic mass ratio distribution or three populations with different ages. Independent high-resolution spectroscopy supports the presence of discrete episodes of star formation, each separated by about a million years. The stars from the two putative youngest populations rotate faster than the older ones, in agreement with the evolution of stellar rotation observed in pre-main sequence stars younger than 4 Myr in several star forming regions. Whatever the final explanation, our results prompt for a revised look at the formation mode and early evolution of stars in clusters.

We present the first high-resolution sub-mm survey of both dust and gas for a large population of protoplanetary disks. Characterizing fundamental properties of protoplanetary disks on a statistical level is critical to understanding how disks evolve into the diverse exoplanet population. We use ALMA to survey 89 protoplanetary disks around stars with $M_{\ast}>0.1~M_{\odot}$ in the young (1--3~Myr), nearby (150--200~pc) Lupus complex. Our observations cover the 890~$\mu$m continuum and the $^{13}$CO and C$^{18}$O 3--2 lines. We use the sub-mm continuum to constrain $M_{\rm dust}$ to a few Martian masses (0.2--0.4~$M_{\oplus}$) and the CO isotopologue lines to constrain $M_{\rm gas}$ to roughly a Jupiter mass (assuming ISM-like $\rm {[CO]/[H_2]}$ abundance). Of 89 sources, we detect 62 in continuum, 36 in $^{13}$CO, and 11 in C$^{18}$O at $>3\sigma$ significance. Stacking individually undetected sources limits their average dust mass to $\lesssim6$ Lunar masses (0.03~$M_{\oplus}$), indicating rapid evolution once disk clearing begins. We find a positive correlation between $M_{\rm dust}$ and $M_{\ast}$, and present the first evidence for a positive correlation between $M_{\rm gas}$ and $M_{\ast}$, which may explain the dependence of giant planet frequency on host star mass. The mean dust mass in Lupus is 3$\times$ higher than in Upper Sco, while the dust mass distributions in Lupus and Taurus are statistically indistinguishable. Most detected disks have $M_{\rm gas}\lesssim1~M_{\rm Jup}$ and gas-to-dust ratios $<100$, assuming ISM-like $\rm {[CO]/[H_2]}$ abundance; unless CO is very depleted, the inferred gas depletion indicates that planet formation is well underway by a few Myr and may explain the unexpected prevalence of super-Earths in the exoplanet population.

The $\sigma$ Orionis cluster is important for studying protoplanetary disk evolution, as its intermediate age ($\sim$3-5 Myr) is comparable to the median disk lifetime. We use ALMA to conduct a high-sensitivity survey of dust and gas in 92 protoplanetary disks around $\sigma$ Orionis members with $M_{\ast}\gtrsim0.1 M_{\odot}$. Our observations cover the 1.33 mm continuum and several CO $J=2-1$ lines: out of 92 sources, we detect 37 in the mm continuum and six in $^{12}$CO, three in $^{13}$CO, and none in C$^{18}$O. Using the continuum emission to estimate dust mass, we find only 11 disks with $M_{\rm dust}\gtrsim10 M_{\oplus}$, indicating that after only a few Myr of evolution most disks lack sufficient dust to form giant planet cores. Stacking the individually undetected continuum sources limits their average dust mass to 5$\times$ lower than that of the faintest detected disk, supporting theoretical models that indicate rapid dissipation once disk clearing begins. Comparing the protoplanetary disk population in $\sigma$ Orionis to those of other star-forming regions supports the steady decline in average dust mass and the steepening of the $M_{\rm dust}$-$M_{\ast}$ relation with age; studying these evolutionary trends can inform the relative importance of different disk processes during key eras of planet formation. External photoevaporation from the central O9 star is influencing disk evolution throughout the region: dust masses clearly decline with decreasing separation from the photoionizing source, and the handful of CO detections exist at projected separations $>1.5$ pc. Collectively, our findings indicate that giant planet formation is inherently rare and/or well underway by a few Myr of age.

A homogeneous determination of basic stellar parameters of young stellar object (YSO) candidates is needed to confirm their evolutionary stage, membership to star forming regions (SFRs), and to get reliable values of the quantities related to chromospheric activity and accretion. We used the code ROTFIT and synthetic BT-Settl spectra for the determination of the atmospheric parameters (Teff and logg), the veiling, the radial (RV) and projected rotational velocity (vsini), from X-Shooter spectra of 102 YSO candidates in the Lupus SFR. We have shown that 13 candidates can be rejected as Lupus members based on their discrepant RV with respect to Lupus and/or the very low logg values. At least 11 of them are background giants. The spectral subtraction of inactive templates enabled us to measure the line fluxes for several diagnostics of both chromospheric activity and accretion. We found that all Class-III sources have H$\alpha$ fluxes compatible with a pure chromospheric activity, while objects with disks lie mostly above the boundary between chromospheres and accretion. YSOs with transitional disks displays both high and low H$\alpha$ fluxes. We found that the line fluxes per unit surface are tightly correlated with the accretion luminosity ($L_{\rm acc}$) derived from the Balmer continuum excess. This rules out that the relationships between $L_{\rm acc}$ and line luminosities found in previous works are simply due to calibration effects. We also found that the CaII-IRT flux ratio, $F_{8542}/F_{8498}$, is always small, indicating an optically thick emission source. The latter can be identified with the accretion shock near the stellar photosphere. The Balmer decrement reaches instead, for several accretors, high values typical of optically thin emission, suggesting that the Balmer emission originates in different parts of the accretion funnels with a smaller optical depth.

After a quiescence period of about 10 years, the classical EXor source V1118 Ori has undergone an accretion outburst in 2015 September. The maximum brightness (DV > 4 mag) was reached in 2015 December and was maintained for several months. Since 2016 September, the source is in a declining phase. Photometry and low/ high-resolution spectroscopy were obtained with MODS and LUCI2 at the {\it Large Binocular Telescope}, with the facilities at the Asiago 1.22 and 1.82 m telescopes, and with GIANO at the {\it Telescopio Nazionale Galileo}. The spectra are dominated by emission lines of \hi\ and neutral metallic species. From line and continuum analysis we derive the mass accretion rate and its evolution during the outburst. Considering that extinction may vary between 1.5 and 2.9 mag, we obtain m_acc= 0.3$-$2.0 10$^{-8}$ m_sun/yr, in quiescence and m_acc= 0.2$-$1.9 10$^{-6}$ m_sun/yr, at the outburst peak. The Balmer decrement shape has been interpreted by means of line excitation models, finding that from quiescence to outburst peak, the electron density has increased from $\sim$ 2 10$^9$ cm$^{-3}$ to $\sim$ 4 10$^{11}$ cm$^{-3}$. The profiles of the metallic lines are symmetric and narrower than 100 km s$^{-1}$, while \hi\, and \hei\,\,lines show prominent wings extending up to $\pm$ 500 km s$^{-1}$. The metallic lines likely originate at the base of the accretion columns, where neutrals are efficiently shielded against the ionizing photons, while faster ionized gas is closer to the star. Outflowing activity is testified by the detection of a variable P Cyg-like profile of the H$\alpha$ and \hei\, 1.08\,$\mu$m lines.

We have measured astrometry for members of the Orion Nebula Cluster with images obtained in 2015 with the Wide Field Camera 3 on board the Hubble Space Telescope. By comparing those data to previous measurements with NICMOS on Hubble in 1998, we have discovered that a star in the Kleinmann-Low Nebula, source x from Lonsdale et al. (1982), is moving with an unusually high proper motion of 29 mas/yr, which corresponds to 55 km/s at the distance of Orion. Previous radio observations have found that three other stars in the Kleinmann-Low Nebula (BN and sources I and n) have high proper motions (5-14 mas/yr) and were near a single location ~540 years ago, and thus may have been members of a multiple system that dynamically decayed. The proper motion of source x is consistent with ejection from that same location 540 years ago, which provides strong evidence that the dynamical decay did occur and that the runaway star BN originated in the Kleinmann-Low Nebula rather than the nearby Trapezium cluster. However, our constraint on the motion of source n is significantly smaller than the most recent radio measurement, which indicates that it did not participate in the event that ejected the other three stars.

We explore how measurements of protoplanetary disc masses and accretion rates provided by surveys of star forming regions can be analysed via the dimensionless accretion parameter, which we define as the product of the accretion rate and stellar age divided by the disc mass. By extending and generalising the study of Jones et al (2012), we demonstrate that this parameter should be less than or of order unity for a wide range of evolutionary scenarios, rising above unity only during the final stages of outside in clearing by external photoevaporation. We use this result to assess the reliability of disc mass estimates derived from CO isotopologues and submm continuum emission by examining the distribution of accretion efficiencies in regions which are not subject to external photoevaporation. We find that while dust based mass estimates produce results compatible with theoretical expectations assuming canonical dust/gas ratio, the systematically lower CO based estimates yield accretion efficiencies significantly above unity in contrast with the theory. This finding provides additional evidence that CO based disc masses are an under-estimate, in line with arguments that have been made on the basis of chemical modelling of relatively small samples. On the other hand, we demonstrate that dust based mass estimates are sufficiently accurate to reveal distinctly higher accretion efficiencies in the Trapezium cluster, where this result is expected given the evident importance of external photoevaporation. We therefore propose the dimensionless accretion parameter as a new diagnostic of external photoevaporation in other star forming regions.

We present the results of a study of the stellar and accretion properties of the (almost) complete sample of class II and transitional YSOs in the Lupus I, II, III and IV clouds, based on spectroscopic data acquired with the VLT/X-Shooter spectrograph. Our study combines the dataset from our previous work with new observations of 55 additional objects. We have investigated 92 YSO candidates in total, 11 of which have been definitely identified with giant stars unrelated to Lupus. The stellar and accretion properties of the 81 bona fide YSOs, which represent more than 90% of the whole class~II and transition disc YSO population in the aforementioned Lupus clouds, have been homogeneously and self-consistently derived, allowing for an unbiased study of accretion and its relationship with stellar parameters. The accretion luminosity, Lacc, increases with the stellar luminosity, Lstar, with an overall slope of ~1.6, similar but with a smaller scatter than in previous studies. There is a significant lack of strong accretors below Lstar~0.1Lsun, where Lacc is always lower than 0.01Lstar. We argue that the Lacc-Lstar slope is not due to observational biases, but is a true property of the Lupus YSOs. The logMacc-logMstar correlation shows a statistically significant evidence of a break, with a steeper relation for Mstar<0.2Msun and a flatter slope for higher masses. The bimodality of the Macc-Mstar relation is confirmed with four different evolutionary models used to derive the stellar mass. The bimodal behaviour of the observed relationship supports the importance of modelling self-gravity in the early evolution of the more massive discs, but other processes, such as photo evaporation and planet formation during the YSO's lifetime, may also lead to disc dispersal on different timescales depending on the stellar mass. We also refined the empirical Lacc vs. Lline relationships.

An era has started in which gas and dust can be observed independently in protoplanetary disks, thanks to the recent surveys with ALMA. The first near-complete high-resolution disk survey in both dust and gas in a single star-forming region has been carried out in Lupus, finding surprisingly low gas/dust ratios. The goal of this work is to fully exploit CO isotopologues observations in Lupus, comparing them with physical-chemical model results, in order to obtain gas masses for a large number of disks. We have employed physical-chemical models to analyze continuum and CO isotopologues observations of Lupus disks, including isotope-selective processes and freeze-out. Employing also the ALMA 13CO-only detections, disk gas masses have been calculated for a total of 34 sources, expanding the sample of 10 disks studied by Ansdell et al. (2016), where also C18O was detected. We confirm that overall gas-masses are very low, often smaller than 1 $M_{\rm J}$, if volatile carbon is not depleted. Accordingly, global gas/dust ratios predominantly between 1 and 10. Low CO-based gas masses and gas/dust ratios may indicate rapid loss of gas, or alternatively chemical evolution, e.g. via sequestering of carbon from CO to more complex molecules, or carbon locked up in larger bodies. Current ALMA observations cannot distinguish between these two hypotheses. We have simulated both scenarios, but chemical model results do not allow us to rule out one of the two. Assuming that all Lupus disks have evolved mainly due to viscous processes over the past few Myr, the observed correlation between the current mass accretion rate and dust mass found by Manara et al. (2016) implies a constant gas-to-dust ratio, which is close to 100 based on the observed $M_{\rm disk}/\dot{M}_{\rm acc}$ ratio. This in turn points to a scenario in which carbon depletion is responsible for the low CO isotopologue line luminosities.

In this work, we report the discovery of a candidate planetary-mass object with a photoevaporating protoplanetary disk, Proplyd 133-353, which is near the massive star $\theta^{1}$ Ori C at the center of the Orion Nebula Cluster (ONC). The object was known to have extended emission pointing away from $\theta^{1}$ Ori C, indicating ongoing external photoevaporation. Our near-infrared spectroscopic data suggests that the central source of Proplyd 133-353 is substellar ($\sim$M9.5), might have a mass probably less than 13 Jupiter mass and an age younger than 0.5 Myr. Proplyd 133-353 shows a similar ratio of X-ray luminosity to stellar luminosity to other young stars in the ONC with a similar stellar luminosity, and has a similar proper motion to the mean one of confirmed ONC members. We propose that Proplyd 133-353 was formed in a very low-mass dusty cloud near $\theta^{1}$ Ori C as a second-generation of star formation, which can explain both its young age and the presence of its disk.

RW Aur is a young binary system showing strong signatures of a recent tidal encounter between the circumprimary disk and the secondary star. The primary star has recently undergone two major dimming events ($\Delta$mag $\approx$ 2 in V-band), whose origin is still under debate. To shed light on the mechanism leading to the dimming events, we study the extinction properties, accretion variability, and gas kinematics using absorption lines from the material obscuring star RW Aur A. We compare our moderate resolution X-Shooter spectra of the dim state of RW Aur A with other spectral observations. In particular, we analyse archival high resolution UVES spectra obtained during the bright state of the system, in order to track the evolution of the spectral properties across the second dimming event. The spectrum obtained during the dim state shows narrow absorption lines in the Na and K optical doublets, where the former is saturated. With a velocity of -60 km/s these lines indicate that during the dim state the disk wind is either enhanced, or significantly displaced into the line of sight. The photometric evolution across the dimming event shows a gray extinction, and is correlated with a significant reduction of the EW of all photospheric lines. Emission lines tracing accretion do not vary significantly across the dimming. We conclude that the dimming event is related to a major perturbation on the inner disk. We suggest that the inner disk is occulting (most of) the star, and thus its photosphere, but is not occulting the accretion regions within a few stellar radii. Since observations of the outer disk indicate that the disk is modestly inclined (45 - 60 deg), we propose that the inner disk might be warped by a yet unseen (sub-)stellar companion, which may also explain the 2.77 day periodic variability of the spectral lines.

EPIC 204278916 has been serendipitously discovered from its K2 light curve which displays irregular dimmings of up to 65% for ~25 consecutive days out of 78.8 days of observations. For the remaining duration of the observations, the variability is highly periodic and attributed to stellar rotation. The star is a young, low-mass (M-type) pre-main-sequence star with clear evidence of a resolved tilted disk from ALMA observations. We examine the K2 light curve in detail and hypothesise that the irregular dimmings are caused by either a warped inner-disk edge or transiting cometary-like objects in either circular or eccentric orbits. The explanations discussed here are particularly relevant for other recently discovered young objects with similar absorption dips.

The disk mass is among the most important input parameter for every planet formation model to determine the number and masses of the planets that can form. We present an ALMA 887micron survey of the disk population around objects from 2 to 0.03Msun in the nearby 2Myr-old Chamaeleon I star-forming region. We detect thermal dust emission from 66 out of 93 disks, spatially resolve 34 of them, and identify two disks with large dust cavities of about 45AU in radius. Assuming isothermal and optically thin emission, we convert the 887micron flux densities into dust disk masses, hereafter Mdust. We find that the Mdust-Mstar relation is steeper than linear with power law indices 1.3-1.9, where the range reflects two extremes of the possible relation between the average dust temperature and stellar luminosity. By re-analyzing all millimeter data available for nearby regions in a self-consistent way, we show that the 1-3 Myr-old regions of Taurus, Lupus, and Chamaeleon I share the same Mdust-Mstar relation, while the 10Myr-old Upper Sco association has a steeper relation. Theoretical models of grain growth, drift, and fragmentation reproduce this trend and suggest that disks are in the fragmentation-limited regime. In this regime millimeter grains will be located closer in around lower-mass stars, a prediction that can be tested with deeper and higher spatial resolution ALMA observations.

A relation between the mass accretion rate onto the central young star and the mass of the surrounding protoplanetary disk has long been theoretically predicted and observationally sought. For the first time, we have accurately and homogeneously determined the photospheric parameters, mass accretion rate, and disk mass for an essentially complete sample of young stars with disks in the Lupus clouds. Our work combines the results of surveys conducted with VLT/X-Shooter and ALMA. With this dataset we are able to test a basic prediction of viscous accretion theory, the existence of a linear relation between the mass accretion rate onto the central star and the total disk mass. We find a correlation between the mass accretion rate and the disk dust mass, with a ratio that is roughly consistent with the expected viscous timescale when assuming an interstellar medium (ISM) gas-to-dust ratio. This confirms that mass accretion rates are related to the properties of the outer disk. We find no correlation between mass accretion rates and the disk mass measured by CO isotopologues emission lines, possibly owing to the small number of measured disk gas masses. This suggests that the mm-sized dust mass better traces the total disk mass and that masses derived from CO may be underestimated, at least in some cases.

Given that a majority of stars form in multiple systems, in order to fully understand the star- and planet-formation processes we must seek to understand them in multiple stellar systems. With this in mind, we present an analysis of the enigmatic binary T-Tauri system VV Corona Australis, in which both components host discs, but only one is visible at optical wavelengths. We seek to understand the peculiarities of this system by searching for a model for the binary which explains all the available continuum observations of the system. We present new mid-infrared interferometry and near-infrared spectroscopy along with archival millimetre-wave observations, which resolve the binary at 1.3mm for the first time. We compute a grid of pre-main-sequence radiative transfer models and calculate their posterior probabilities given the observed spectral energy distributions and mid-infrared interferometric visibilities of the binary components, beginning with the assumption that the only differences between the two components are their inclination and position angles. Our best-fitting solution corresponds to a relatively low luminosity T-tauri binary, with each component's disc having a large scale height and viewed at moderate inclination ($\sim 50^\circ$), with the infrared companion inclined by $\sim 5^\circ$ degrees more than the primary. Comparing the results of our model to evolutionary models suggests stellar masses $\sim 1.7$M$_\odot$ and an age for the system of 3.5Myr, towards the upper end of previous estimates. Combining these results with accretion indicators from near-IR spectroscopy, we determine an accretion rate of 4.0$\times 10^{-8}$M$_\odot$ yr$^{-1}$ for the primary. We suggest that future observations of VV~CrA and similar systems should prioritise high angular resolution sub-mm and near-IR imaging of the discs and high resolution optical/NIR spectroscopy of the central stars.

Depending on whether a T Tauri star accretes material from its circumstellar disk or not, different X-ray emission properties can be found. The accretion shocks produce cool heating of the plasma, contributing to the soft X-ray emission from the star. Using X-ray data from the Chandra Orion Ultra-deep Project and accretion rates that were obtained with the Hubble Space Telescope/WFPC2 photometric measurements in the Orion Nebula Cluster, we studied the relation between the accretion processes and the X-ray emissions of a coherent sample of T Tauri sources in the region. We performed regression and correlation analyses of our sample of T Tauri stars between the X-ray parameters, stellar properties, and the accretion measurements. We find that a clear anti-correlation is present between the residual X-ray luminosity and the accretion rates in our samples in Orion that is consistent with that found on the XMM-Newton Extended Survey of the Taurus molecular cloud (XEST) study. We provide a catalog with X-ray luminosities (corrected from distance) and accretion measurements of an Orion Nebula Cluster (ONC) T Tauri stars sample. Although Orion and Taurus display strong differences in their properties (total gas and dust mass, star density, strong irradiation from massive stars), we find that a similar relation between the residual X-ray emission and accretion rate is present in the Taurus molecular cloud and in the accreting samples from the Orion Nebula Cluster. The spread in the data suggests dependencies of the accretion rates and the X-ray luminosities other than the stellar mass, but the similarity between Orion and Taurus hints at the environment not being one of them. The anti-correlation between the residual X-ray luminosity and mass accretion rate is inherent to the T Tauri stars in general, independent of their birthplace and environment, and intrinsic to early stellar evolution.

We present the analysis of 34 new VLT/X-Shooter spectra of young stellar objects in the Chamaeleon I star forming region, together with four more spectra of stars in Taurus and two in Chamaeleon II. The broad wavelength coverage and accurate flux calibration of our spectra allow us to estimate stellar and accretion parameters for our targets by fitting the photospheric and accretion continuum emission from the Balmer continuum down to 700 nm. The dependence of accretion with stellar properties for this sample is consistent with previous results from the literature. The accretion rates for transitional disks are consistent with those of full disks in the same region. The spread of mass accretion rates at any given stellar mass is found to be smaller than in many studies, but is larger than that derived in the Lupus clouds using similar data and techniques. Differences in the stellar mass range and in the environmental conditions between our sample and that of Lupus may account for the discrepancy in scatter between Chamaeleon I and Lupus. Complete samples in Chamaeleon I and Lupus are needed to determine whether the difference in scatter of accretion rates and the lack of evolutionary trends are robust to sample selection.

The binary system RW Aur consists of two classical T~Tauri stars (CTTSs). The primary recently underwent its second observed major dimming event ($\Delta V\,\sim2\,$mag). We present new, resolved Chandra X-ray and UKIRT near-IR (NIR) data as well as unresolved optical photometry obtained in the dim state to study the gas and dust content of the absorber causing the dimming. The X-ray data show that the absorbing column density increased from $N_H<0.1\times10^{22}\,$cm$^{-2}$ during the bright state to $\approx2\times10^{22}\,$cm$^{-2}$ in the dim state. The brightness ratio between dim and bright state at optical to NIR wavelengths shows only a moderate wavelength dependence and the NIR color-color diagram suggests no substantial reddening. Taken together, this indicates gray absorption by large grains ($\gtrsim1\,\mu$m) with a dust mass column density of $\gtrsim2\times10^{-4}\,$g$\,$cm$^{-2}$. Comparison with $N_H$ shows that an absorber responsible for the optical/NIR dimming and the X-ray absorption is compatible with the ISM's gas-to-dust ratio, i.e., that grains grow in the disk surface layers without largely altering the gas-to-dust ratio. Lastly, we discuss a scenario in which a common mechanism can explain the long-lasting dimming in RW Aur and recently in AA Tau.

This booklet contains a collection of contributions to the meeting of the JEts and Disks at INAF (JEDI) group, which took place at the Capodimonte Observatory during 9-10 April 2015. Scope of the meeting was to bring together the JEDI researchers of the Italian Istituto Nazionale di Astrofisica (INAF) working in the field of circumstellar disks and jets in young stars, to discuss together the different agents affecting the structure and the evolution of disks, namely accretion, jets and winds. More information on the JEDI group and its activities can be found at \texttt{http://www.oa-roma.inaf.it/irgroup/JEDI}.

We present new VLT/X-Shooter optical and NIR spectra of a sample of 17 candidate young low-mass stars and BDs in the rho-Ophiucus cluster. We derived SpT and Av for all the targets, and then we determined their physical parameters. All the objects but one have M*<0.6 Msun, and 8 have mass below or close to the hydrogen-burning limit. Using the intensity of various emission lines present in their spectra, we determined the Lacc and Macc for all the objects. When compared with previous works targeting the same sample, we find that, in general, these objects are not as strongly accreting as previously reported, and we suggest that the reason is our more accurate estimate of the photospheric parameters. We also compare our findings with recent works in other slightly older star-forming regions to investigate possible differences in the accretion properties, but we find that the accretion properties for our targets have the same dependence on the stellar and substellar parameters as in the other regions. This leads us to conclude that we do not find evidence for a different dependence of Macc with M* when comparing low-mass stars and BDs. Moreover, we find a similar small (1 dex) scatter in the Macc-M* relation as in some of our recent works in other star-forming regions, and no significant differences in Macc due to different ages or properties of the regions. The latter result suffers, however, from low statistics and sample selection biases in the current studies. The small scatter in the Macc-M* correlation confirms that Macc in the literature based on uncertain photospheric parameters and single accretion indicators, such as the Ha width, can lead to a scatter that is unphysically large. Our studies show that only broadband spectroscopic surveys coupled with a detailed analysis of the photospheric and accretion properties allows us to properly study the evolution of disk accretion rates.

Photoevaporation due to high-energy stellar photons is thought to be one of the main drivers of protoplanetary disk dispersal. The fully or partially ionized disk surface is expected to produce free-free continuum emission at centimeter (cm) wavelengths that can be routinely detected with interferometers such as the upgraded Very Large Array (VLA). We use deep (rms noise down to 8 $\mu$Jy beam$^{-1}$ in the field of view center) 3.5 cm maps of the nearby (130 pc) Corona Australis (CrA) star formation (SF) region to constrain disk photoevaporation models. We find that the radio emission from disk sources in CrA is surprisingly faint. Only 3 out of 10 sources within the field of view are detected, with flux densities of order $10^2$ $\mu$Jy. However, a significant fraction of their emission is non-thermal. Typical upper limits for non-detections are $3\sigma\sim 60~\mu$Jy beam$^{-1}$. Assuming analytic expressions for the free-free emission from extreme-UV (EUV) irradiation, we derive stringent upper limits to the ionizing photon luminosity impinging on the disk surface $\Phi_\mathrm{EUV}<1-4\times10^{41}$ s$^{-1}$. These limits constrain $\Phi_\mathrm{EUV}$ to the low end of the values needed by EUV-driven photoevaporation models to clear protoplanetary disks in the observed few Myr timescale. Therefore, at least in CrA, EUV-driven photoevaporation is not likely to be the main agent of disk dispersal. We also compare the observed X-ray luminosities $L_X$ of disk sources with models in which photoevaporation is driven by such photons. Although predictions are less specific than for the EUV case, most of the observed fluxes (upper limits) are roughly consistent with the (scaled) predictions. Deeper observations, as well as predictions spanning a wider parameter space, are needed to properly test X-ray driven photoevaporation.

Transitional disks (TDs) are thought to be a late evolutionary stage of protoplanetary disks with dust depleted inner regions. The mechanism responsible for this depletion is still under debate. To constrain the models it is mandatory to have a good understanding of the properties of the gas content of the inner disk. Using X-Shooter broad band -UV to NIR- medium resolution spectroscopy we derive the stellar, accretion, and wind properties of a sample of 22 TDs. The analysis of these properties allows us to put strong constraints on the gas content in a region very close to the star (<0.2 AU) which is not accessible with any other observational technique. We fit the spectra with a self-consistent procedure to derive simultaneously SpT,Av,and mass accretion rates (Macc) of the targets. From forbidden emission lines we derive the wind properties of the targets. Comparing our findings to values for cTTs, we find that Macc and wind properties of 80% of the TDs in our sample, which is strongly biased towards strongly accreting objects, are comparable to those of cTTs. Thus, there are (at least) some TDs with Macc compatible with those of cTTs, irrespective of the size of the dust inner hole.Only in 2 cases Macc are much lower, while the wind properties are similar. We do not see any strong trend of Macc with the size of the dust depleted cavity, nor with the presence of a dusty optically thick disk close to the star. In the TDs in our sample there is a gas rich inner disk with density similar to that of cTTs disks. At least for some TDs, the process responsible of the inner disk clearing should allow for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanism producing the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.

The temporal decay of mass accretion in young stars is a fundamental tracer of the early evolution of circumstellar disks. Through population syntheses, we study how correlated uncertainties between the estimated parameters of young stars (luminosity, temperature, mass, age) and mass accretion rates Mdot, as well as observational selection effects, can bias the temporal decay of mass accretion rates (Mdot t^-eta) inferred from a comparison of measured Mdot with isochronal ages in young stellar clusters. We find that the presence of realistic uncertainties reduces the measured value of eta by up to a factor of 3, leading to the inference of shallower decays than the true value. This suggests a much faster temporal decay of Mdot than generally assumed. When considering the minimum uncertainties in ages affecting the Orion Nebula Cluster, the observed value eta~1.4, typical of Galactic star forming regions, can only be reproduced if the real decay exponent is eta>4. This effect becomes more severe if one assumes that observational uncertainties are larger, as required by some fast star formation scenarios. Our analysis shows that while selection effects due to sample incompleteness do bias eta, they can not alter this main result and strengthen it in many cases. A remaining uncertainty in our work is that it applies to the most commonly used and simple relationship between Mdot, the accretion luminosity and the stellar parameters. We briefly explore how a more complex interplay between these quantities might change the results.

We analyse current measurements of accretion rates onto pre-main sequence stars as a function of stellar mass, and conclude that the steep dependance of accretion rates on stellar mass is real and not driven by selection/detection threshold, as has been previously feared. These conclusions are reached by means of statistical tests including a survival analysis which can account for upper limits. The power-law slope of the Mdot-M* relation is found to be in the range of 1.6-1.9 for young stars with masses lower than 1 Msun. The measured slopes and distributions can be easily reproduced by means of a simple disc model which includes viscous accretion and X-ray photoevaporation. We conclude that the Mdot-M* relation in pre-main sequence stars bears the signature of disc dispersal by X-ray photoevaporation, suggesting that the relation is a straight- forward consequence of disc physics rather than an imprint of initial conditions.

We report JVLA 8-10 GHz ($\lambda$=3.0-3.7 cm) monitoring observations toward the YSO cluster R Coronae Australis (R\,CrA), taken in 2012, from March 15 to September 12. These observations were planned to measure the radio flux variabilities in timescales from 0.5 hours to several days, to tens of days, and up to $\sim$200 days. We found that among the YSOs detectable in individual epochs, in general, the most reddened objects in the \textit{Spitzer} observations show the highest mean 3.5 cm Stokes \textit{I} emission, and the lowest fractional variabilities on $<$200-day timescales. The brightest radio flux emitters in our observations are the two reddest sources IRS7W and IRS7E. In addition, by comparing with observations taken in 1996-1998 and 2005, we found that the radio fluxes of these two sources have increased by a factor $\sim$1.5. The mean 3.5-cm fluxes of the three Class I/II sources IRSI, IRS2, and IRS6 appear to be correlated with their accretion rates derived by a previous near infrared line survey. The weakly accreting Class I/II YSOs, or those in later evolutionary stages, present radio flux variability on $<$0.5-hour timescales. Some YSOs were detected only during occasional flaring events. The source R\,CrA went below our detection limit during a few fading events.

Young stars are known to show variability due to non-steady mass accretion rate from their circumstellar disks. Accretion flares can produce strong energetic irradiation and heating that may affect the disk in the planet formation region, close to the central star. During an extreme accretion outburst in the young star EX Lupi, the prototype of EXor variables, remarkable changes in molecular gas emission from $\sim1$ AU in the disk have recently been observed (Banzatti et al. 2012). Here, we focus on water vapor and explore how it is affected by variable accretion luminosity in T Tauri stars. We monitored a young highly variable solar-mass star, DR Tau, using simultaneously two high/medium-resolution ESO-VLT spectrographs: VISIR at 12.4 $\mu$m to observe water lines from the disk, and X-shooter covering from 0.3 to 2.5 $\mu$m to constrain the stellar accretion. Three epochs spanning timescales from several days to several weeks were obtained. Accretion luminosity was estimated to change within a factor $\sim2$, and no change in water emission was detected at a significant level. In comparison to EX Lupi and EXor outbursts, DR Tau suggests that the less long-lived and weaker variability phenomena typical of T Tauri stars may leave water at planet-forming radii in the disk mostly unaffected. We propose that these systems may provide evidence for two processes that act over different timescales: UV photochemistry in the disk atmosphere (faster) and heating of the disk deeper layers (slower).

We present X-Shooter/VLT observations of a sample of 36 accreting low-mass stellar and sub-stellar objects (YSOs) in the Lupus star forming region, spanning a range in mass from ~0.03 to ~1.2Msun, but mostly with 0.1Msun < Mstar < 0.5Msun. Our aim is twofold: firstly, analyse the relationship between excess-continuum and line emission accretion diagnostics, and, secondly, to investigate the accretion properties in terms of the physical properties of the central object. The accretion luminosity (Lacc), and from it the accretion rate (Macc), is derived by modelling the excess emission, from the UV to the near-IR, as the continuum emission of a slab of hydrogen. The flux and luminosity (Ll) of a large number of emission lines of H, He, CaII, etc., observed simultaneously in the range from ~330nm to 2500nm, were computed. The luminosity of all the lines is well correlated with Lacc. We provide empirical relationships between Lacc and the luminosity of 39 emission lines, which have a lower dispersion as compared to previous relationships in the literature. Our measurements extend the Pab and Brg relationships to Lacc values about two orders of magnitude lower than those reported in previous studies. We confirm that different methodologies to measure Lacc and Macc yield significantly different results: Ha line profile modelling may underestimate Macc by 0.6 to 0.8dex with respect to Macc derived from continuum-excess measures. Such differences may explain the likely spurious bi-modal relationships between Macc and other YSOs properties reported in the literature. We derive Macc in the range 2e-12 -- 4e-8 Msun/yr and conclude that Macc is proportional to Mstar^1.8(+/-0.2), with a dispersion lower by a factor of about 2 than in previous studies. A number of properties indicate that the physical conditions of the accreting gas are similar over more than 5 orders of magnitude in Macc.

We analyzed X-Shooter/VLT spectra of 24 ClassIII sources from three nearby star-forming regions (sigmaOrionis, LupusIII, and TWHya). We determined the effective temperature, surface gravity, rotational velocity, and radial velocity by comparing the observed spectra with synthetic BT-Settl model spectra. We investigated in detail the emission lines emerging from the stellar chromospheres and combined these data with archival X-ray data to allow for a comparison between chromospheric and coronal emissions. Both X-ray and Halpha luminosity as measured in terms of the bolometric luminosity are independent of the effective temperature for early-M stars but decline toward the end of the spectral M sequence. For the saturated early-M stars the average emission level is almost one dex higher for X-rays than for Halpha: log(L_x/L_bol) = -2.85 +- 0.36 vs. log(L_Halpha/L_bol) = -3.72 +- 0.21. When all chromospheric emission lines (including the Balmer series up to H11, CaII HK, the CaII infrared triplet, and several HeI lines) are summed up the coronal flux still dominates that of the chromosphere, typically by a factor 2-5. Flux-flux relations between activity diagnostics that probe different atmospheric layers (from the lower chromosphere to the corona) separate our sample of active pre-main sequence stars from the bulk of field M dwarfs studied in the literature. Flux ratios between individual optical emission lines show a smooth dependence on the effective temperature. The Balmer decrements can roughly be reproduced by an NLTE radiative transfer model devised for another young star of similar age. Future, more complete chromospheric model grids can be tested against this data set.

Current planet formation models are largely based on the observational constraint that protoplanetary disks have lifetime 3Myr. Recent studies, however, report the existence of PMS stars with signatures of accretion (strictly connected with the presence of circumstellar disks)and photometrically determined ages of 30 Myr, or more. Here we present a spectroscopic study of two major age outliers in the ONC. We use broad band, intermediate resolution VLT/X-Shooter spectra combined with an accurate method to determine the stellar parameters and the related age of the targets to confirm their peculiar age estimates and the presence of ongoing accretion.The analysis is based on a multi-component fitting technique, which derives simultaneously SpT, extinction, and accretion properties of the objects. With this method we confirm and quantify the ongoing accretion. From the photospheric parameters of the stars we derive their position on the HRD, and the age given by evolutionary models. Together with other age indicators like the lithium equivalent width we estimate with high accuracy the age of the objects. Our study shows that the two objects analyzed are not older than the typical population of the ONC. Our results show that, while photometric determination of the photospheric parameters are an accurate method to estimate the parameters of the bulk of young stellar populations, those of individual objects with high accretion rates and extinction may be affected by large uncertainties. Broad band spectroscopic determinations should thus be used to confirm the nature of individual objects. Our analysis shows that this method allows us to obtain an accurate determination of the photospheric parameters of accreting YSOs in any nearby star-forming region. We suggest that our detailed, broad- band spectroscopy method should be used to derive accurate properties of candidate old and accreting YSOs.

The Hubble Space Telescope (HST) Treasury Program on the Orion Nebula Cluster has used 104 orbits of HST time to image the Great Orion Nebula region with the Advanced Camera for Surveys (ACS), the Wide-Field/Planetary Camera 2 (WFPC2) and the Near Infrared Camera and Multi Object Spectrograph (NICMOS) instruments in 11 filters ranging from the U-band to the H-band equivalent of HST. The program has been intended to perform the definitive study of the stellar component of the ONC at visible wavelengths, addressing key questions like the cluster IMF, age spread, mass accretion, binarity and cirumstellar disk evolution. The scanning pattern allowed to cover a contiguous field of approximately 600 square arcminutes with both ACS and WFPC2, with a typical exposure time of approximately 11 minutes per ACS filter, corresponding to a point source depth AB(F435W) = 25.8 and AB(F775W)=25.2 with 0.2 magnitudes of photometric error. We describe the observations, data reduction and data products, including images, source catalogs and tools for quick look preview. In particular, we provide ACS photometry for 3399 stars, most of them detected at multiple epochs, WFPC2 photometry for 1643 stars, 1021 of them detected in the U-band, and NICMOS JH photometry for 2116 stars. We summarize the early science results that have been presented in a number of papers. The final set of images and the photometric catalogs are publicly available through the archive as High Level Science Products at the STScI Multimission Archive hosted by the Space Telescope Science Institute.

Context. The lack of knowledge of photospheric parameters and the level of chromospheric activity in young low-mass pre-main sequence stars introduces uncertainties when measuring mass accretion rates in accreting (Class II) Young Stellar Objects. A detailed investigation of the effect of chromospheric emission on the estimates of mass accretion rate in young low-mass stars is still missing. This can be undertaken using samples of young diskless (Class III) K and M-type stars. Aims. Our goal is to measure the chromospheric activity of Class III pre main sequence stars to determine its effect on the estimates of accretion luminosity (Lacc) and mass accretion rate (Macc) in young stellar objects with disks. Methods. Using VLT/X-Shooter spectra we have analyzed a sample of 24 non-accreting young stellar objects of spectral type between K5 and M9.5. We identify the main emission lines normally used as tracers of accretion in Class II objects, and we determine their fluxes in order to estimate the contribution of the chromospheric activity to the line luminosity. Results. We have used the relationships between line luminosity and accretion luminosity derived in the literature for Class II objects to evaluate the impact of chromospheric activity on the accretion rate measurements. We find that the typical chromospheric activity would bias the derived accretion luminosity by Lacc,noise< 10-3Lsun, with a strong dependence with the Teff of the objects. The noise on Macc depends on stellar mass and age, and the typical values of log(Macc,noise) range between -9.2 to -11.6Msun/yr. Conclusions. Values of Lacc< 10-3Lsun obtained in accreting low-mass pre main sequence stars through line luminosity should be treated with caution as the line emission may be dominated by the contribution of chromospheric activity.

The present observational understanding of the evolution of the mass accretion rates (Macc) in pre-main sequence stars is limited by the lack of accurate measurements of Macc over homogeneous and large statistical samples of young stars. Such observational effort is needed to properly constrain the theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study of Macc for a sample of \sim 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning limit to M\ast \sim 2M\odot. We derive Macc from both the U-band excess and the H{\alpha} luminosity (LH{\alpha}), after determining empirically both the shape of the typical accretion spectrum across the Balmer jump and the relation between the accretion luminosity (Lacc) and LH{\alpha}, that is Lacc/L\odot = (1.31\pm0.03)\cdotLH{\alpha}/L\odot + (2.63\pm 0.13). Given our large statistical sample, we are able to accurately investigate relations between Macc and the parameters of the central star such as mass and age. We clearly find Macc to increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that the decay of Macc with stellar age occurs over longer timescales for more massive PMS stars. Our best fit relation between these parameters is given by: log(Macc/M\odot\cdotyr)=(-5.12 \pm 0.86) -(0.46 \pm 0.13) \cdot log(t/yr) -(5.75 \pm 1.47)\cdot log(M\ast/M\odot) + (1.17 \pm 0.23)\cdot log(t/yr) \cdot log(M\ast/M\odot). These results also suggest that the similarity solution model could be revised for sources with M\ast > 0.5M\odot. Finally, we do not find a clear trend indicating environmental effects on the accretion properties of the sources.