Nearby Star-forming Regions Research Articles

Effects of Ringed Structures and Dust Size Growth on Millimeter Observations of Protoplanetary Disks

Abstract The growth of solids from submicron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relatively small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall submillimeter and millimeter flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity (∼10 m s−1) and a high dust surface density (∼10 g cm−2 in the ring) can have millimeter spectral indices as low as ∼2.0, consistent with millimeter observations of faint disks in nearby star-forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than ∼2.5. Future multiwavelength high-resolution observations of these low-spectral-index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth.

Dynamical evolution of star-forming regions: III. Unbound stars and predictions for Gaia

We use $N$-body simulations to probe the early phases of the dynamical evolution of star-forming regions and focus on mass and velocity distributions of unbound stars. In this parameter space study, we vary the initial virial ratio and degree of spatial and kinematic substructure and analyse the fraction of stars that become unbound in two different mass classes (above and below 8 M$_{\odot}$). We find that the fraction of unbound stars differs depending on the initial conditions. After 10 Myr, in initially highly subvirial, substructured simulations, the high-mass and lower-mass unbound fractions are similar at $\sim$23 per cent. In initially virialised, substructured simulations, we find only $\sim$16 per cent of all high-mass stars are unbound, whereas $\sim$37 per cent of all lower-mass stars are. The velocity distributions of unbound stars only show differences for extremely different initial conditions. The distributions are dominated by large numbers of lower-mass stars becoming unbound just above the escape velocity of $\sim$3 km s$^{-1}$ with unbound high-mass stars moving faster on average than lower-mass unbound stars. We see no high-mass runaway stars (velocity > 30 km s$^{-1}$) from any of our initial conditions and only an occasional lower-mass runaway star from initially subvirial/substructured simulations. In our simulations, we find a small number of lower-mass walkaway stars (with velocity 5-30 km s$^{-1}$) from all of our initial conditions. These walkaway stars should be observable around many nearby star-forming regions with Gaia.

A centrally concentrated sub-solar-mass starless core in the Taurus L1495 filamentary complex

