Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms.

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

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.

Transition disks with large inner dust cavities are thought to host massive companions. However, the disk structure inside the companion orbit and how material flows toward an actively accreting star remain unclear. We present a high resolution continuum study of inner disks in the cavities of 38 transition disks. Measurements of the dust mass from archival Atacama Large Millimeter/Submillimeter Array observations are combined with stellar properties and spectral energy distributions to assemble a detailed picture of the inner disk. An inner dust disk is detected in 18 of 38 disks in our sample. Of the 14 resolved disks, 9 are significantly misaligned with the outer disk. The near-infrared excess is uncorrelated with the mm dust mass of the inner disk. The size-luminosity correlation known for protoplanetary disks is recovered for the inner disks as well, consistent with radial drift. The inner disks are depleted in dust relative to the outer disk and their dust mass is uncorrelated with the accretion rates. This is interpreted as the result of radial drift and trapping by planets in a low $\alpha$ ($\sim 10^{-3}$) disk, or a failure of the $\alpha$-disk model to describe angular momentum transport and accretion. The only disk in our sample with confirmed planets in the gap, PDS 70, has an inner disk with a significantly larger radius and lower inferred gas-to-dust ratio than other disks in the sample. We hypothesize that these inner disk properties and the detection of planets are due to the gap having only been opened recently by young, actively accreting planets.

Our aim was to directly detect the thermal emission of the putative exoring system responsible for the complex deep transits observed in the light curve for the young Sco-Cen star 1SWASP J140747.93-394542.6 (V1400 Cen, hereafter J1407), confirming it as the occulter seen in May 2007, and to determine its orbital parameters with respect to the star. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the field centred on J1407 in the 340 GHz (Band 7) continuum in order to determine the flux and astrometric location of the ring system relative to the star. We used the VLT/NACO camera to observe the J1407 system in March 2019 and to search for the central planetary mass object at thermal infrared wavelengths. We detect no point source at the expected location of J1407, and derive an upper limit $3\sigma$ level of $57.6~\mu\rm{Jy}$. There is a point source detected at an angular separation consistent with the expected location for a free-floating ring system that occulted J1407 in May 2007, with a flux of $89~\mu\rm{Jy}$ consistent with optically thin dust surrounding a massive substellar companion. At 3.8 microns with the NACO camera, we detect the star J1407 but no other additional point sources within 1.3 arcseconds of the star, with a lower bound on the sensitivity of $6M_{Jup}$ at the location of the ALMA source, and down to $4M_{Jup}$ in the sky background limit. The ALMA upper limit at the location of J1407 implies that a hypothesised bound ring system is composed of dust smaller than $1\rm{~mm}$ in size, implying a young ring structure. The detected ALMA source has multiple interpretations, including: (i) it is an unbound substellar object surrounded by warm dust in Sco-Cen with an upper mass limit of $6M_{Jup}$, or (ii) it is a background galaxy.

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

