The recent ALMA DSHARP survey provided illuminating results on the diversity of substructures in planet forming disks. These substructures trace pebble-sized grains accumulated at local pressure maxima, possibly due to planet-disk interactions or other planet formation processes. DSHARP sources are heavily biased to large and massive disks that only represent the high (dust flux) tail end of the disk population. Thus it is unclear whether similar substructures and corresponding physical processes also occur in the majority of disks which are fainter and more compact. Here we explore the presence and characteristics of features in a compact disk around GQ Lup A, the effective radius of which is 1.5 to 10 times smaller than those of DSHARP disks. We present our analysis of ALMA 1.3mm continuum observations of the GQ Lup system. By fitting visibility profiles of the continuum emission, we find substructures including a gap at ~ 10 au. The compact disk around GQ Lup exhibits similar substructures to those in the DSHARP sample, suggesting that mechanisms of trapping pebble-sized grains are at work in small disks as well. Characteristics of the feature at ~ 10 au, if due to a hidden planet, are evidence of planet formation at Saturnian distances. Our results hint at a rich world of substructures to be identified within the common population of compact disks, and subsequently a population of solar system analogs within these disks. Such study is critical to understanding the formation mechanisms and planet populations in the majority of protoplanetary disks.

We present high resolution millimeter continuum ALMA observations of the disks around the T Tauri stars LkCa 15 and J1610. These disks host dust-depleted inner regions, possibly carved by massive planets, and are of prime interest to study the imprints of planet-disk interactions. While at moderate angular resolution they appear as a broad ring surrounding a cavity, the continuum emission resolves into multiple rings at a resolution of ~60$\times$40 mas (~7.5 au for LkCa 15, ~6 au for J1610) and ~$7\,\mu$Jy beam$^{-1}$ rms at 1.3 mm. In addition to a broad extended component, LkCa 15 and J1610 host 3 and 2 narrow rings, respectively, with two bright rings in LkCa 15 being radially resolved. The rings look marginally optically thick, with peak optical depths of ~0.5 (neglecting scattering), in agreement with high angular resolution observations of full disks. We perform hydrodynamical simulations with an embedded, sub-Jovian-mass planet and show that the observed multi-ringed substructure can be qualitatively explained as the outcome of the planet-disk interaction. We note however that the choice of the disk cooling timescale alone can significantly impact the resulting gas and dust distributions around the planet, leading to different numbers of rings and gaps and different spacings between them. We propose that the massive outer disk regions of transition disks are favorable places for planetesimals and possibly second generation planet formation of objects with a lower mass than the planets carving the inner cavity (typically few $M_{\rm Jup}$), and that the annular substructures observed in LkCa 15 and J1610 may be indicative of planetary core formation within dust-rich pressure traps. Current observations are compatible with other mechanisms being at the origin of the observed substructures, in particular with narrow rings generated at the edge of the CO and N$_2$ snowlines.