Abstract The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka the Morita array) and the IRAM 30-m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H2 volume density of ∼106 cm−3, MC5-N, connected to a narrow (width ∼0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim {0.2\!-\!0.4}\, M_{\odot }$, which is an order of magnitude smaller than typical low-mass pre-stellar cores. Taking into account a typical core to star formation efficiency for pre-stellar cores (∼20%–40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N2H+ and N2D+ observations using the Nobeyama 45-m telescope have confirmed the high-deuterium fractionation (∼30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

Discovery of a New Class I Methanol Maser Transition at 266.8 GHz

Abstract We report the detection of a new class I methanol maser candidate from the 52−41 E transition (266.8 GHz). This methanol transition has been detected toward a nearby high-mass star-forming region G352.630-1.067 (distance ∼0.7 kpc), in Submillimeter Array (SMA) observations. The new candidate transition has a similar spatial distribution as the 42−31 E (218.4 GHz) and 8−1−70 E (229.7 GHz) transitions, which are known class I maser transitions. Thermal methanol emission in this source is confined to a central hot core, while the three class I maser transitions are detected in two additional regions. These two maser-only emission regions are clearly associated with shocked gas traced by 2 μm Ks-band and thermal v = 0, J = 5−4 SiO molecular emission. In contrast to the thermal methanol emission from the hot core, the three class I maser transitions show an positive trend in the rotation diagram for the two maser regions. Large velocity gradient modeling of the 266.8, 218.4, and 229.7 GHz transitions shows that the 266.8 GHz transition can be a maser for a wide range of conditions. The intensity ratios for the three methanol transitions detected in maser regions can be reproduced under conditions that are typical for class I methanol maser sites. These facts all support the hypothesis that the detected emission from the 266.8 GHz methanol (52−41 E) transition is masing.

Close Companions around Young Stars

Abstract Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from Apache Point Observatory Galactic Evolution Experiment-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double-lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of ∼0.05–1.5 M ⊙, we find 399 binaries, including both radial velocity (RV) variables and SB2s. The mass ratio distribution of SB2s is consistent with being uniform for q < 0.95 with an excess of twins for q > 0.95. The period distribution is consistent with what has been observed in close binaries (<10 au) in the evolved populations. Three systems are found to have q ∼ 0.1, with a companion located within the brown dwarf desert. There are no strong trends in the multiplicity fraction as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at Σ* ∼ 30 stars/pc−2 and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of ∼2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal-mass systems that is less effective in binaries with lower q.

Driving conditions of protostellar outflows in different star-forming environments

The evolution of collapsing clouds embedded in different star-forming environments is investigated using three-dimensional non-ideal magnetohydrodynamics simulations considering different cloud metallicities ($Z/\thinspace Z_\odot$ = 0, 10$^{-5}$, 10$^{-4}$, 10$^{-3}$, 10$^{-2}$, 10$^{-1}$ and 1) and ionisation strengths ($C_\zeta$=0, 0.01, 1 and 10, where $C_\zeta$ is a coefficient controlling the ionisation intensity and $C_\zeta=1$ corresponds to the ionisation strength of nearby star-forming regions). With all combinations of these considered values of $Z/\thinspace Z_\odot$ and $C_\zeta$, 28 different star-forming environments are prepared and simulated. The cloud evolution in each environment is calculated until the central density reaches $n\approx10^{16}\,{\rm cm}^{-3}$ just before protostar formation, and the outflow driving conditions are derived. An outflow appears when the (first) adiabatic core forms in a magnetically active region where the magnetic field is well coupled with the neutral gas. In cases where outflows are driven, their momentum fluxes are always comparable to the observations of nearby star-forming regions. Thus, these outflows should control the mass growth of the protostars as in the local universe. Roughly, an outflow appears when $Z/\thinspace Z_\odot>10^{-4}$ and $C_\zeta \ge 0.01$. It is expected that the transition of the star formation mode from massive stars to normal solar-type stars occurs when the cloud metallicity is enhanced to the range of $Z/\thinspace Z_\odot\approx 10^{-4}$--$10^{-3}$, above which relatively low-mass stars would preferentially appear as a result of strong mass ejection.

The Green Bank Ammonia Survey: A Virial Analysis of Gould Belt Clouds in Data Release 1

Abstract We perform a virial analysis of starless dense cores in three nearby star-forming regions: L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which is used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used Planck- and Herschel-based column density maps for external cloud weight pressure and Five College Radio Astronomy Observatory 13CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.

Distances to molecular clouds at high galactic latitudes based on Gaia DR2

We report the distances of molecular clouds at high Galactic latitudes (|b| > 10°) derived from parallax and G-band extinction (AG) measurements in the second Gaia data release, Gaia DR2. Aided by Bayesian analyses, we determined distances by identifying the breakpoint in the extinction AG toward molecular clouds and using the extinction AG of Gaia stars around molecular clouds to confirm the breakpoint. We used nearby star-forming regions, such as Orion, Taurus, Cepheus, and Perseus, whose distances are well known to examine the reliability of our method. By comparing with previous results, we found that the molecular cloud distances derived from this method are reliable. The systematic error in the distances is approximately 5%. In total, 52 molecular clouds have well-determined distances, most of which are at high Galactic latitudes, and we provide reliable distances for 13 molecular clouds for the first time.

The potential of combining MATISSE and ALMA observations: constraining the structure of the innermost region in protoplanetary discs

Context. In order to study the initial conditions of planet formation, it is crucial to obtain spatially resolved multi-wavelength observations of the innermost region of protoplanetary discs. Aims. We evaluate the advantage of combining observations with MATISSE/VLTI and ALMA to constrain the radial and vertical structure of the dust in the innermost region of circumstellar discs in nearby star-forming regions. Methods. Based on a disc model with a parameterized dust density distribution, we apply 3D radiative-transfer simulations to obtain ideal intensity maps. These are used to derive the corresponding wavelength-dependent visibilities we would obtain with MATISSE as well as ALMA maps simulated with CASA. Results. Within the considered parameter space, we find that constraining the dust density structure in the innermost 5 au around the central star is challenging with MATISSE alone, whereas ALMA observations with reasonable integration times allow us to derive significant constraints on the disc surface density. However, we find that the estimation of the different disc parameters can be considerably improved by combining MATISSE and ALMA observations. For example, combining a 30-min ALMA observation (at 310 GHz with an angular resolution of 0.03′′) for MATISSE observations in the L and M bands (with visibility accuracies of about 3%) allows the radial density slope and the dust surface density profile to be constrained to within Δα = 0.3 and Δ(α − β) = 0.15, respectively. For an accuracy of ~1% even the disc flaring can be constrained to within Δβ = 0.1. To constrain the scale height to within 5 au, M band accuracies of 0.8% are required. While ALMA is sensitive to the number of large dust grains settled to the disc midplane we find that the impact of the surface density distribution of the large grains on the observed quantities is small.

The Stellar Membership of the Taurus Star-forming Region*

Abstract The high-precision astrometry from the second data release of the Gaia mission has made it possible to greatly improve the census of members of nearby clusters and associations. I have applied the Gaia data to the Taurus star-forming region, refining the sample of known members and identifying candidates for undiscovered members. The resulting samples of members and candidates provide the best constraints to date on the distribution of ages and the initial mass function (IMF) in Taurus. Several studies over the last 30 years have proposed the existence of a population of older stars (≳10 Myr) that is associated with the Taurus clouds. The data from Gaia demonstrate that such a population does not exist. Meanwhile, previous IMF estimates for small fields surrounding the Taurus aggregates have exhibited a surplus of K7–M0 stars (0.7–0.8 M ⊙) relative to star-forming clusters such as IC 348 and the Orion Nebula Cluster. However, that difference disappears when the new census of the entire region is considered, which should be complete for spectral types earlier than M6–M7 at A J < 1. Thus, there is little variation in the stellar IMF across the 3–4 orders of magnitude in stellar density that are present in nearby star-forming regions. Finally, I note that the proper motions of two previously known members, KPNO 15 and 2MASS J04355209 + 2255039, indicate that they may have been ejected from the same location within the L1536 cloud ∼7200 years ago.

Erratum: Physical properties and H-ionizing-photon production rates of extreme nearby star-forming regions

Erratum: Physical properties and H-ionizing-photon production rates of extreme nearby star-forming regions

An ALMA Survey of Faint Disks in the Chamaeleon I Star-forming Region: Why Are Some Class II Disks so Faint?

Abstract 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 signal-to-noise (S/N) detections or non-detections in our Cycle 2 survey of the ∼2 Myr old Chamaeleon I star-forming region. With five 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 Pascucci et al. 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 ⊕ and three very faint disks (CHXR20, ISO91, and T51) with dust masses ∼0.5 M ⊕. 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.

Physical properties and H-ionizing-photon production rates of extreme nearby star-forming regions

Measurements of the galaxy UV luminosity function at z>6 suggest that young stars hosted in low-mass star-forming galaxies produced the bulk of hydrogen-ionizing photons necessary to reionize the intergalactic medium (IGM) by redshift z~6. Whether star-forming galaxies dominated cosmic reionization, however, also depends on their stellar populations and interstellar medium properties, which set, among other things, the production rate of H-ionizing photons, $\xi_\text{ion}^\star$, and the fraction of these escaping into the IGM. Given the difficulty of constraining with existing observatories the physical properties of z>6 galaxies, in this work we focus on a sample of ten nearby objects showing UV spectral features comparable to those observed at z>6. We use the new-generation Beagle tool to model the UV-to-optical photometry and UV/optical emission lines of these Local 'analogues' of high-redshift galaxies, finding that our relatively simple, yet fully self-consistent, physical model can successfully reproduce the different observables considered. Our galaxies span a broad range of metallicities and are characterised by high ionization parameters, low dust attenuation, and very young stellar populations. Through our analysis, we derive a novel diagnostic of the production rate of H-ionizing photons per unit UV luminosity, $\xi_\text{ion}^\star$, based on the equivalent width of the bright $[\text{OIII}] \lambda 4959,5007$ doublet, which does not require measurements of H-recombination lines. This new diagnostic can be used to estimate $\xi_\text{ion}^\star$ from future direct measurements of the $[\text{OIII}] \lambda 4959,5007$ line using JWST/NIRSpec (out to z~9.5), and by exploiting the contamination by $\text{H}\beta + [\text{OIII}] \lambda 4959,5007$ of photometric observations of distant galaxies, for instance from existing Spitzer/IRAC data and from future ones with JWST/NIRCam.

The Pan-STARRS1 Proper-motion Survey for Young Brown Dwarfs in Nearby Star-forming Regions. I. Taurus Discoveries and a Reddening-free Classification Method for Ultracool Dwarfs

Abstract We are conducting a proper-motion survey for young brown dwarfs in the Taurus-Auriga molecular cloud based on the Pan-STARRS1 3π Survey. Our search uses multi-band photometry and astrometry to select candidates, and is wider (370 deg2) and deeper (down to ≈3 M Jup) than previous searches. We present here our search methods and spectroscopic follow-up of our high-priority candidates. Since extinction complicates spectral classification, we have developed a new approach using low-resolution (R ≈ 100) near-infrared spectra to quantify reddening-free spectral types, extinctions, and gravity classifications for mid-M to late-L ultracool dwarfs (≲100–3 M Jup in Taurus). We have discovered 25 low-gravity (vl-g) and the first 11 intermediate-gravity (int-g) substellar (M6–L1) members of Taurus, constituting the largest single increase of Taurus brown dwarfs to date. We have also discovered 1 new Pleiades member and 13 new members of the Perseus OB2 association, including a candidate very wide separation (58 kau) binary. We homogeneously reclassify the spectral types and extinctions of all previously known Taurus brown dwarfs. Altogether our discoveries have thus far increased the substellar census in Taurus by ≈40% and added three more L-type members (≲5–10 M Jup). Most notably, our discoveries reveal an older (>10 Myr) low-mass population in Taurus, in accord with recent studies of the higher-mass stellar members. The mass function appears to differ between the younger and older Taurus populations, possibly due to incompleteness of the older stellar members or different star formation processes.

The unexpectedly large proportion of high-mass star-forming cores in a Galactic mini-starburst

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

Circumstellar disc lifetimes in numerous galactic young stellar clusters

Photometric detections of dust circumstellar disks around pre-main sequence (PMS) stars, coupled with estimates of stellar ages, provide constraints on the time available for planet formation. Most previous studies on disk longevity, starting with Haisch, Lada & Lada (2001), use star samples from PMS clusters but do not consider datasets with homogeneous photometric sensitivities and/or ages placed on a uniform timescale. Here we conduct the largest study to date of the longevity of inner dust disks using X-ray and 1--8 micrometre infrared photometry from the MYStIX and SFiNCs projects for 69 young clusters in 32 nearby star-forming regions with ages t<=5 Myr. Cluster ages are derived by combining the empirical AgeJX method with PMS evolutionary models, which treat dynamo-generated magnetic fields in different ways. Leveraging X-ray data to identify disk-free objects, we impose similar stellar mass sensitivity limits for disk-bearing and disk-free YSOs while extending the analysis to stellar masses as low as M~0.1 Mo. We find that the disk longevity estimates are strongly affected by the choice of PMS evolutionary model. Assuming a disk fraction of 100% at zero age, the inferred disk half-life changes significantly, from t1/2 ~ 1.3--2 Myr to t1/2 ~ 3.5 Myr when switching from non-magnetic to magnetic PMS models. In addition, we find no statistically significant evidence that disk fraction varies with stellar mass within the first few Myr of life for stars with masses <2 Mo, but our samples may not be complete for more massive stars. The effects of initial disk fraction and star-forming environment are also explored.

A Spatially Resolved Study of the GRB 020903 Host Galaxy

Abstract GRB 020903 is a long-duration gamma-ray burst with a host galaxy close enough and extended enough for spatially resolved observations, making it one of less than a dozen GRBs where such host studies are possible. GRB 020903 lies in a galaxy host complex that appears to consist of four interacting components. Here we present the results of spatially resolved spectroscopic observations of the GRB 020903 host. By taking observations at two different position angles, we were able to obtain optical spectra (3600–9000 Å) of multiple regions in the galaxy. We confirm redshifts for three regions of the host galaxy that match that of GRB 020903. We measure the metallicity of these regions, and find that the explosion site and the nearby star-forming regions both have comparable subsolar metallicities. We conclude that, in agreement with past spatially resolved studies of GRBs, the GRB explosion site is representative of the host galaxy as a whole rather than localized in a metal-poor region of the galaxy.

Emission-line diagnostics of nearby H ii regions including interacting binary populations

We present numerical models of the nebular emission from H II regions around young stellar populations over a range of compositions and ages. The synthetic stellar pop- ulations include both single stars and interacting binary stars. We compare these models to the observed emission lines of 254 HII regions of 13 nearby spiral galax- ies and 21 dwarf galaxies drawn from archival data. The models are created using the combination of the Binary Population and Spectral Synthesis (bpass) code with the photoionization code cloudy to study the differences caused by the inclusion of interacting binary stars in the stellar population. We obtain agreement with the observed emission line ratios from the nearby star-forming regions and discuss the effect of binary star evolution pathways on the nebular ionization of H II regions. We find that at population ages above 10 Myr, single-star models rapidly decrease in flux and ionization strength, while binary-star model still produce strong flux and high O III/H\b{eta} ratios. Our models can reproduce the metallicity of H II regions from spiral galaxies but we find higher metallicities than previously estimated for the H II regions from dwarf galaxies. Comparing the equivalent width of H\b{eta} emission between mod- els and observations, we find that accounting for ionizing photon leakage can affect age estimates for HII regions. When it is included, the typical age derived for HII regions is 5Myr from single-star models, and up to 10Myr with binary-star models. This is due to the existence of binary-star evolution pathways which produce more hot WR and helium stars at older ages. For future reference, we calculate new BPASS binary maximal starburst lines as a function of metallicity, and for the total model population, and present these in Appendix A.

Constraints on observing brightness asymmetries in protoplanetary disks at solar system scale

We have quantified the potential capabilities of detecting local brightness asymmetries in circumstellar disks with the Very Large Telescope Interferometer (VLTI) in the mid-infrared wavelength range. The study is motivated by the need to evaluate theoretical models of planet formation by direct observations of protoplanets at early evolutionary stages, when they are still embedded in their host disk. Up to now, only a few embedded candidate protoplanets have been detected with semi-major axes of 20–50 au. Due to the small angular separation from their central star, only long-baseline interferometry provides the angular resolving power to detect disk asymmetries associated to protoplanets at solar system scales in nearby star-forming regions. In particular, infrared observations are crucial to observe scattered stellar radiation and thermal re-emission in the vicinity of embedded companions directly. For this purpose we performed radiative transfer simulations to calculate the thermal re-emission and scattered stellar flux from a protoplanetary disk hosting an embedded companion. Based on that, visibilities and closure phases are calculated to simulate observations with the future beam combiner MATISSE, operating at the L, M and N bands at the VLTI. We find that the flux ratio of the embedded source to the central star can be as low as 0.5 to 0.6% for a detection at a feasible significance level due to the heated dust in the vicinity of the embedded source. Furthermore, we find that the likelihood for detection is highest for sources at intermediate distances r ≈ 2–5 au and disk masses not higher than ≈10−4 M⊙.

New Insights into the Nature of Transition Disks from a Complete Disk Survey of the Lupus Star-forming Region

Abstract 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 toward 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 to 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 , 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.