We test the hypothesis that the disc cavity in the `transition disc' Oph IRS 48 is carved by an unseen binary companion. We use 3D dust-gas smoothed-particle hydrodynamics simulations to demonstrate that marginally coupled dust grains concentrate in the gas over-density that forms in in the cavity around a low binary mass ratio binary. This produces high contrast ratio dust asymmetries at the cavity edge similar to those observed in the disc around IRS 48 and other transition discs. This structure was previously assumed to be a vortex. However, we show that the observed velocity map of IRS 48 displays a peculiar asymmetry that is not predicted by the vortex hypothesis. We show the unusual kinematics are naturally explained by the non-Keplerian flow of gas in an eccentric circumbinary cavity. We further show that perturbations observed in the isovelocity curves of IRS 48 may be explained as the product of the dynamical interaction between the companion and the disc. The presence of a $\sim$0.4 M$_{\odot}$ companion at a $\sim$10 au separation can qualitatively explain these observations. High spatial resolution line and continuum imaging should be able to confirm this hypothesis.

ALMA surveys have shown that CO emission in protoplanetary disks is much fainter than expected. Accordingly, CO-based gas masses and gas/dust ratios are orders of magnitude lower than previously thought. This may be explained either as fast gas dispersal, or as chemical evolution and locking up of volatiles in larger bodies leading to the low observed CO fluxes. The latter processes lead to enhanced C/O ratios in the gas, which may be reflected in enhanced abundances of molecules like C2H. The goal of this work is to employ C2H observations to understand whether low CO fluxes are caused by volatile depletion, or by fast gas dissipation. We present ALMA Cycle 4 C2H observations of a subsample of nine sources in Lupus. The integrated C2H emission is determined and compared to previous CO observations and model predictions. Seven out of nine disks are detected in C2H, whose line emission is almost as bright as 13CO. All detections are significantly brighter than the typical sensitivity of the observations, hinting at a bimodal distribution of the C2H line intensities. When compared with DALI models, the observed C2H fluxes can be reproduced only if some level of volatile carbon and oxygen depletion is allowed and C/O>1 in the gas. Models with reduced gas/dust ratios fail instead to reproduce the observations. A steeper than linear correlation between C2H and CN emission line is found. This is linked to the fact that C2H emission lines are affected more strongly by C/O variations than CN lines. Ring-like structures are detected both in C2H and in continuum emission but, as for CN, they do not seem to be connected. Sz 71 shows ring shaped emission in both C2H and CN with the location of the peak intensity coinciding. Our new ALMA C2H observations favour volatile carbon and oxygen depletion rather than fast gas dispersal to explain the faint CO observations for most of the disks.

We present new 0.6 - 4 $\mu$m imaging of the SR 21 transition disk from Keck/NIRC2 and Magellan/MagAO. The protoplanetary disk around SR 21 has a large (~ 30 - 40 AU) clearing first inferred from its spectral energy distribution and later detected in sub-millimeter imaging. Both the gas and small dust grains are known to have a different morphology, with an inner truncation in CO at ~7 AU, and micron-sized dust detected within the millimeter clearing. Previous near-infrared imaging could not distinguish between an inner dust disk with a truncation at ~7 AU or one that extended to the sublimation radius. The imaging data presented here require an inner dust disk radius of a few AU, and complex structure such as a warp or spiral. We present a parametric warped disk model that can reproduce the observations. Reconciling the images with the spectral energy distribution gathered from the literature suggests grain growth to ~2 - 5 $\mu$m within the sub-millimeter clearing. The complex disk structure and possible grain growth can be connected to dynamical shaping by a giant-planet mass companion, a scenario supported by previous observational and theoretical studies.

The wide, deep cavities of transition disks are often believed to have been hollowed out by nascent planetary systems. PDS 70, a ${\sim}5$ Myr old transition disk system in which a multi-Jupiter-mass planet candidate at 22 au coexists with a ${\sim}30$ au gas and ${\sim}60$ au dust-continuum gap, provides a valuable case study for this hypothesis. Using the PEnGUIn hydrodynamics code, we simulate the orbital evolution and accretion of PDS 70b in its natal disk over the lifetime of the system. When the accreting planet reaches about 2.5 Jupiter masses, it spontaneously grows in eccentricity and consumes material from a wide swathe of the PDS 70 disk; radiative transfer post-processing with DALI shows that this accurately reproduces the observed gap profile. Our results demonstrate that super-Jupiter planets can single-handedly carve out transition disk cavities, and indicate that the high eccentricities measured for such giants may be a natural consequence of disk-planet interaction.

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.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

Planet formation is thought to begin with the growth of dust particles in protoplanetary disks from micrometer to millimeter and centimeter sizes. Dust growth is hindered by a number of growth barriers, according to dust evolution theory, while observational evidence indicates that somehow these barriers must have been overcome. The observational evidence of dust traps, in particular the Oph IRS 48 disk, with the Atacama Large Millimeter/submillimeter Array (ALMA) has changed our view of the dust growth process. In this article I review the history of dust trapping in models and observations.

Since the discovery of the multi-ring structure of the HL Tau disk, ALMA data suggest that the dust continuum emission of many, if not all, protoplanetary disks consists of rings and gaps, no matter their spectral type or age. The origin of these gaps so far remains unclear. We present a sample study of 16 disks with multiple ring-like structures in the continuum, using published ALMA archival data, to compare their morphologies and gap locations in a systematic way. The 16 targets range from early to late type stars, from <0.5 Myr to >10 Myr, from ~0.2 to 40 L_Sun and include both full and transitional disks with cleared inner dust cavities. Stellar ages are revised using new Gaia distances. Gap locations are derived using a simple radial fit to the intensity profiles. Using a radiative transfer model, the temperature profiles are computed. The gap radii generally do not correspond to the orbital radii of snow lines of the most common molecules. A snow line model can likely be discarded as a common origin of multi-ring systems. In addition, there are no systematic trends in the gap locations that could be related to resonances of planets. Finally, the outer radius of the disks decreases for the oldest disks in the sample, indicating that if multi-ring disks evolve in a similar way, outer dust rings either dissipate with the gas or grow into planetesimal belts.

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.

ALMA has revolutionized our view of protoplanetary disks, revealing structures such as gaps, rings and asymmetries that indicate dust trapping as an important mechanism in the planet formation process. However, the high resolution images have also shown that the optically thin assumption for millimeter continuum emission may not be valid and the low values of the spectral index may be related to optical depth rather than dust growth. Longer wavelength observations are essential to properly disentangle these effects. The high sensitivity and spatial resolution of the next-generation Very Large Array (ngVLA) will open up the possibilities to spatially resolve disk continuum emission at centimeter wavelengths and beyond, which allows the study of dust growth in disks in the optically thin regime and further constrain models of planet formation.

High resolution ALMA observations of protoplanetary disks have revealed that many, if not all primordial disks consist of ring-like dust structures. The origin of these dust rings remains unclear, but a common explanation is the presence of planetary companions that have cleared gaps along their orbit and trapped the dust at the gap edge. A signature of this scenario is a decrease of gas density inside these gaps. In recent work, Isella et al. 2016 derived drops in gas density consistent with Saturn-mass planets inside the gaps in the HD163296 disk through spatially resolved CO isotopologue observations. However, as CO abundance and temperature depends on a large range of factors, the interpretation of CO emission is non-trivial. We use the physical-chemical code DALI to show that the gas temperature increases inside dust density gaps, implying that any gaps in the gas, if present, would have to be much deeper, consistent with planet masses higher than a Jupiter mass. Furthermore, we show that a model with increased grain growth at certain radii, as expected at a snowline, can reproduce the dust rings in HD163296 equally well without the need for companions. This scenario can explain both younger and older disks with observed gaps, as gaps have been seen in systems as young <1 Myr. While the origin of the rings in HD163296 remains unclear, these modeling results demonstrate that care has to be taken when interpreting CO emission in protoplanetary disk observations.

Spiral arms, rings and large scale asymmetries are structures observed in high resolution observations of protoplanetary disks, and it appears that some of the disks showing spiral arms in scattered light also show asymmetries in millimeter-sized dust. HD~135344B is one of these disks. Planets are invoked as the origin of these structures, but no planet has been observed so far. We investigate the nature of the asymmetric structure in the HD~135344B disk in order to understand the origin of the spirals and of the asymmetry seen in this disk. Ultimately, we aim at understanding whether or not one or more planets are needed to explain such structures. We present new ALMA sub-0.1" resolution observations in Band 3 and 4. The high spatial resolution allows us to characterize the mm-dust morphology of the disk. The low optical depth of continuum emission probes the bulk of the dust in vortex. Moreover, we combine the new observations with archival data to perform a multi-wavelength analysis and to obtain information about the dust distribution and properties inside the asymmetry. We resolve the asymmetric disk into a symmetric ring + asymmetric crescent, and observe that: (1) the spectral index strongly decreases at the center of the vortex, consistent with the presence of large grains; (2) for the first time, an azimuthal shift of the peak of the vortex with wavelength is observed; (3) the azimuthal width of the vortex decreases at longer wavelengths, as expected for dust traps. These features allow to confirm the nature of the asymmetry as a vortex. Finally a lower limit to the total mass of the vortex is $0.3 M_{\rm Jupiter}$. Considering the uncertainties involved in this estimate, it is possible that the actual mass of the vortex is higher and possibly within the required values ($\sim 4\,\rm M_{\rm Jupiter}$) to launch spiral arms similar to those observed in scattered light.

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.

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.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

We present a characterization of the protostar embedded within the BHR7 dark cloud, based on both photometric measurements from the near-infrared to millimeter and interferometric continuum and molecular line observations at millimeter wavelengths. We find that this protostar is a Class 0 system, the youngest class of protostars, measuring its bolometric temperature to be 50.5~K, with a bolometric luminosity of 9.3~L$_{\odot}$. The near-infrared and \textit{Spitzer} imaging show a prominent dark lane from dust extinction separating clear bipolar outflow cavities. Observations of $^{13}$CO ($J=2\rightarrow1$), C$^{18}$O ($J=2\rightarrow1$), and other molecular lines with the Submillimeter Array (SMA) exhibit a clear rotation signature on scales $<$1300~AU. The rotation can be traced to an inner radius of $\sim$170~AU and the rotation curve is consistent with an R$^{-1}$ profile, implying that angular momentum is being conserved. Observations of the 1.3~mm dust continuum with the SMA reveal a resolved continuum source, extended in the direction of the dark lane, orthogonal to the outflow. The deconvolved size of the continuum indicates a radius of $\sim$100~AU for the continuum source at the assumed distance of 400~pc. The visibility amplitude profile of the continuum emission cannot be reproduced by an envelope alone and needs a compact component. Thus, we posit that the resolved continuum source could be tracing a Keplerian disk in this very young system. If we assume that the continuum radius traces a Keplerian disk (R$\sim$120~AU) the observed rotation profile is consistent with a protostar mass of 1.0~$M_{\odot}$.

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.

We used the Robo-AO laser adaptive optics system to image 99 main sequence and subgiant stars that have Kepler-detected asteroseismic signals. Robo-AO allows us to resolve blended secondary sources at separations as close as 0.15" that may contribute to the measured Kepler light curves and affect asteroseismic analysis and interpretation. We report 8 new secondary sources within 4.0" of these Kepler asteroseismic stars. We used Subaru and Keck adaptive optics to measure differential infrared photometry for these candidate companion systems. Two of the secondary sources are likely foreground objects and at least 6 of the secondaries are background sources; however we cannot exclude the possibility that three of the objects may be physically associated. We measured a range of i'-band amplitude dilutions for the candidate companion systems from 0.43% to 15.4%. We find that the measured amplitude dilutions are insufficient to explain the previously identified excess scatter in the relationship between asteroseismic oscillation amplitude and the frequency of maximum power.

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.

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.

We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the $^{13}$CO and C$^{18}$O (J=2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is $\lesssim 0.12\,M_{\bigoplus}$, which gives a disk mass ratio in dust between the two components of $M_{\mathrm{dust, SR\,24S}}/M_{\mathrm{dust, SR\,24N}}\gtrsim840$. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to mm-sizes is inhibited there and that only warm gas, as seen by ro-vibrational CO emission inside the truncation radii of the binary, is present.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and $^{12}$CO, $^{13}$CO, and C$^{18}$O J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ~0".25 (35 AU). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-$\mu$m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius $\sim$70 AU that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 AU, or an edge that transitions smoothly over an annular width of 10 AU near 60 AU. In gas, the data are consistent with a cavity that is smaller, about 15 AU in radius, and whose surface density at 15 AU is $10^{3\pm1}$ times smaller than the surface density at 70 AU; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 AU or a double-drop model as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities, e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation.

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.

The recent progress in instrumentation and telescope development has brought us different ways to observe protoplanetary disks, including interferometers, space missions, adaptive optics, polarimetry, and time- and spectrally-resolved data. While the new facilities have changed the way we can tackle the existing open problems in disk structure and evolution, there is a substantial lack of interconnection between different observing techniques and their user communities. Here, we explore the complementarity of some of the state-of-the-art observing techniques, and how they can be brought together in a collective effort to understand how disks evolve and disperse at the time of planet formation. This paper was born at the "Protoplanetary Discussions" meeting in Edinburgh, 2016. Its goal is to clarify where multi-wavelength observations of disks converge in unveiling disk structure and evolution, and where they diverge and challenge our current understanding. We discuss caveats that should be considered when linking results from different observations, or when drawing conclusions based on limited datasets (in terms of wavelength or sample). We focus on disk properties that are currently being revolutionized by multi-wavelength observations. Specifically: the inner disk radius, holes and gaps and their link to large-scale disk structures, the disk mass, and the accretion rate. We discuss how the links between them, as well as the apparent contradictions, can help us to disentangle the disk physics and to learn about disk evolution.

In recent years spiral structures have been seen in scattered light observations and signs of vortices in millimeter images of protoplanetary disks, both probably linked with the presence of planets. We present ALMA Band 7 (335 GHz or 0.89 mm) continuum observations of the transition disk HD135344B at unprecedented spatial resolution of 0.16", using superuniform weighting. The data show that the asymmetric millimeter dust ring seen in previous work actually consists of an inner ring and an outer asymmetric structure. The outer feature is cospatial with the end of one of the spiral arms seen in scattered light, but the feature itself is not consistent with a spiral arm due to its coradiance. We propose a new possible scenario to explain the observed structures at both wavelengths. Hydrodynamical simulations show that a massive planet can generate a primary vortex (which dissipates at longer timescales, becoming an axisymmetric ring) and trigger the formation of a second generation vortex further out. Within this scenario the two spiral arms observed at scattered light originate from a planet at ~30 AU and from the secondary vortex at ~75 AU rather than a planet further out as previously reported.

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.

Understanding disk evolution and dissipation is essential for studies of planet formation. Transition disks, i.e., disks with large dust cavities and gaps, are promising candidates of active evolution. About two dozen SED-selected candidates have been confirmed to have dust cavities through millimeter interferometric imaging, but this sample is biased towards the brightest disks. The Spitzer surveys of nearby low-mass star forming regions have resulted in more than 4000 Young Stellar Objects (YSOs). Using color criteria we have selected a sample of ~150 candidates, and an additional 40 candidates and known transition disks from the literature. The Spitzer data were complemented by new observations at longer wavelengths, including new JCMT and APEX submillimeter photometry, and WISE and Herschel-PACS mid and far-infrared photometry. Furthermore, optical spectroscopy was obtained and stellar types were derived for 85% of the sample, including information from the literature. The SEDs were fit to a grid of RADMC-3D disk models with a limited number of parameters: disk mass, inner disk mass, scale height and flaring, and disk cavity radius, where the latter is the main parameter of interest. A large fraction of the targets possibly have dust cavities based on the SED. The derived cavity sizes are consistent with imaging/modeling results in the literature, where available. Trends are found with Ldisk/Lstar and stellar mass and a possible connection with exoplanet orbital radii. A comparison with a previous study where color observables are used (Cieza et al. 2010) reveals large overlap between their category of planet-forming disks and our transition disks with cavities. A large number of the new transition disk candidates are suitable for follow-up observations with ALMA.

Using Keck/NIRC2 $L^\prime$ (3.78 $\mu m$) data, we report the direct imaging discovery of a scattered light-resolved, solar system-scale residual protoplanetary disk around the young A-type star HD 141569A, interior to and concentric with the two ring-like structures at wider separations. The disk is resolved down to $\sim$ 0\farcs{}25 and appears as an arc-like rim with attached hook-like features. It is located at an angular separation intermediate between that of warm CO gas identified from spatially-resolved mid-infrared spectroscopy and diffuse dust emission recently discovered with the \textit{Hubble Space Telescope}. The inner disk has a radius of $\sim$ 39 AU, a position angle consistent with north-up, an inclination of $i$ $\sim$ 56$^{o}$, and has a center offset from the star. Forward-modeling of the disk favors a thick torus-like emission sharply truncated at separations beyond the torus' photocenter and heavily depleted at smaller separations. In particular, the best-fit density power law for the dust suggests that the inner disk dust and gas (as probed by CO) are radially segregated, a feature consistent with the dust trapping mechanism inferred from observations of "canonical" transitional disks. However, the inner disk component may instead be explained by radiation pressure-induced migration in optically-thin conditions, in contrast to the two stellar companion/planet-influenced ring-like structures at wider separations. HD 141569A's circumstellar environment --- with three nested, gapped, concentric dust populations --- is an excellent laboratory for understanding the relationship between planet formation and the evolution of both dust grains and disk architecture.

The gas-solid budget of carbon in protoplanetary disks is related to the composition of the cores and atmospheres of the planets forming in them. The key gas-phase carbon carriers CO, C$^{0}$ and C$^{+}$ can now be observed in disks. The gas-phase carbon abundance in disks has not yet been well characterized, we aim to obtain new constraints on the [C]/[H] ratio in a sample of disks, and to get an overview of the strength of [CI] and warm CO emission. We carried out a survey of the CO$\,6$--$5$ and [CI]$\,1$--$0$ and $2$--$1$ lines towards $37$ disks with APEX, and supplemented it with [CII] data from the literature. The data are interpreted using a grid of models produced with the DALI code. We also investigate how well the gas-phase carbon abundance can be determined in light of parameter uncertainties. The CO$\,6$--$5$ line is detected in $13$ out of $33$ sources, the [CI]$\,1$--$0$ in $6$ out of $12$, and the [CI]$\,2$--$1$ in $1$ out of $33$. With deep integrations, the first unambiguous detections of [CI]~$1$--$0$ in disks are obtained, in TW~Hya and HD~100546. Gas-phase carbon abundance reductions of a factor $5$--$10$ or more can be identified robustly based on CO and [CI] detections. The atomic carbon detection in TW~Hya confirms a factor $100$ reduction of [C]/[H]$_{\rm gas}$ in that disk, while the data are consistent with an ISM-like carbon abundance for HD~100546. In addition, BP~Tau, T~Cha, HD~139614, HD~141569, and HD~100453 are either carbon-depleted or gas-poor disks. The low [CI]~$2$--$1$ detection rates in the survey mostly reflect insufficient sensitivity to detect T~Tauri disks. The Herbig~Ae/Be disks with CO and [CII] upper limits below the models are debris disk like systems. A roughly order of magnitude increase in sensitivity compared to our survey is required to obtain useful constraints on the gas-phase [C]/[H] ratio in most of the targeted systems.

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336GHz of two transition disks, SR21 and HD135344B. In combination with previous ALMA observations from Cycle 0 at 689GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($\alpha_{\rm{mm}}$). The observations of SR21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR21 with values of $\alpha_{\rm{mm}}{\sim}3.8-4.2$ in the inner region ($r<35$ AU) and $\alpha_{\rm{mm}}{\sim}2.6-3.0$ outside. An axisymmetric ring (which we call the ring model) or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex model), is assumed for fitting the disk morphology. For SR21, the ring model better fits the emission at 336GHz, conversely the vortex model better fits the 689GHz emission. For HD135344B, neither a significant shift in the null of the visibilities nor radial variations of $\alpha_{\rm{mm}}$ are detected. Furthermore, for HD135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $\alpha_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD135344B, although it cannot be excluded either.

Azimuthally asymmetric dust distributions observed with ALMA in transition disks have been interpreted as dust traps. We present VLA Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 GHz or 0.45 mm) observations at 0.2" resolution of the Oph IRS 48 disk, which suggest that larger particles could be more azimuthally concentrated than smaller dust grains, assuming an axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 $\pm0.9$ times as wide as the centimeter emission, marginally consistent with the particle trapping mechanism under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized gas emission from the star. Both the morphology and the spectral index variations are consistent with an increase of large particles in the center of the trap, but uncertainties remain due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet formation models to allow dust particles to grow beyond millimeter sizes in the outer regions of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the dust trapping mechanism for centimeter-sized grains, although additional data is required for definitive confirmation.

We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

ALMA observations of a small sample of transitional disks with large dust cavities observed in Cycle 0 and 1 are summarized. The gas and dust surface density structures are inferred from the continuum and 12CO, 13CO and C18O line data using the DALI physical-chemical code. Thanks to its ability to self-shield, CO can survive inside dust cavities in spite of being exposed to intense UV radiation and can thus be used as a probe of the gas structure. Modeling of the existing data shows that gas is present inside the dust cavities in all cases, but at a reduced level compared with the gas surface density profile of the outer disk. The gas density decrease inside the dust cavity radius by factors of up to 10^4 suggests clearing by one or more planetary-mass companions. The accompanying pressure bumps naturally lead to trapping of the mm-sized dust grains observed in the ALMA images.

Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$\mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

Simple molecules like H2CO and CH3OH in protoplanetary disks are the starting point for the production of more complex organic molecules. So far, the observed chemical complexity in disks has been limited due to freeze out of molecules onto grains in the bulk of the cold outer disk. Complex molecules can be studied more directly in transitional disks with large inner holes, as these have a higher potential of detection, through UV heating of the outer disk and the directly exposed midplane at the wall. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 9 (~680 GHz) line data of the transitional disk Oph IRS 48, previously shown to have a large dust trap, to search for complex molecules in regions where planetesimals are forming. We report the detection of the H2CO 9(1,8)-8(1,7) line at 674 GHz, which is spatially resolved as a semi-ring at ~60 AU radius centered south from the star. The inferred H2CO abundance is ~10^{-8} derived by combining a physical disk model of the source with a non-LTE excitation calculation. Upper limits for CH3OH lines in the same disk give an abundance ratio H2CO/CH3OH>0.3, which points to both ice formation and gas-phase routes playing a role in the H2CO production. Upper limits on the abundances of H13CO+, CN and several other molecules in the disk are also derived and found to be consistent with full chemical models. The detection of the H2CO line demonstrates the start of complex organic molecules in a planet-forming disk. Future ALMA observations should be able to push down the abundance detection limits of other molecules by 1-2 orders of magnitude and test chemical models of organic molecules in (transitional) disks.

(Abridged) Transition disks are recognized by the absence of emission of small dust grains inside a radius of up to several 10s of AUs. Due to the lack of angular resolution and sensitivity, the gas content of such dust holes has not yet been determined, but is of importance to constrain the mechanism leading to the dust holes. Transition disks are thought to currently undergo the process of dispersal, setting an end to the giant planet formation process. We present new high-resolution observations with the Atacama Large Millimeter/ submillimeter Array (ALMA) of gas lines towards the transition disk Oph IRS 48 previously shown to host a large dust trap. ALMA has detected the $J=6-5$ line of $^{12}$CO and C$^{17}$O around 690 GHz (434 $\mu$m) at a resolution of $\sim$0.25$''$ corresponding to $\sim$30 AU (FWHM). The observed gas lines are used to set constraints on the gas surface density profile. New models of the physical-chemical structure of gas and dust in Oph IRS 48 are developed to reproduce the CO line emission together with the spectral energy distribution (SED) and the VLT-VISIR 18.7 $\mu$m dust continuum images. Integrated intensity cuts and the total spectrum from models having different trial gas surface density profiles are compared to observations. Using the derived surface density profiles, predictions for other CO isotopologues are made, which can be tested by future ALMA observations of the object. The derived gas surface density profile points to the clearing of the cavity by one or more massive planet/companion rather than just photoevaporation or grain-growth.

The statistics of discovered exoplanets suggest that planets form efficiently. However, there are fundamental unsolved problems, such as excessive inward drift of particles in protoplanetary disks during planet formation. Recent theories invoke dust traps to overcome this problem. We report the detection of a dust trap in the disk around the star Oph IRS 48 using observations from the Atacama Large Millimeter/submillimeter Array (ALMA). The 0.44-millimeter-wavelength continuum map shows high-contrast crescent-shaped emission on one side of the star originating from millimeter-sized grains, whereas both the mid-infrared image (micrometer-sized dust) and the gas traced by the carbon monoxide 6-5 rotational line suggest rings centered on the star. The difference in distribution of big grains versus small grains/gas can be modeled with a vortex-shaped dust trap triggered by a companion.

The outflow force of molecular bipolar outflows is a key parameter in theories of young stellar feedback on their surroundings. The focus of many outflow studies is the correlation between the outflow force, bolometric luminosity and envelope mass. However, it is difficult to combine the results of different studies in large evolutionary plots over many orders of magnitude due to the range of data quality, analysis methods and corrections for observational effects such as opacity and inclination. We aim to determine the outflow force for a sample of low luminosity embedded sources. We will quantify the influence of the analysis method and the assumptions entering the calculation of the outflow force. We use the James Clerk Maxwell Telescope to map 12CO J=3-2 over 2'x2' regions around 16 Class I sources of a well-defined sample in Ophiuchus at 15" resolution. The outflow force is then calculated using seven different methods differing e.g. in the use of intensity-weighted emission and correction factors for inclination. The results from the analysis methods differ from each other by up to a factor of 6, whereas observational properties and choices in the analysis procedure affect the outflow force by up to a factor of 4. For the sample of Class I objects, bipolar outflows are detected around 13 sources including 5 new detections, where the three non-detections are confused by nearby outflows from other sources. When combining outflow forces from different studies, a scatter by up to a factor of 5 can be expected. Although the true outflow force remains unknown, the separation method (separate calculation of dynamical time and momentum) is least affected by the uncertain observational parameters. The correlations between outflow force, bolometric luminosity and envelope mass are further confirmed down to low luminosity sources.

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

Gas-phase complex organic molecules are commonly detected toward high-mass protostellar hot cores. Detections toward low-mass protostars and outflows are comparatively rare, and a larger sample is key to investigate how the chemistry responds to its environment. Guided by the prediction that complex organic molecules form in CH3OH-rich ices and thermally or non-thermally evaporate with CH3OH, we have identified three sight-lines in the Serpens core - SMM1, SMM4 and SMM4-W - which are likely to be rich in complex organics. Using the IRAM 30m telescope, narrow lines (FWHM of 1-2 km s-1) of CH3CHO and CH3OCH3 are detected toward all sources, HCOOCH3 toward SMM1 and SMM4-W, and C2H5OH not at all. Beam-averaged abundances of individual complex organics range between 0.6 and 10% with respect to CH3OH when the CH3OH rotational temperature is applied. The summed complex organic abundances also vary by an order of magnitude, with the richest chemistry toward the most luminous protostar SMM1. The range of abundances compare well with other beam-averaged observations of low-mass sources. Complex organic abundances are of the same order of magnitude toward low-mass protostars and high-mass hot cores, but HCOOCH3 is relatively more important toward low-mass protostars. This is consistent with a sequential ice photochemistry, dominated by CHO-containing products at low temperatures and early times.