Unraveling the 3D physical structure, the temperature and density distribution, of protoplanetary discs is an essential step if we are to confront simulations of embedded planets or dynamical instabilities. In this paper we focus on Submillimeter Array observations of the edge-on source, Gomez's Hamburger, believed to host an over-density hypothesised to be a product of gravitational instability in the disc, GoHam b. We demonstrate that, by leveraging the well characterised rotation of a Keplerian disc to deproject observations of molecular lines in position-position-velocity space into disc-centric coordinates, we are able to map out the emission distribution in the (r, z) plane and (x, y, z) space. We show that 12CO traces an elevated layer of $z\,/\,r \sim 0.3$, while 13CO traces deeper in the disc at $z\,/\,r \lesssim 0.2$. We localize emission associated with GoHam b, finding it at deprojected radius of approximately 500 au and at polar angle of +\- 30 degrees from the disc major axis. At the spatial resolution of $\sim 1.5^{\prime\prime}$, GoHam b is spatially unresolved, with an upper limit to its radius of $<190$~au.

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $\alpha/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($\alpha/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $\alpha$ turbulent parameter ($\alpha \sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

The chemical composition of planets is inherited by the distribution of the various molecular species in the protoplanetary disk at thetime of their formation. As of today, only a handful of disks has been imaged in multiple spectral lines with high spatial resolution. As part of a small campaign devoted to the chemical characterization of disk-outflow sources in Taurus, we report on new ALMA Band 6 observations with 20 au resolution toward the embedded young star DG Tau B. Images of the continuum emission reveals a dust disk with rings and, putatively, a leading spiral arm. The disk, as well as the prominent outflow cavities, are detected in CO, H2CO, CS, and CN while they remain undetected in SO2, HDO, and CH3OH. From the absorption of the back-side outflow, we inferred that the disk emission is optically thick in the inner 50 au. This morphology explains why no line emission is detected from this inner region and poses some limitations toward the calculation of the dust mass and the characterization of the inner gaseous disk. The H2CO and CS emission from the inner 200 au is mostly from the disk and their morphology is very similar. The CN emission significantly differs from the other two molecules as it is observed only beyond 150 au. This ring-like morphology is consistent with previous observations and the predictions of thermochemical disk models. Finally, we constrained the disk-integrated column density of all molecules. In particular, we found that the CH3OH/H2CO ratio must be smaller than 2, making the methanol non-detection still consistent with the only such a ratio available from the literature (1.27 in TW Hya).

Protoplanetary disks around young stars are the birth sights of planetary systems like our own. Disks represent the gaseous dusty matter left after the formation of their central stars. The mass and luminosity of the star, initial disk mass and angular momentum, and gas viscosity govern disk evolution and accretion. Protoplanetary disks are the cosmic nurseries where microscopic dust grains grow into pebbles, planetesimals, and planets.

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $\alpha$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

Understanding the physical structure of the planet formation environment, the protoplanetary disk, is essential for the interpretation of high resolution observations of the dust and future observations of the magnetic field structure. Observations of multiple transitions of molecular species offers a unique view of the underlying physical structure through excitation analyses. Here we describe a new method to extract high-resolution spectra from low-resolution observations, then provide two case studies of how molecular excitation analyses were used to constrain the physical structure in TW Hya, the closest protoplanetary disk to Earth.

Scattered light high-resolution imaging of the proto-planetary disc orbiting HD100453 shows two symmetric spiral arms, possibly launched by an external stellar companion. In this paper we present new, sensitive high-resolution ($\sim$30 mas) Band 7 ALMA observations of this source. This is the first source where we find counterparts in the sub-mm continuum to both scattered light spirals. The CO J=3-2 emission line also shows two spiral arms; in this case they can be traced over a more extended radial range, indicating that the southern spiral arm connects to the companion position. This is clear evidence that the companion is responsible for launching the spirals. The pitch angle of the sub-millimeter continuum spirals ($\sim 6 ^{\circ}$) is lower than the one in scattered light ($\sim 16 ^{\circ}$). We show that hydrodynamical simulations of binary-disc interaction can account for the difference in pitch angle only if one takes into account that the midplane is colder than the upper layers of the disc, as expected for the case of externally irradiated discs.

Protoplanetary disks are known to posses a stunning variety of substructure in the distribution of their mm~sized grains, predominantly seen as rings and gaps (Andrews et al. 2018), which are frequently interpreted as due to the shepherding of large grains by either hidden, still-forming planets within the disk (Zhang et al. 2018) or (magneto-)hydrodynamic instabilities (Flock et al. 2015). The velocity structure of the gas offers a unique probe of both the underlying mechanisms driving the evolution of the disk, the presence of embedded planets and characterising the transportation of material within the disk, such as following planet-building material from volatile-rich regions to the chemically-inert midplane, or detailing the required removal of angular momentum. Here we present the radial profiles of the three velocity components of gas in upper disk layers in the disk of HD 163296 as traced by 12CO molecular emission. These velocities reveal significant flows from the disk surface towards the midplane of disk at the radial locations of gaps argued to be opened by embedded planets (Isella et al. 2016, 2018, Teague et al. 2018, Pinte et al. 2018), bearing striking resemblance to meridional flows, long predicted to occur during the early stages of planet formation (Szulagyi et al. 2014, Morbidelli et al. 2014, Fung & Chiang 2016, Dong et al. 2019). In addition, a persistent radial outflow is seen at the outer edge of the disk, potentially the base of a wind associated with previously detected extended emission (Klaassen et al. 2013).

We report the detection of spiral substructure in both the gas velocity and temperature structure of the disk around TW~Hya, suggestive of planet-disk interactions with an unseen planet. Perturbations from Keplerian rotation tracing out a spiral pattern are observed in the SE of the disk, while significant azimuthal perturbations in the gas temperature are seen in the outer disk, outside 90~au, extending the full azimuth of the disk. The deviation in velocity is either $\Delta v_{\phi} \, / \, v_{\rm kep} \sim 0.1$ or $\Delta v_{z} \, / \, v_{\rm kep} \sim 0.01$ depending on whether the perturbation is in the rotational or vertical direction, while radial perturbations can be ruled out. Deviations in the gas temperature are $\pm 4$ K about the azimuthally averaged profile, equivalent to deviations of $\Delta T_{\rm gas} \, / \, T_{\rm gas} \sim 0.05$. Assuming all three structures can be described by an Archimedean spiral, measurements of the pitch angles of both velocity and temperature spirals show a radially decreasing trend for all three, ranging from 9$^{\circ}$ at 70 au, dropping to 3$^{\circ}$ at 200 au. Such low pitch-angled spirals are not readily explained through the wake of an embedded planet in the location of previously reported at 94 au, but rather require a launching mechanism which results in much more tightly wound spirals. Molecular emission tracing distinct heights in the disk is required to accurately distinguish between spiral launching mechanisms.

While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$\alpha$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$\mu$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$\mu$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$\mu$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems.

We present the discovery of a spatially unresolved source of sub-millimeter continuum emission ($\lambda=855$ $\mu$m) associated with a young planet, PDS 70 c, recently detected in H$\alpha$ emission around the 5 Myr old T Tauri star PDS 70. We interpret the emission as originating from a dusty circumplanetary disk with a dust mass between $2\times10^{-3}$ and $4.2 \times 10^{-3}$ Earth masses. Assuming a standard gas-to-dust ratio of 100, the ratio between the total mass of the circumplanetary disk and the mass of the central planet would be between $10^{-4}-10^{-5}$. Furthermore, we report the discovery of another compact continuum source located $0.074''\pm0.013''$ South-West of a second known planet in this system, PDS 70 b, that was previously detected in near-infrared images. We speculate that the latter source might trace dust orbiting in proximity of the planet, but more sensitive observations are required to unveil its nature.

We report on ALMA continuum observations of the black hole X-ray binary A0620-00, at an X-ray luminosity nine orders of magnitude sub-Eddington. The system was significantly detected at 98 GHz (at $44 \pm 7~\mu{\rm Jy}$) and only marginally at 233 GHz ($20 \pm 8~\mu{\rm Jy}$), about 40 days later. These results suggest either an optically thin sub-mm synchrotron spectrum, or highly variable sub-mm jet emission on month timescales. Although the latter appears more likely, we note that, at the time of the ALMA observations, A0620-00 was in a somewhat less active optical-IR state than during all published multi-wavelength campaigns when a flat-spectrum, partially self-absorbed jet has been suggested to extend from the radio to the mid-IR regime. Either interpretation is viable in the context of an internal shock model, where the jet's spectral shape and variability are set by the power density spectrum of the shells' Lorentz factor fluctuations. While strictly simultaneous radio-mm-IR observations are necessary to draw definitive conclusions for A0620-00, the data presented here, in combination with recent radio and sub-mm results from higher luminosity systems, demonstrate that jets from black hole X-ray binaries exhibit a high level of variability - either in flux density or intrinsic spectral shape, or both - across a wide spectrum of Eddington ratios. This is not in contrast with expectations from an internal shock model, where lower jet power systems can be expected to exhibit larger fractional variability owing to an overall decrease in synchrotron absorption.

Line ratios for different transitions of the same molecule have long been used as a probe of gas temperature. Here we use ALMA observations of the N2H+ J~=~1-0 and J~=~4-3 lines in the protoplanetary disk around TW Hya to derive the temperature at which these lines emit. We find an averaged temperature of 39~K with a one sigma uncertainty of 2~K for the radial range 0.8-2'', significantly warmer than the expected midplane temperature beyond 0.5'' in this disk. We conclude that the N2H+ emission in TW Hya is not emitting from near the midplane, but rather from higher in the disk, in a region likely bounded by processes such as photodissociation or chemical reprocessing of CO and N2 rather than freeze out.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns.

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.

We present an improved method to measure the rotation curves for disks with non-axisymmetric brightness profiles initially published in Teague et al. (2018a). Application of this method to the well studied AS$~$209 system shows substantial deviations from Keplerian rotation of up to $\pm 5\%$. These deviations are most likely due to perturbations in the gas pressure profile, including a perturbation located at $\approx 250~$au and spanning up to $\approx 50~$au which is only detected kinematically. Modelling the required temperature and density profiles required to recover the observed rotation curve we demonstrate that the rings observed in $\mu$m scattered light are coincident with the pressure maxima, and are radially offset from the rings observed in mm continuum emission. This suggests that if rings in the NIR are due to sub-$\mu$m grains trapped in pressure maxima that there is a vertical dependence on the radius of the pressure minima.

Measuring the centroid of a spectral line is a common problem in astronomy. Many methods have been devised to overcome limitations due to either noise in the spectra or asymmetric profiles, the most common of which are the intensity weighted averages (first moment) or fits of analytical (typically Gaussian) profiles. If the spectral line can be considered a single component, we demonstrate that a simple quadratic fit to the pixel of maximum intensity and its two neighboring pixels provides a robust measure of the line centroid. This approach allows for a sub-velocity resolution precision on the line centroid, without be biases by noise or asymmetric features in the line profile and outperforming traditional methods in most situations.

Observations of multiple rotational transitions from a single molecule allow for unparalleled constraints on the physical conditions of the emitting region. We present an analysis of CS in TW~Hya using the $J=7-6$, $5-4$ and $3-2$ transitions imaged at $\sim 0.5^{\prime\prime}$ spatial resolution, resulting in a temperature and column density profile of the CS emission region extending out to 230~au, far beyond previous measurements. In addition, the 15~kHz resolution of the observations and the ability to directly estimate the temperature of the CS emitting gas, allow for one of the most sensitive searches for turbulent broadening in a disk to date. Limits of $v_{\rm turb} \lesssim 0.1 c_s$ can be placed across the entire radius of the disk. We are able to place strict limits of the local H$_2$ density due to the collisional excitations of the observed transitions. From these we find that a minimum disk mass of $3 \times 10^{-4}~M_{\rm sun}$ is required to be consistent with the CS excitation conditions and can uniquely constrain the gas surface density profile in the outer disk.

Context. Several sulfur-bearing molecules are observed in the interstellar medium and in comets, in strong contrast to protoplanetary disks where only CS, H$_2$CS and SO have been detected so far. Aims. We combine observations and chemical models to constrain the sulfur abundances and their sensitivity to physical and chemical conditions in the DM Tau protoplanetary disk. Methods. We obtained $0.5^{"}$ ALMA observations of DM Tau in Bands 4 and 6 in lines of CS, SO, SO$_2$, OCS, CCS, H$_2$CS and H$_2$S, achieving a $\sim 5$ mJy sensitivity. Using the non-LTE radiative transfer code RADEX and the forward-modeling tool DiskFit, disk-averaged CS column densities and upper limits for the other species were derived. Results. Only CS was detected with a derived column density of $\sim 2-6 \times 10^{12}$ cm$^{-2}$. We report a first tentative detection of SO$_2$ in DM Tau. The upper limits range between $\sim 10^{11}$ and $10^{14}$ cm$^{-2}$ for the other S-bearing species. The best-fit chemical model matching these values requires a gas-phase C/O ratio of > 1 at $r \sim 50-100$ au. With chemical modeling we demonstrate that sulfur-bearing species could be robust tracers of the gas-phase C/O ratio, surface reaction rates, grain size and UV intensities. Conclusions. The lack of detections of a variety of sulfur-bearing molecules in DM Tau other than CS implies a dearth of reactive sulfur in the gas phase, either through efficient freeze-out or because most of the elemental sulfur is in other large species, as found in comets. The inferred high CS/SO and CS/SO$_2$ ratios require a non-solar C/O gas-phase ratio of > 1, consistent with the recent observations of hydrocarbon rings in DM Tau. The stronger depletion of oxygen-bearing S-species compared to CS is likely linked to the low observed abundances of gaseous water in DM Tau and points to a removal mechanism of oxygen from the gas.

We present the first kinematical detection of embedded protoplanets within a protoplanetary disk. Using archival ALMA observations of HD 163296, we demonstrate a new technique to measure the rotation curves of CO isotopologue emission to sub-percent precision relative to the Keplerian rotation. These rotation curves betray substantial deviations caused by local perturbations in the radial pressure gradient, likely driven by gaps carved in the gas surface density by Jupiter-mass planets. Comparison with hydrodynamic simulations shows excellent agreement with the gas rotation profile when the disk surface density is perturbed by two Jupiter mass planets at 83 au and 137 au. As the rotation of the gas is dependent on the pressure of the total gas component, this method provides a unique probe of the gas surface density profile without incurring significant uncertainties due to gas-to-dust ratios or local chemical abundances which plague other methods. Future analyses combining both methods promise to provide the most accurate and robust measures of embedded planetary mass. Furthermore, this method provides a unique opportunity to explore wide-separation planets beyond the mm continuum edge and to trace the gas pressure profile essential in modelling grain evolution in disks.

Turbulence is a fundamental parameter in models of grain growth during the early stages of planet formation. As such, observational constraints on its magnitude are crucial. Here we self-consistently analyze ALMA CO(2-1), SMA CO(3-2), and SMA CO(6-5) observations of the disk around TW Hya and find an upper limit on the turbulent broadening of $<$0.08c$_s$ ($\alpha<$0.007 for $\alpha$ defined only within 2-3 pressure scale heights above the midplane), lower than the tentative detection previously found from an analysis of the CO(2-1) data. We examine in detail the challenges of image plane fitting vs directly fitting the visibilities, while also considering the role of the vertical temperature gradient, systematic uncertainty in the amplitude calibration, and assumptions about the CO abundance, as potential sources of the discrepancy in the turbulence measurements. These tests result in variations of the turbulence limit between $<$0.04c$_s$ and $<$0.13c$_s$, consistently lower than the 0.2-0.4c$_s$ found previously. Having ruled out numerous factors, we restrict the source of the discrepancy to our assumed coupling between temperature and density through hydrostatic equilibrium in the presence of a vertical temperature gradient and/or the confinement of CO to a thin molecular layer above the midplane, although further work is needed to quantify the influence of these prescriptions. Assumptions about hydrostatic equilibrium and the CO distribution are physically motivated, and may have a small influence on measuring the kinematics of the gas, but they become important when constraining small effects such as the strength of the turbulence within a protoplanetary disk.

The paper presents new high angular resolution ALMA 1.3 mm dust continuum observations of the protoplanetary system AS 209 in the Ophiuchus star forming region. The dust continuum emission is characterized by a main central core and two prominent rings at $r = 75\,$au and $r = 130\,$au intervaled by two gaps at at $r = 62\,$au and $r = 103\,$au. The two gaps have different widths and depths, with the inner one being narrower and shallower. We determined the surface density of the millimeter dust grains using the 3D radiative transfer disk code \textsc{dali}. According to our fiducial model the inner gap is partially filled with millimeter grains while the outer gap is largely devoid of dust. The inferred surface density is compared to 3D hydrodynamical simulations (FARGO-3D) of planet-disk interaction. The outer dust gap is consistent with the presence of a giant planet ($M_{\rm planet} \sim 0.8\,M_{\rm Staturn}$); the planet is responsible for the gap opening and for the pile-up of dust at the outer edge of the planet orbit. The simulations also show that the same planet can give origin to the inner gap at $r = 62\,$au. The relative position of the two dust gaps is close to the 2:1 resonance and we have investigated the possibility of a second planet inside the inner gap. The resulting surface density (including location, width and depth of the two dust gaps) are in agreement with the observations. The properties of the inner gap pose a strong constraint to the mass of the inner planet ($M_{\rm planet} < 0.1\,M_{\rm J}$). In both scenarios (single or pair of planets), the hydrodynamical simulations suggest a very low disk viscosity ($\alpha < 10^{-4}$). Given the young age of the system (0.5 - 1 Myr), this result implies that the formation of giant planets occurs on a timescale of $\lesssim$ 1\,Myr.

Planets are born in protostellar disks, which are now observed with enough resolution to address questions about internal gas flows. Candidates for driving the flows include magnetic forces, but ionization state estimates suggest much of the gas mass decouples from magnetic fields. Thus, hydrodynamical instabilities could play a major role. We investigate disk dynamics under conditions typical for a T Tauri system, using global 3D radiation hydrodynamics simulations with embedded particles and a resolution of 70 cells per scale height. Stellar irradiation heating is included with realistic dust opacities. The disk starts in joint radiative balance and hydrostatic equilibrium. The vertical shear instability (VSI) develops into turbulence that persists up to at least 1600 inner orbits (143 outer orbits). Turbulent speeds are a few percent of the local sound speed at the midplane, increasing to 20%, or 100 m/s, in the corona. These are consistent with recent upper limits on turbulent speeds from optically thin and thick molecular line observations of TW Hya and HD 163296. The predominantly vertical motions induced by the VSI efficiently lift particles upwards. Grains 0.1 and 1 mm in size achieve scale heights greater than expected in isotropic turbulence. We conclude that while kinematic constraints from molecular line emission do not directly discriminate between magnetic and nonmagnetic disk models, the small dust scale heights measured in HL Tau and HD 163296 favor turbulent magnetic models, which reach lower ratios of the vertical kinetic energy density to the accretion stress.

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions.

Context: High-mass stars have a high degree of multiplicity and most likely form via disk accretion processes. The detailed physics of the binary and disk formation are still poorly constrained. Methods: Using the VLA in its most extended configuration at ~24GHz toward the prototypical high-mass star-forming region NGC7538IRS1 has allowed us to study the NH3 and thermal CH3OH emission and absorption as well as the cm continuum emission at an unprecedented spatial resolution of 0.06"x0.05", corresponding to a linear resolution of ~150AU at a distance of 2.7kpc. Results: A comparison of these new cm continuum data with previous VLA observations from 23yrs ago reveals no recognizable proper motions. If the emission were caused by a protostellar jet, proper motion signatures should have been easily identified. In combination with the high spectral indices S~nu^{alpha} (alpha between 1 and 2), this allows us to conclude that the continuum emission is from two hypercompact HII regions separated in projection by about 430AU. The NH3 spectral line data reveal a common rotating envelope indicating a bound high-mass binary system. In addition to this, the thermal CH3OH data show two separate velocity gradients across the two hypercompact HII regions. This indicates two disk-like structures within the same rotating circumbinary envelope. Disk and envelope structures are inclined by ~33deg, which can be explained by initially varying angular momentum distributions within the natal, turbulent cloud. Conclusions: Studying high-mass star formation at sub-0.1" resolution allows us to isolate multiple sources as well as to separate circumbinary from disk-like rotating structures.

The recent detection of gas-phase methanol (CH$_3$OH) lines in the disc of TW Hya by Walsh et al. provided the first observational constraints on the complex O-bearing organic content in protoplanetary discs. The emission has a ring-like morphology, with a peak at $\sim 30-50$ au and an inferred column density of $\sim 3-6\times10^{12}$ cm$^{-2}$. A low CH$_3$OH fractional abundance of $\sim 0.3-4\times 10^{-11}$ (with respect to H$_2$) is derived, depending on the assumed vertical location of the CH$_3$OH molecular layer. In this study, we use a thermo-chemical model of the TW Hya disc, coupled with the ALCHEMIC gas-grain chemical model, assuming laboratory-motivated, fast diffusivities of the surface molecules to interpret the CH$_3$OH detection. Based on this disc model, we performed radiative transfer calculations with the LIME code and simulations of the observations with the CASA simulator. We found that our model allows to reproduce the observations well. The CH$_3$OH emission in our model appears as a ring with radius of $\sim60$ au. Synthetic and observed line flux densities are equal within the rms noise level of observations. The synthetic CH$_3$OH spectra calculated assuming local thermodynamic equilibrium (LTE) can differ by up to a factor of 3.5 from the non-LTE spectra. For the strongest lines, the differences between LTE and non-LTE flux densities are very small and practically negligible. Variations in the diffusivity of the surface molecules can lead to variations of the CH$_3$OH abundance and, therefore, line flux densities by an order of magnitude.

We present ALMA Cycle~2 observations at $0.5^{\prime\prime}$ resolution of TW Hya of CS $J=5-4$ emission. The radial profile of the integrated line emission displays oscillatory features outwards of $1.5^{\prime\prime}$ ($\approx 90$ au). A dip-like feature at $1.6^{\prime\prime}$ is coincident in location, depth and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modelling and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modelling performed for scattered light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at $3.1^{\prime\prime}$ is identified as a region of intensified desorption due to enhanced penetration of the interstellar FUV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices.

We present scattered light images of the TW Hya disk performed with SPHERE in PDI mode at 0.63, 0.79, 1.24 and 1.62 micron. We also present H2/H3-band ADI observations. Three distinct radial depressions in the polarized intensity distribution are seen, around 85, 21, and 6~au. The overall intensity distribution has a high degree of azimuthal symmetry; the disk is somewhat brighter than average towards the South and darker towards the North-West. The ADI observations yielded no signifiant detection of point sources in the disk. Our observations have a linear spatial resolution of 1 to 2au, similar to that of recent ALMA dust continuum observations. The sub-micron sized dust grains that dominate the light scattering in the disk surface are strongly coupled to the gas. We created a radiative transfer disk model with self-consistent temperature and vertical structure iteration and including grain size-dependent dust settling. This method may provide independent constraints on the gas distribution at higher spatial resolution than is feasible with ALMA gas line observations. We find that the gas surface density in the "gaps" is reduced by 50% to 80% relative to an unperturbed model. Should embedded planets be responsible for carving the gaps then their masses are at most a few 10 Mearth. The observed gaps are wider, with shallower flanks, than expected for planet-disk interaction with such low-mass planets. If forming planetary bodies have undergone collapse and are in the "detachted phase" then they may be directly observable with future facilities such as METIS at the E-ELT.

We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

To unambiguously diagnose the detailed physical mechanisms and the evolutionary status of high-mass star-forming regions (HMSFRs), we have performed 0.4" ($\sim1000$ AU) resolution observations towards NGC7538S and IRS1, using the Plateau de Bure Interferometer (PdBI). This paper presents a joint analysis of the 1.37 mm continuum emission and the line intensity of 15 molecular species. Assuming local thermal equilibrium, we derived molecular column densities and molecular abundances for each internal gas substructure which is spatially resolved. These derived quantities are compared with a suite of 1-D gas-grain models. NGC7538S is resolved into at least three dense gas condensations. Despite the comparable continuum intensity of these condensations, their differing molecular line emission is suggestive of an overall chemical evolutionary trend from the northeast to the southeast. Since these condensations are embedded within the same parent gas core, their differing chemical properties are the most likely due to the different warm-up histories, rather than the different dynamic timescales. Despite remaining spatially unresolved, in IRS1 we detect abundant complex organic molecules, indicating that IRS1 is the most chemically evolved hot molecular core presented here. We observe a continuum that is dominated by absorption features with at least three strong emission lines, potentially from $\rm CH_3OH$. The $\rm CH_3OH$ lines which are purely in emission have higher excitation than the ones showing purely absorption. Potential reasons for that difference are discussed. This is the first comprehensive comparison of observations of the high-mass cores NGC7538S and IRS1 and a chemical model. We have found that different chemical evolutionary stages can coexist in the same natal gas core. Our achievement illustrates the strength of chemical analysis for understanding HMSFRs.

We present resolved Plateau de Bure Array observations of DM Tau in lines of HCO+ (3-2), (1-0) and DCO+ (3-2). A power-law fitting approach allowed a derivation of column densities of these two molecules. A chemical inner hole of ~50 AU was found in both HCO+ and DCO+ with DCO+ emission extending to only 450 AU. An isotopic ratio of R_D = N(DCO+) / N(HCO+) was found to range from 0.1 at 50 AU and 0.2 at 450 AU. Chemical modeling allowed an exploration of the sensitivity of these molecular abundances to physical parameters out with temperature, finding that X-rays were the domination ionization source in the HCO+ molecular region and that R_D also is sensitive to the CO depletion. The ionization fraction, assuming a steady state system, was found to be x(e-) ~ 10$^{-7}$. Modeling suggests that HCO+ is the dominant charged molecule in the disk but its contribution to ionization fraction is dwarfed by atmoic ions such as C+, S+ and H+.