The volatile contents of protoplanetary disks both set the potential for planetary chemistry and provide valuable probes of defining disk system characteristics such as stellar mass, gas mass, ionization, and temperature structure. Current disk molecular inventories are fragmented, however, giving an incomplete picture: unbiased spectral line surveys are needed to assess the volatile content. We present here an overview of such a survey of the protoplanetary disks around the Herbig Ae star MWC 480 and the T Tauri star LkCa 15 in ALMA Band 7, spanning $\sim$36 GHz from 275--317 GHz and representing an order of magnitude increase in sensitivity over previous single-dish surveys. We detect 14 molecular species (including isotopologues), with 5 species (C$^{34}$S, $^{13}$CS, H$_{2}$CS, DNC, and C$_2$D) detected for the first time in protoplanetary disks. Significant differences are observed in the molecular inventories of MWC~480 and LkCa~15, and we discuss how these results may be interpreted in light of the different physical conditions of these two disk systems.

The volatile composition of a planet is determined by the inventory of gas and ice in the parent disk. The volatile chemistry in the disk is expected to evolve over time, though this evolution is poorly constrained observationally. We present ALMA observations of C18O, C2H, and the isotopologues H13CN, HC15N, and DCN towards five Class 0/I disk candidates. Combined with a sample of fourteen Class II disks presented in Bergner et al. (2019b), this data set offers a view of volatile chemical evolution over the disk lifetime. Our estimates of C18O abundances are consistent with a rapid depletion of CO in the first ~0.5-1 Myr of the disk lifetime. We do not see evidence that C2H and HCN formation are enhanced by CO depletion, possibly because the gas is already quite under-abundant in CO. Further CO depletion may actually hinder their production by limiting the gas-phase carbon supply. The embedded sources show several chemical differences compared to the Class II stage, which seem to arise from shielding of radiation by the envelope (impacting C2H formation and HC15N fractionation) and sublimation of ices from infalling material (impacting HCN and C18O abundances). Such chemical differences between Class 0/I and Class II sources may affect the volatile composition of planet-forming material at different stages in the disk lifetime.

Turbulence is a crucial factor in many models of planet formation, but it has only been directly constrained among a small number of planet forming disks. Building on the upper limits on turbulence placed in disks around HD 163296 and TW Hya, we present ALMA CO J=2-1 line observations at $\sim0.3"$ (20-50 au) resolution and 80 ms$^{-1}$ channel spacing of the disks around DM Tau, MWC 480, and V4046 Sgr. Using parametric models of disk structure, we robustly detect non-thermal gas motions around DM Tau of between 0.25 c$_s$ and 0.33 c$_s$, with the range dominated by systematic effects, making this one of the only systems with directly measured non-zero turbulence. Using the same methodology, we place stringent upper limits on the non-thermal gas motion around MWC 480 ($<$0.08c$_s$) and V4046 Sgr ($<$0.12 c$_s$). The preponderance of upper limits in this small sample, and the modest turbulence levels consistent with dust studies, suggest that weak turbulence ($\alpha\lesssim10^{-3}$) may be a common, albeit not universal, feature of planet-forming disks. We explore the particular physical conditions around DM Tau that could lead this system to be more turbulent than the others.

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.

H$_2$CO is one of the most abundant organic molecules in protoplanetary disks and can serve as a precursor to more complex organic chemistry. We present an ALMA survey of H$_2$CO towards 15 disks covering a range of stellar spectral types, stellar ages, and dust continuum morphologies. H$_2$CO is detected towards 13 disks and tentatively detected towards a 14th. We find both centrally-peaked and centrally-depressed emission morphologies, and half of the disks show ring-like structures at or beyond expected CO snowline locations. Together these morphologies suggest that H$_2$CO in disks is commonly produced through both gas-phase and CO-ice-regulated grain-surface chemistry. We extract disk-averaged and azimuthally-averaged H$_2$CO excitation temperatures and column densities for four disks with multiple H$_2$CO line detections. The temperatures are between 20-50K, with the exception of colder temperatures in the DM Tau disk. These temperatures suggest that H$_2$CO emission in disks is generally emerging from the warm molecular layer, with some contributions from the colder midplane. Applying the same H$_2$CO excitation temperatures to all disks in the survey, we find that H$_2$CO column densities span almost three orders of magnitude ($\sim 5 \times 10^{11} - 5 \times 10^{14} \mathrm{cm}^{-2}$). The column densities appear uncorrelated with disk size and stellar age, but Herbig Ae disks may have less H$_2$CO compared to T Tauri disks, possibly because of less CO freeze-out. More H$_2$CO observations towards Herbig Ae disks are needed to confirm this tentative trend, and to better constrain under which disk conditions H$_2$CO and other oxygen-bearing organics efficiently form during planet formation.

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."

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

Snowlines of major volatiles regulate the gas and solid C/N/O ratios in the planet-forming midplanes of protoplanetary disks. Snow surfaces are the 2D extensions of snowlines in the outer disk regions, where radiative heating results in a decreasing temperature with disk height. CO and N$_2$ are two of the most abundant carriers of C, N and O. N$_2$H$^+$ can be used to probe the snow surfaces of both molecules, because it is destroyed by CO and formed from N$_2$. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations of N$_2$H$^+$ at 0.2$''$-0.4$''$ resolution in the disks around LkCa 15, GM Aur, DM Tau, V4046 Sgr, AS 209, and IM Lup. We find two distinctive emission morphologies: N$_2$H$^+$ is either present in a bright, narrow ring surrounded by extended tenuous emission, or in a broad ring. These emission patterns can be explained by two different kinds of vertical temperature structures. Bright, narrow N$_2$H$^+$ rings are expected in disks with a thick Vertically Isothermal Region above the Midplane (VIRaM) layer (LkCa 15, GM Aur, DM Tau) where the N$_2$H$^+$ emission peaks between the CO and N$_2$ snowlines. Broad N$_2$H$^+$ rings come from disks with a thin VIRaM layer (V4046 Sgr, AS 209, IM Lup). We use a simple model to extract the first sets of CO and N$_2$ snowline pairs and corresponding freeze-out temperatures towards the disks with a thick VIRaM layer. The results reveal a range of N$_2$ and CO snowline radii towards stars of similar spectral type, demonstrating the need for empirically determined snowlines in disks.

Molecular lines observed towards protoplanetary disks carry information about physical and chemical processes associated with planet formation. We present ALMA Band 6 observations of C2H, HCN, and C18O in a sample of 14 disks spanning a range of ages, stellar luminosities, and stellar masses. Using C2H and HCN hyperfine structure fitting and HCN/H13CN isotopologue analysis, we extract optical depth, excitation temperature, and column density radial profiles for a subset of disks. C2H is marginally optically thick (\tau ~1-5) and HCN is quite optically thick (\tau ~ 5-10) in the inner 200 AU. The extracted temperatures of both molecules are low (10-30K), indicative of either sub-thermal emission from the warm disk atmosphere or substantial beam dilution due to chemical substructure. We explore the origins of C2H morphological diversity in our sample using a series of toy disk models, and find that disk-dependent overlap between regions with high UV fluxes and high atomic carbon abundances can explain a wide range of C2H emission features (e.g. compact vs. extended and ringed vs. ringless emission). We explore the chemical relationship between C2H, HCN, and C18O and find a positive correlation between C2H and HCN fluxes, but no relationship between C2H or HCN with C18O fluxes. We also see no evidence that C2H and HCN are enhanced with disk age. C2H and HCN seem to share a common driver, however more work remains to elucidate the chemical relationship between these molecules and the underlying evolution of C, N, and O chemistries in disks.

Methanol (CH$_3$OH) is at the root of organic ice chemistry in protoplanetary disks. However, its weak emission has made detections difficult. To date, gas-phase CH$_3$OH has been detected in only one Class II disk, TW Hya. We use the Atacama Large Millimeter/submillimeter Array (ALMA) to search for a total of four CH$_3$OH emission lines in bands 6 and 7 toward the disk around the young Herbig Ae star HD 163296. The disk-averaged column density of methanol and its related species formaldehyde (H$_2$CO) are estimated assuming optically thin emission in local thermodynamic equilibrium. We compare these results to the gas-phase column densities of the TW Hya disk. No targeted methanol lines were detected individually nor after line stacking. The 3$\sigma$ disk-integrated intensity upper limits are $< 51$ mJy km s$^{-1}$ for the band 6 lines and $< 26$ mJy km s$^{-1}$ for the band 7 lines. The band 7 lines provide the strictest 3$\sigma$ upper limit on disk-averaged column density with $N_{\mathrm{avg}} < 5.0 \times 10^{11}$ cm$^{-2}$. The methanol-to-formaldehyde ratio is CH$_3$OH/H$_2$CO $< 0.24$ in the HD 163296 disk compared to a ratio of $1.27$ in the TW Hya disk. Differences in the stellar irradiation of Herbig disks compared to T Tauri disks likely influence the gaseous methanol and formaldehyde content. Possible reasons for the lower HD 163296 methanol-to-formaldehyde ratio include: a higher than expected gas-phase formation of H$_2$CO in the HD 163296 disk, uncertainties in the grain surface formation efficiency of CH$_3$OH and H$_2$CO, and differences in the disk structure and/or CH$_3$OH and H$_2$CO desorption processes that release the molecules from ice mantles back into the gas phase. These results provide observational evidence that the gas-phase chemical complexity found in disks may be strongly influenced by the spectral type of the host star.

Detailed mapping of the distributions and kinematics of gases in cometary comae at radio wavelengths can provide fundamental advances in our understanding of cometary activity and outgassing mechanisms. Furthermore, the measurement of molecular abundances in comets provides new insights into the chemical composition of some of the Solar System's oldest and most primitive materials. Here we investigate the opportunities for significant progress in cometary science using a very large radio interferometer. The ngVLA concept will enable detection and mapping of a range of key coma species in the 1.2-116 GHz range, and will allow for the first time, high-resolution mapping of the fundamental cometary molecules OH and NH$_3$. The extremely high angular resolution and continuum sensitivity of the proposed ngVLA will also allow the possibility of imaging thermal emission from the nucleus itself, as well as large dust/ice grains in the comae, of comets passing within $\sim1$ au of Earth.

Ortho-to-para (o/p) ratios of species like water, ammonia and formaldehyde (H2CO) are believed to encode information about the formation history of the molecule. Measurements of o/p ratios in protoplanetary disks could thus be used to constrain their physical and chemical histories. We present the first measurement of the H2CO o/p ratio in a protoplanetary disk, using three ortho and two para lines observed with the Sub-millimeter Array (SMA) combined with one highly resolved measurement of a single H2CO line with ALMA toward the disk around Herbig Ae star HD 163296. We find a disk-averaged H2CO o/p ratio of 1.8-2.8 (depending on the assumed disk structure), corresponding to a spin temperature of 11-22 K. We also derive a rotational temperature of 24 K from the flux ratio of the three ortho lines. The observed spatial distribution, as seen by ALMA, as well as the rotational temperature and the o/p ratio, at the large scales the SMA is most sensitive to, are consistent with a low-temperature formation pathway, most likely grain surface chemistry, of H2CO in this disk.

The formation of asteroids, comets and planets occurs in the interior of protoplanetary disks during the early phase of star formation. Consequently, the chemical composition of the disk might shape the properties of the emerging planetary system. In this context, it is crucial to understand whether and what organic molecules are synthesized in the disk. In this Letter, we report the first detection of formic acid (HCOOH) towards the TW Hydrae protoplanetary disk. The observations of the trans-HCOOH 6$_{(1,6)-5(1,5)}$ transition were carried out at 129~GHz with ALMA. We measured a disk-averaged gas-phase t-HCOOH column density of $\sim$ (2-4)$\times$10$^{12}$~cm$^{-2}$, namely as large as that of methanol. HCOOH is the first organic molecules containing two oxygen atoms detected in a protoplanetary disk, a proof that organic chemistry is very active even though difficult to observe in these objects. Specifically, this simplest acid stands as the basis for synthesis of more complex carboxylic acids used by life on Earth.

We present a suite of ALMA interferometric molecular line and continuum images that elucidate, on linear size scales of $\sim$30--40 AU, the chemical structure of the nearby, evolved, protoplanetary disk orbiting the close binary system V4046 Sgr. The observations were undertaken in the 1.1--1.4 mm wavelength range (ALMA Bands 6 and 7) with antenna configurations involving maximum baselines of several hundred meters, yielding subarcsecond-resolution images in more than a dozen molecular species and isotopologues. Isotopologues of CO and HCN display centrally peaked morphologies of integrated emission line intensity, whereas the line emission from complex nitrile group molecules (HC$_3$N, CH$_3$CN), deuterated molecules (DCN, DCO$^+$), hydrocarbons (as traced by C$_2$H), and potential CO ice line tracers (N$_2$H$^+$, and H$_2$CO) appears as a sequence of sharp and diffuse rings. The dimensions and morphologies of HC$_3$N and CH$_3$CN emission are suggestive of photodesorption of organic ices from the surfaces of dust grains, while the sequence of increasing radius of peak intensity represented by DCN (smallest), DCO$^+$, N$_2$H$^+$, and H$_2$CO (largest) is qualitatively consistent with the expected decline of midplane gas temperature with increasing disk radius. Empirical modeling indicates that the sharp-edged C$_2$H emission ring lies at relatively deep disk layers, leaving open the question of the origin of C$_2$H abundance enhancements in evolved disks. This study of the "molecular anatomy" of V4046 Sgr should serve as motivation for additional subarcsecond ALMA molecular line imaging surveys of nearby, evolved protoplanetary disks aimed at addressing major uncertainties in protoplanetary disk physical and chemical structure and molecular production pathways.

We aim to reproduce the DCO$^+$ emission in the disk around HD163296 using a simple 2D chemical model for the formation of DCO$^+$ through the cold deuteration channel and a parametric treatment of the warm deuteration channel. We use data from ALMA in band 6 to obtain a resolved spectral imaging data cube of the DCO$^+$ $J$=3--2 line in HD163296 with a synthesized beam of 0."53$\times$ 0."42. We adopt a physical structure of the disk from the literature that reproduces the spectral energy distribution. We then apply a simplified chemical network for the formation of DCO$^+$ that uses the physical structure of the disk as parameters along with a CO abundance profile, a constant HD abundance and a constant ionization rate. Finally, from the resulting DCO$^+$ abundances, we calculate the non-LTE emission using the 3D radiative transfer code LIME. The observed DCO$^+$ emission is reproduced by a model with cold deuteration producing abundances up to $1.6\times 10^{-11}$. Warm deuteration, at a constant abundance of $3.2\times 10^{-12}$, becomes fully effective below 32 K and tapers off at higher temperatures, reproducing the lack of DCO$^+$ inside 90 AU. Throughout the DCO$^+$ emitting zone a CO abundance of $2\times 10^{-7}$ is found, with $\sim$99\% of it frozen out below 19 K. At radii where both cold and warm deuteration are active, warm deuteration contributes up to 20\% of DCO$^+$, consistent with detailed chemical models. The decrease of DCO$^+$ at large radii is attributed to a temperature inversion at 250 AU, which raises temperatures above values where cold deuteration operates. Increased photodesorption may also limit the radial extent of DCO$^+$. The corresponding return of the DCO$^+$ layer to the midplane, together with a radially increasing ionization fraction, reproduces the local DCO$^+$ emission maximum at $\sim$260 AU.

We present \emph{Herschel} PACS observations of 8 Classical T Tauri Stars in the $\sim 7-10$ Myr old OB1a and the $\sim 4-5$ Myr old OB1b Orion sub-asscociations. Detailed modeling of the broadband spectral energy distributions, particularly the strong silicate emission at 10 $\mu$m, shows that these objects are (pre)transitional disks with some amount of small optically thin dust inside their cavities, ranging from $\sim 4$ AU to $\sim 90$ AU in size. We analyzed \emph{Spitzer} IRS spectra for two objects in the sample: CVSO-107 and CVSO-109. The IRS spectrum of CVSO-107 indicates the presence of crystalline material inside its gap while the silicate feature of CVSO-109 is characterized by a pristine profile produced by amorphous silicates; the mechanisms creating the optically thin dust seem to depend on disk local conditions. Using millimeter photometry we estimated dust disk masses for CVSO-107 and CVSO-109 lower than the minimum mass of solids needed to form the planets in our Solar System, which suggests that giant planet formation should be over in these disks. We speculate that the presence and maintenance of optically thick material in the inner regions of these pre-transitional disks might point to low-mass planet formation.

We present observations with the Submillimeter Array of the continuum emission at $\lambda = 1.3$ mm from 62 young stars surrounded by a protoplanetary disk in the Serpens star-forming region. The typical angular resolution for the survey in terms of beam size is $3.5^{\prime \prime}\times2.5^{\prime \prime}$ with a median rms noise level of 1.6 mJy beam$^{-1}$. These data are used to infer the dust content in disks around low-mass stars $(0.1$-$2.5\,M_{\odot})$ at a median stellar age of $1$-$3$ Myr. Thirteen sources were detected in the 1.3 mm dust continuum with inferred dust masses of ${\approx} 10$-$260\,M_{\oplus}$ and an upper limit to the median dust mass of $5.1_{-4.3}^{+6.1}\,M_{\oplus}$, derived using survival analysis. Comparing the protoplanetary disk population in Serpens to those of other nearby star-forming regions, we find that the populations of dust disks in Serpens and Taurus, which have a similar age, are statistically indistinguishable. This is potentially surprising since Serpens has a stellar surface density two orders of magnitude in excess of Taurus. Hence, we find no evidence that dust disks in Serpens have been dispersed as a result of more frequent and/or stronger tidal interactions due its elevated stellar density. We also report that the fraction of Serpens disks with $M_{\rm{dust}} \geq 10\,M_{\oplus}$ is less than 20%, which supports the notion that the formation of giant planets is likely inherently rare or has substantially progressed by a few Myrs.

The formation pathways of deuterated species trace different regions of protoplanetary disks and may shed light into their physical structure. We aim to constrain the radial extent of main deuterated species; we are particularly interested in spatially characterizing the high and low temperature pathways for enhancing deuteration of these species. We observed the disk surrounding the Herbig Ae star HD 163296 using ALMA in Band 6 and obtained resolved spectral imaging data of DCO$^+$ ($J$=3-2), DCN ($J$=3-2) and N$_2$D$^+$ ($J$=3-2). We model the radial emission profiles of DCO$^+$, DCN and N$_2$D$^+$, assuming their emission is optically thin, using a parametric model of their abundances and radial excitation temperature estimates. DCO$^+$ can be described by a three-region model, with constant-abundance rings centered at 70 AU, 150 AU and 260 AU. The DCN radial profile peaks at about ~60 AU and N$_2$D$^+$ is seen in a ring at ~160 AU. Simple models of both molecules using constant abundances reproduce the data. Assuming reasonable average excitation temperatures for the whole disk, their disk-averaged column densities (and deuterium fractionation ratios) are 1.6-2.6$\times 10^{12}$ cm$^{-2}$ (0.04-0.07), 2.9-5.2$\times 10^{12}$ cm$^{-2}$ ($\sim$0.02) and 1.6-2.5 $\times 10^{11}$ cm$^{-2}$ (0.34-0.45) for DCO$^+$, DCN and N$_2$D$^+$, respectively. Our simple best-fit models show a correlation between the radial location of the first two rings in DCO$^+$ and the DCN and N$_2$D$^+$ abundance distributions that can be interpreted as the high and low temperature deuteration pathways regimes. The origin of the third DCO$^+$ ring at 260 AU is unknown but may be due to a local decrease of ultraviolet opacity allowing the photodesorption of CO or due to thermal desorption of CO as a consequence of radial drift and settlement of dust grains.

Isotopic ratios are keys to understanding the origin and early evolution of the solar system in the context of Galactic nucleosynthesis. The large range of measured $^{14}$N/$^{15}$N isotopic ratios in the solar system reflects distinct reservoirs of nitrogen whose origins remain to be determined. We have directly measured a C$^{14}$N/C$^{15}$N abundance ratio of 323$\pm$30 in the disk orbiting the nearby young star TW Hya. This value, which is in good agreement with nitrogen isotopic ratios measured for prestellar cores, likely reflects the primary present-day reservoir of nitrogen in the solar neighbourhood. These results support models invoking novae as primary $^{15}$N sources as well as outward migration of the Sun over its lifetime, and suggest that comets sampled a secondary, $^{15}$N-rich reservoir during solar system formation.

Gas kinematics are an important part of the planet formation process. Turbulence influences planetesimal growth and migration from the scale of sub-micron dust grains through gas-giant planets. Radio observations of resolved molecular line emission can directly measure this non-thermal motion and, taking advantage of the layered chemical structure of disks, different molecular lines can be combined to map the turbulence throughout the vertical extent of a protoplanetary disk. Here we present ALMA observations of three molecules (DCO$^+$(3-2), C$^{18}$O(2-1) and CO(2-1)) from the disk around HD 163296. We are able to place stringent upper limits ($v_{\rm turb}<$0.06c$_s$, $<$0.05c$_s$ and $<$0.04c$_s$ for CO(2-1), C$^{18}$O(2-1) and DCO$^+$(3-2) respectively), corresponding to $\alpha\lesssim$3$\times$10$^{-3}$, similar to our prior limit derived from CO(3-2). This indicates that there is little turbulence throughout the vertical extent of the disk, contrary to theoretical predictions based on the magneto-rotational instability and gravito-turbulence. In modeling the DCO$^+$ emission we also find that it is confined to three concentric rings at 65.7$\pm$0.9 au, 149.9$^{+0.5}_{-0.7}$ au and 259$\pm$1 au, indicative of a complex chemical environment.

Three formaldehyde lines were observed (H$_2$CO 3$_{03}$--2$_{02}$, H$_2$CO 3$_{22}$--2$_{21}$, and H$_2$CO 3$_{21}$--2$_{20}$) in the protoplanetary disk around the Herbig Ae star HD 163296 with ALMA at 0.5 arcsecond (60 AU) spatial resolution. H$_2$CO 3$_{03}$--2$_{02}$ was readily detected via imaging, while the weaker H$_2$CO 3$_{22}$--2$_{21}$ and H$_2$CO 3$_{21}$--2$_{20}$ lines required matched filter analysis to detect. H$_2$CO is present throughout most of the gaseous disk, extending out to 550 AU. An apparent 50 AU inner radius of the H$_2$CO emission is likely caused by an optically thick dust continuum. The H$_2$CO radial intensity profile shows a peak at 100 AU and a secondary bump at around 300 AU, suggesting increased production in the outer disk. Different parameterizations of the H$_2$CO abundance were compared to the observed visibilities with $\chi^2$ minimization, using either a characteristic temperature, a characteristic radius or a radial power law index to describe the H$_2$CO chemistry. Similar models were applied to ALMA Science Verification data of C$^{18}$O. In all modeling scenarios, fits to the H$_2$CO data show an increased abundance in the outer disk. The overall best-fit H$_2$CO model shows a factor of two enhancement beyond a radius of 270$\pm$20 AU, with an inner abundance of $2\!-\!5 \times 10^{-12}$. The H$_2$CO emitting region has a lower limit on the kinetic temperature of $T > 20$ K. The C$^{18}$O modeling suggests an order of magnitude depletion in the outer disk and an abundance of $4\!-\!12 \times 10^{-8}$ in the inner disk. The increase in H$_2$CO outer disk emission could be a result of hydrogenation of CO ices on dust grains that are then sublimated via thermal desorption or UV photodesorption, or more efficient gas-phase production beyond about 300 AU if CO is photodisocciated in this region.

H2CO is one of the most readily detected organic molecules in protoplanetary disks. Yet its distribution and dominant formation pathway(s) remain largely unconstrained. To address these issues, we present ALMA observations of two H2CO lines (3_{12}-2_{11} and 5_{15}-4_{14}) at 0".5 (~30 au) spatial resolution toward the disk around the nearby T Tauri star TW Hya. Emission from both lines is spatially resolved, showing a central depression, a peak at 0".4 radius, and a radial decline at larger radii with a bump at ~1", near the millimeter continuum edge. We adopt a physical model for the disk and use toy models to explore the radial and vertical H2CO abundance structure. We find that the observed emission implies the presence of at least two distinct H2CO gas reservoirs: (1) a warm and unresolved inner component (<10 au), and (2) an outer component that extends from ~15 au to beyond the millimeter continuum edge. The outer component is further constrained by the line ratio to arise in a more elevated disk layer at larger radii. The inferred H2CO abundance structure agrees well with disk chemistry models, which predict efficient H2CO gas-phase formation close to the star, and cold H2CO grain surface formation, through H additions to condensed CO, followed by non-thermal desorption in the outer disk. The implied presence of active grain surface chemistry in the TW Hya disk is consistent with the recent detection of CH3OH emission, and suggests that more complex organic molecules are formed in disks, as well.

The deuterium enrichment of molecules is sensitive to their formation environment. Constraining patterns of deuterium chemistry in protoplanetary disks is therefore useful for probing how material is inherited or reprocessed throughout the stages of star and planet formation. We present ALMA observations at $\sim0.6"$ resolution of DCO$^+$, H$^{13}$CO$^+$, DCN, and H$^{13}$CN in the full disks around T Tauri stars AS 209 and IM Lup, the transition disks around T Tauri stars V4046 Sgr and LkCa 15, and the full disks around Herbig Ae stars MWC 480 and HD 163296. We also present ALMA observations of HCN in the IM Lup disk. DCN, DCO$^+$, and H$^{13}$CO$^+$ are detected in all disks, and H$^{13}$CN in all but the IM Lup disk. We find efficient deuterium fractionation for the sample, with estimates of disk-averaged DCO$^+$/HCO$^+$ and DCN/HCN abundance ratios ranging from $\sim0.02-0.06$ and $\sim0.005-0.08$, respectively, which is comparable to values reported for other ISM environments. The relative distributions of DCN and DCO$^+$ vary between disks, suggesting that multiple formation pathways may be needed to explain the diverse emission morphologies. In addition, gaps and rings observed in both H$^{13}$CO$^+$ and DCO$^+$ emission provide new evidence that DCO$^+$ bears a complex relationship with the location of the midplane CO snowline.

Nitrogen fractionation is commonly used to assess the thermal history of Solar System volatiles. With ALMA it is for the first time possible to directly measure 14N/15N ratios in common molecules during the assembly of planetary systems. We present ALMA observations of the H13CN and HC15N J=3-2 lines at 0".5 angular resolution, toward a sample of six protoplanetary disks, selected to span a range of stellar and disk structure properties. Adopting a typical 12C/13C ratio of 70, we find comet-like 14N/15N ratios of 80-160 in 5/6 of the disks (3 T Tauri and 2 Herbig Ae disks) and lack constraints for one of the T Tauri disks (IM Lup). There are no systematic differences between T Tauri and Herbig Ae disks, or between full and transition disks within the sample. In addition, no correlation is observed between disk-averaged D/H and 14N/15N ratios in the sample. One of the disks, V4046 Sgr, presents unusually bright HCN isotopologue emission, enabling us to model the radial profiles of H13CN and HC15N. We find tentative evidence of an increasing 14N/15N ratio with radius, indicating that selective photodissociation in the inner disk is important in setting the 14N/15N ratio during planet formation.

CO is widely used as a tracer of molecular gas. However, there is now mounting evidence that gas phase carbon is depleted in the disk around TW Hya. Previous efforts to quantify this depletion have been hampered by uncertainties regarding the radial thermal structure in the disk. Here we present resolved ALMA observations of 13CO 3-2, C18O 3-2, 13CO 6-5, and C18O 6-5 emission in TW Hya, which allow us to derive radial gas temperature and gas surface density profiles, as well as map the CO abundance as a function of radius. These observations provide a measurement of the surface CO snowline at ~30 AU and show evidence for an outer ring of CO emission centered at 53 AU, a feature previously seen only in less abundant species. Further, the derived CO gas temperature profile constrains the freeze-out temperature of CO in the warm molecular layer to < 21 K. Combined with the previous detection of HD 1-0, these data constrain the surface density of the warm H2 gas in the inner ~30 AU. We find that CO is depleted by two orders of magnitude from R=10-60 AU, with the small amount of CO returning to the gas phase inside the surface CO snowline insufficient to explain the overall depletion. Finally, this new data is used in conjunction with previous modeling of the TW Hya disk to constrain the midplane CO snowline to 17-23 AU.

The disk around HD 141569 is one of a handful of systems whose weak infrared emission is consistent with a debris disk, but still has a significant reservoir of gas. Here we report spatially resolved mm observations of the CO(3-2) and CO(1-0) emission as seen with the SMA and CARMA. We find that the excitation temperature for CO is lower than expected from cospatial blackbody grains, similar to previous observations of analogous systems, and derive a gas mass that lies between that of gas-rich primordial disks and gas-poor debris disks. The data also indicate a large inner hole in the CO gas distribution and an outer radius that lies interior to the outer scattered light rings. This spatial distribution, with the dust rings just outside the gaseous disk, is consistent with the expected interactions between gas and dust in an optically thin disk. This indicates that gas can have a significant effect on the location of the dust within debris disks.

Grain growth in planet-forming disks is the first step toward the formation of planets. The growth of grains and their inward drift leaves a distinct imprint on the dust surface-density distribution and the resulting surface-brightness profile of the thermal continuum emission. We determine the surface-brightness profile of the continuum emission using resolved observations at millimeter wavelengths of the disk around TW Hya, and infer the signature of dust evolution on the surface density and dust opacity. Archival ALMA observations at 820 micron on baselines up to 410 kilolambda are compared to parametrized disk models to determine the surface-brightness profile. Under the assumption of a constant dust opacity, a broken radial power law best describes the dust surface density, with a slope of -0.53 +/- 0.01 from the 4.1 au radius of the (already known) inner hole to a turn-over radius of 47.1 +/- 0.2 au, steepening to -8.0 +/- 0.1 at larger radii. The emission drops below the detection limit beyond ~60 au. The shape of the dust surface density is consistent with theoretical expectations for grain growth, fragmentation, and drift, but its total dust content and its turn-over radius are too large for TW Hya's age of 8-10 Myr even when taking into account a radially varying dust opacity. Higher resolution imaging with ALMA of TW Hya and other disks is required to establish if unseen gaps associated with, e.g., embedded planets trap grains at large radii or if locally enhanced grain growth associated with the CO snow line explains the extent of the millimeter-continuum surface brightness profile. In the latter case, population studies should reveal a correlation between the location of the CO snow line and the extent of the millimeter continuum. In the former case, and if CO freeze out promotes planet formation, this correlation should extend to the location of gaps as well.

HCN is a commonly observed molecule in Solar System bodies and in interstellar environments. Its abundance with respect to CN is a proposed tracer of UV exposure. HCN is also frequently used to probe the thermal history of objects, by measuring its degree of nitrogen fractionation. To address the utility of HCN as a probe of disks, we present ALMA observations of CN, HCN, H13CN and HC15N toward the protoplanetary disk around Herbig Ae star MWC480, and of CN and HCN toward the disk around T Tauri star DM Tau. Emission from all molecules is clearly detected and spatially resolved, including the first detection of HC15N in a disk. Toward MWC 480, CN emission extends radially more than 1" exterior to the observed cut-off of HCN emission. Quantitative modeling further reveals very different radial abundance profiles for CN and HCN, with best-fit outer cut-off radii of >300 AU and 110+-10 AU, respectively. This result is in agreement with model predictions of efficient HCN photodissociation into CN in the outer-part of the disk where the vertical gas and dust column densities are low. No such difference in CN and HCN emission profiles are observed toward DM Tau, suggestive of different photochemical structures in Herbig Ae and T Tauri disks. We use the HCN isotopologue data toward the MWC 480 disk to provide the first measurement of the 14N/15N ratio in a disk. We find a low disk averaged 14N/15N ratio of 200+-100, comparable to what is observed in cloud cores and comets, demonstrating interstellar inheritance and/or efficient nitrogen fractionation in this disk.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

The condensation fronts (snow lines) of H2O, CO and other abundant volatiles in the midplane of a protoplanetary disk affect several aspects of planet formation. Locating the CO snow line, where the CO gas column density is expected to drop substantially, based solely on CO emission profiles is challenging. This has prompted an exploration of chemical signatures of CO freeze-out. We present ALMA Cycle 1 observations of the N2H+ J=3-2 and DCO+ J=4-3 emission lines toward the disk around the Herbig Ae star HD~163296 at ~0.5" (60 AU) resolution, and evaluate their utility as tracers of the CO snow line location. The N2H+ emission is distributed in a ring with an inner radius at 90 AU, corresponding to a midplane temperature of 25 K. This result is consistent with a new analysis of optically thin C18O data, which implies a sharp drop in CO abundance at 90 AU. Thus N2H+ appears to be a robust tracer of the midplane CO snow line. The DCO+ emission also has a ring morphology, but neither the inner nor the outer radius coincides with the CO snow line location of 90 AU, indicative of a complex relationship between DCO+ emission and CO freeze-out in the disk midplane. Compared to TW Hya, CO freezes out at a higher temperature in the disk around HD 163296 (25 vs. 17 K in the TW Hya disk), perhaps due to different ice compositions. This highlights the importance of actually measuring the CO snow line location, rather than assuming a constant CO freeze-out temperature for all disks.

We have used the Submillimeter Array to image, at ~1" resolution, C2H(3-2) emission from the molecule-rich circumstellar disks orbiting the nearby, classical T Tauri star systems TW Hya and V4046 Sgr. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology within each disk, the inner hole radius of the C2H ring within the V4046 Sgr disk (~70 AU) is somewhat larger than than of its counterpart within the TW Hya disk (~45 AU). We suggest that, in each case, the C2H emission likely traces irradiation of the tenuous surface layers of the outer disks by high-energy photons from the central stars.

In a protoplanetary disk, a combination of thermal and non-thermal desorption processes regulate where volatiles are liberated from icy grain mantles into the gas phase. Non-thermal desorption should result in volatile-enriched gas in disk-regions where complete freeze-out is otherwise expected. We present ALMA observations of the disk around the young star IM Lup in 1.4 mm continuum, C18O 2-1, H13CO+ 3-2 and DCO+ 3-2 emission at ~0".5 resolution. The images of these dust and gas tracers are clearly resolved. The DCO+ line exhibits a striking pair of concentric rings of emission that peak at radii of ~0".6 and 2" (~90 and 300 AU, respectively). Based on disk chemistry model comparison, the inner DCO+ ring is associated with the balance of CO freeze-out and thermal desorption due to a radial decrease in disk temperature. The outer DCO+ ring is explained by non-thermal desorption of CO ice in the low-column-density outer disk, repopulating the disk midplane with cold CO gas. The CO gas then reacts with abundant H2D+ to form the observed DCO+ outer ring. These observations demonstrate that spatially resolved DCO+ emission can be used to trace otherwise hidden cold gas reservoirs in the outmost disk regions, opening a new window onto their chemistry and kinematics.

We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.10" (13 AU) out to 0.50" (65 AU). This inner working angle is interior to the location of the system's gap inferred by previous studies using SED modeling (15 AU). We detected a candidate point source companion 1.08" northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50 deg), 0.01 Msun disk that has a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to variety of mechanisms.

The distributions and abundances of small organics in protoplanetary disks are potentially powerful probes of disk physics and chemistry. HNC is a common probe of dense interstellar regions and the target of this study. We use the Submillimeter Array (SMA) to observe HNC 3--2 towards the protoplanetary disks around the T Tauri star TW Hya and the Herbig Ae star HD 163296. HNC is detected toward both disks, constituting the first spatially resolved observations of HNC in disks. We also present SMA observations of HCN 3--2, and IRAM 30m observations of HCN and HNC 1--0 toward HD 163296. The disk-averaged HNC/HCN emission ratio is 0.1--0.2 toward both disks. Toward TW Hya, the HNC emission is confined to a ring. The varying HNC abundance in the TW Hya disk demonstrates that HNC chemistry is strongly linked to the disk physical structure. In particular, the inner rim of the HNC ring can be explained by efficient destruction of HNC at elevated temperatures, similar to what is observed in the ISM. To realize the full potential of HNC as a disk tracer requires, however, a combination of high SNR spatially resolved observations of HNC and HCN, and disk specific HNC chemical modeling.

Resolved submillimeter imaging of transitional disks is increasingly revealing the complexity of disk structure. Here we present the first high-resolution submillimeter image of a recently identified transitional disk around IRAS 04125+2902 in the Taurus star-forming region. We measure an inner disk hole of ~20 AU around IRAS 04125+2902 by simultaneously modeling new 880 micron Submillimeter Array (SMA) data along with an existing spectral energy distribution supplemented by new Discovery Channel Telescope (DCT) photometry. We also constrain the outer radius of the dust disk in IRAS~04125+2902 to ~50-60 AU. Such a small dust disk could be attributed to initial formation conditions, outward truncation by an unseen companion, or dust evolution in the disk. Notably, the dust distribution of IRAS 04125+2902 resembles a narrow ring (delta R ~ 35 AU) composed of large dust grains at the location of the disk wall. Such narrow dust rings are also seen in other transitional disks and may be evidence of dust trapping in pressure bumps, possibly produced by planetary companions. More sensitive submillimeter observations of the gas are necessary to further probe the physical mechanisms at work in shaping the spatial distribution of large dust in this disk. Interestingly, the IRAS 04125+2902 disk is significantly fainter than other transitional disks that have been resolved at submillimeter wavelengths, hinting that more objects with large disk holes may exist at the faint end of the submillimeter luminosity distribution that await detection with more sensitive imaging telescopes.

We investigate the chemistry of ion molecules in protoplanetary disks, motivated by the detection of N$_2$H$^+$ ring around TW Hya. While the ring inner radius coincides with the CO snow line, it is not apparent why N$_2$H$^+$ is abundant outside the CO snow line in spite of the similar sublimation temperatures of CO and N$_2$. Using the full gas-grain network model, we reproduced the N$_2$H$^+$ ring in a disk model with millimeter grains. The chemical conversion of CO and N$_2$ to less volatile species (sink effect hereinafter) is found to affect the N$_2$H$^+$ distribution. Since the efficiency of the sink depends on various parameters such as activation barriers of grain surface reactions, which are not well constrained, we also constructed the no-sink model; the total (gas and ice) CO and N$_2$ abundances are set constant, and their gaseous abundances are given by the balance between adsorption and desorption. Abundances of molecular ions in the no-sink model are calculated by analytical formulas, which are derived by analyzing the full-network model. The N$_2$H$^+$ ring is reproduced by the no-sink model, as well. The 2D (R-Z) distribution of N$_2$H$^+$, however, is different among the full-network model and no-sink model. The column density of N$_2$H$^+$ in the no-sink model depends sensitively on the desorption rate of CO and N$_2$, and the flux of cosmic ray. We also found that N$_2$H$^+$ abundance can peak at the temperature slightly below the CO sublimation, even if the desorption energies of CO and N$_2$ are the same.

Observations of comets and asteroids show that the Solar Nebula that spawned our planetary system was rich in water and organic molecules. Bombardment brought these organics to the young Earth's surface, seeding its early chemistry. Unlike asteroids, comets preserve a nearly pristine record of the Solar Nebula composition. The presence of cyanides in comets, including 0.01% of methyl cyanide (CH3CN) with respect to water, is of special interest because of the importance of C-N bonds for abiotic amino acid synthesis. Comet-like compositions of simple and complex volatiles are found in protostars, and can be readily explained by a combination of gas-phase chemistry to form e.g. HCN and an active ice-phase chemistry on grain surfaces that advances complexity[3]. Simple volatiles, including water and HCN, have been detected previously in Solar Nebula analogues - protoplanetary disks around young stars - indicating that they survive disk formation or are reformed in situ. It has been hitherto unclear whether the same holds for more complex organic molecules outside of the Solar Nebula, since recent observations show a dramatic change in the chemistry at the boundary between nascent envelopes and young disks due to accretion shocks[8]. Here we report the detection of CH3CN (and HCN and HC3N) in the protoplanetary disk around the young star MWC 480. We find abundance ratios of these N-bearing organics in the gas-phase similar to comets, which suggests an even higher relative abundance of complex cyanides in the disk ice. This implies that complex organics accompany simpler volatiles in protoplanetary disks, and that the rich organic chemistry of the Solar Nebula was not unique.

We have used the Submillimeter Array to image, at ~1.5" resolution, C2H (3-2) emission from the circumstellar disk orbiting the nearby (D = 54 pc), ~8 Myr-old, ~0.8 Msun classical T Tauri star TW Hya. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology. Based on a model in which the C2H column density follows a truncated radial power-law distribution, we find that the inner edge of the ring lies at ~45 AU, and that the ring extends to at least ~120 AU. Comparison with previous (single-dish) observations of C2H (4-3) emission indicates that the C2H molecules are subthermally excited and, hence, that the emission arises from the relatively warm, tenuous upper atmosphere of the disk. We propose that the C2H emission most likely traces particularly efficient photo-destruction of small grains and/or photodesorption and photodissociation of hydrocarbons derived from grain ice mantles in the surface layers of the outer disk. The presence of a C2H ring in the TW Hya disk hence likely serves as a marker of dust grain processing and radial and vertical grain size segregation within the disk.

In this Letter we model the chemistry of DCO$^{+}$ in protoplanetary disks. We find that the overall distribution of the DCO$^{+}$ abundance is qualitatively similar to that of CO but is dominated by thin layer located at the inner disk surface. To understand its distribution, we investigate the different key gas-phase deuteration pathways that can lead to the formation of DCO$^{+}$. Our analysis shows that the recent update in the exothermicity of the reaction involving CH$_2$D$^{+}$ as a parent molecule of DCO$^{+}$ favors deuterium fractionation in warmer conditions. As a result the formation of DCO$^{+}$ is enhanced in the inner warm surface layers of the disk where X-ray ionization occurs. Our analysis points out that DCO$^{+}$ is not a reliable tracer of the CO snow line as previously suggested. We thus predict that DCO$^{+}$ is a tracer of active deuterium and in particular X-ray ionization of the inner disk.

We present an observational and theoretical study of the primary ionizing agents (cosmic rays and X-rays) in the TW Hya protoplanetary disk. We use a set of resolved and unresolved observations of molecular ions and other molecular species, encompassing eleven lines total, in concert with a grid of disk chemistry models. The molecular ion constraints comprise new data from the Submillimeter Array on HCO$^+$, acquired at unprecedented spatial resolution, and data from the literature, including ALMA observations of N$_2$H$^+$. We vary the model incident CR flux and stellar X-ray spectra and find that TW Hya's HCO$^+$ and N$_2$H$^+$ emission are best fit by a moderately hard X-ray spectra, as would be expected during the "flaring" state of the star, and a low CR ionization rate, $\zeta_{\rm CR}\lesssim10^{-19}$ s$^{-1}$. This low CR rate is the first indication of the presence of CR exclusion by winds and/or magnetic fields in an actively accreting T Tauri disk system. With this new constraint, our best fit ionization structure predicts a low turbulence "dead-zone" extending from the inner edge of the disk out to $50-65$ AU. This region coincides with an observed concentration of millimeter grains, and we propose that the inner region of TW Hya is a dust (and possibly planet) growth factory as predicted by previous theoretical work.

We present high angular resolution Submillimeter Array observations ofthe outbursting Jupiter family comet 17P/Holmes on 2007 October 26-29, achieving a spatial resolution of 2.5", or ~3000 km at the comet distance. The observations resulted in detections of the rotational lines CO 3-2, HCN 4-3, H$^{13}$CN 4-3, CS 7-6, H$_2$CO 3$_{1,2}$-2$_{1,1}$, H$_2$S 2$_{2,0}$-2$_{1,1}$, and multiple CH$_3$OH lines, along with the associated dust continuum at 221 and 349 GHz. The continuum has a spectral index of 2.7$\pm$0.3, slightly steeper than blackbody emission from large dust particles. From the imaging data, we identify two components in the molecular emission. One component is characterized by a relatively broad line width (~1 km s$^{-1}$ FWHM) exhibiting a symmetric outgassing pattern with respect to the nucleus position. The second component has a narrower line width (<0.5 km s$^{-1}$ FWHM) with the line center red-shifted by 0.1-0.2 km s$^{-1}$ (cometocentric frame), and shows a velocity shift across the nucleus position with the position angle gradually changing from 66 to 30 degrees within the four days of observations. We determine distinctly different CO/HCN ratios for each of the components. For the broad-line component we find CO/HCN <7, while in the narrow-line component, CO/HCN = 40$\pm$5. We hypothesize that the narrow-line component originates from the ice grain halo found in near-nucleus photometry, believed to be created by sublimating recently released ice grains around the nucleus during the outburst. In this interpretation, the high CO/HCN ratio of this component reflects the more pristine volatile composition of nucleus material released in the outburst.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

In this Letter we report the CO abundance relative to H2 derived toward the circumstellar disk of the T-Tauri star TW Hya from the HD (1-0) and C18O (2-1) emission lines. The HD (1-0) line was observed by the Herschel Space Observatory Photodetector Array Camera and Spectrometer whereas C18O (2-1) observations were carried out with the Submillimeter Array at a spatial resolution of 2.8" x 1.9" (corresponding to 142 x 97 AU). In the disk's warm molecular layer (T>20 K) we measure a disk-averaged gas-phase CO abundance relative to H2 of $\chi{\rm(CO)}=(0.1-3)x10^{-5}$, substantially lower than the canonical value of $\chi{\rm(CO)}=10^{-4}$. We infer that the best explanation of this low $\chi$(CO) is the chemical destruction of CO followed by rapid formation of carbon chains, or perhaps CO2, that can subsequently freeze-out, resulting in the bulk mass of carbon locked up in ice grain mantles and oxygen in water. As a consequence of this likely time-dependent carbon sink mechanism, CO may be an unreliable tracer of H2 gas mass.

Planets form in the disks around young stars. Their formation efficiency and composition are intimately linked to the protoplanetary disk locations of "snow lines" of abundant volatiles. We present chemical imaging of the CO snow line in the disk around TW Hya, an analog of the solar nebula, using high spatial and spectral resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of N2H+, a reactive ion present in large abundance only where CO is frozen out. The N2H+ emission is distributed in a large ring, with an inner radius that matches CO snow line model predictions. The extracted CO snow line radius of ~ 30 AU helps to assess models of the formation dynamics of the Solar System, when combined with measurements of the bulk composition of planets and comets.

We analyze sensitive, sub-arcsecond resolution ALMA Science Verification observations of CO emission lines in the protoplanetary disk hosted by the young, isolated Ae star HD 163296. The observed spatial morphology of the 12CO J=3-2 emission line is asymmetric across the major axis of the disk; the 12CO J=2-1 line features a much less pronounced, but similar, asymmetry. The J=2-1 emission from 12CO and its main isotopologues have no resolved spatial asymmetry. We associate this behavior as the direct signature of a vertical temperature gradient and layered molecular structure in the disk. This is demonstrated using both toy models and more sophisticated calculations assuming non-local thermodynamic equilibrium (non LTE) conditions. A model disk structure is developed to reproduce both the distinctive spatial morphology of the 12CO J=3-2 line as well as the J=2-1 emission from the CO isotopologues assuming relative abundances consistent with the interstellar medium. This model disk structure has tau=1 emitting surfaces for the 12CO emission lines that make an angle of about 15 degrees with respect to the disk midplane. Furthermore, we show that the spatial and spectral sensitivity of these data can distinguish between models that have sub-Keplerian gas velocities due to the vertical extent of the disk and its associated radial pressure gradient (a fractional difference in the bulk gas velocity field of approximately greater than 5%).

From the masses of planets orbiting our Sun, and relative elemental abundances, it is estimated that at birth our Solar System required a minimum disk mass of ~0.01 solar masses within ~100 AU of the star. The main constituent, gaseous molecular hydrogen, does not emit from the disk mass reservoir, so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally probes the disk surface, while the conversion from dust emission to gas mass requires knowledge of the grain properties and gas-to-dust mass ratio, which likely differ from their interstellar values. Thus, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3--10 Myr) star TW Hya, with estimates ranging from 0.0005 to 0.06 solar masses. Here we report the detection the fundamental rotational transition of hydrogen deuteride, HD, toward TW Hya. HD is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The HD detection, combined with existing observations and detailed models, implies a disk mass >0.05 solar masses, enough to form a planetary system like our own.

We report the first detection of c-C3H2 in a circumstellar disk. The c-C3H2 J=6-5 line (217.882 GHz) is detected and imaged through Atacama Large Millimeter Array (ALMA) Science Verification observations toward the disk around the Herbig Ae star HD 163296 at 0.8" resolution. The emission is consistent with that arising from a Keplerian rotating disk. Two additional c-C3H2 transitions are also tentatively detected, bolstering the identification of this species, but with insufficient signal-to-noise ratio to constrain the spatial distribution. Using a previously developed model for the physical structure of this disk, we fit a radial power-law distribution model to the c-C3H2 6-5 emission and find that c-C3H2 is present in a ring structure from an inner radius of about 30 AU to an outer radius of about 165 AU. The column density is estimated to be 1e12-1e13 cm-2. The clear detection and intriguing ring structure suggest that c-C3H2 has the potential to become a useful probe of radiation penetration in disks.

We present Submillimeter Array observations of H2CO and N2H+ emission in the disks around the T Tauri star TW Hya and the Herbig Ae star HD 163296 at 2"-6" resolution and discuss the distribution of these species with respect to CO freeze-out. The H2CO and N2H+ emission toward HD 163296 does not peak at the continuum emission center that marks the stellar position but is instead significantly offset. Using a previously developed model for the physical structure of this disk, we show that the H2CO observations are reproduced if H2CO is present predominantly in the cold outer disk regions. A model where H2CO is present only beyond the CO snow line (estimated at a radius of 160 AU) matches the observations well. We also show that the average H2CO excitation temperature, calculated from two transitions of H2CO observed in these two disks and a larger sample of disks around T Tauri stars in the DISCS (the Disk Imaging Survey of Chemistry with SMA) program, is consistent with the CO freeze-out temperature of 20 K. In addition, we show that N2H+ and H2CO line fluxes in disks are strongly correlated, indicative of co-formation of these species across the sample. Taken together, these results imply that H2CO and N2H+ are generally present in disks only at low temperatures where CO depletes onto grains, consistent with fast destruction of N2H+ by gas-phase CO, and in situ formation of H2CO through hydrogenation of CO ice. In this scenario H2CO, CH3OH and N2H+ emission in disks should appear as rings with the inner edge at the CO midplane snow line. This prediction can be tested directly using observations from ALMA with higher resolution and better sensitivity.

We present observations of the debris disk around gamma Doradus, an F1V star, from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved at 70, 100 and 160 micron, resolved along its major axis at 250 micron, detected but not resolved at 350 micron, and confused with a background source at 500 micron. It is one of our best resolved targets and we find it to have a radially broad dust distribution. The modelling of the resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring at several AU (best-fit 4 AU) and a cool outer belt extending from ~55 to 400 AU or an arrangement of two cool, narrow rings at ~70 AU and ~190 AU. This suggests that any configuration between these two is also possible. Both models have a total fractional luminosity of ~10^{-5} and are consistent with the disk being aligned with the stellar equator. The inner edge of either possible configuration suggests that the most likely region to find planets in this system would be within ~55 AU of the star. A transient event is not needed to explain the warm dust's fractional luminosity.

We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.

The distributions of deuterated molecules in protoplanetary disks are expected to depend on the molecular formation pathways. We use observations of spatially resolved DCN emission from the disk around TW Hya, acquired during ALMA Science verification with a ~3" synthesized beam, together with comparable DCO+ observations from the Submillimeter Array, to investigate differences in the radial distributions of these species and hence differences in their formation chemistry. In contrast to DCO+, which shows an increasing column density with radius, DCN is better fit by a model that is centrally peaked. We infer that DCN forms at a smaller radii and thus at higher temperatures than DCO+. This is consistent with chemical network model predictions of DCO+ formation from H2D+ at T<30 K and DCN formation from additional pathways involving CH2D+ at higher temperatures. We estimate a DCN/HCN abundance ratio of ~0.017, similar to the DCO+/HCO+ abundance ratio. Deuterium fractionation appears to be efficient at a range of temperatures in this protoplanetary disk. These results suggest caution in interpreting the range of deuterium fractions observed in Solar System bodies, as multiple formation pathways should be taken into account.

Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the disk and we call these objects "transitional" and pre-transitional" disks, respectively. For the latter group, we can fit these SEDs with full disk models and find that millimeter data are necessary to break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only be applied to objects that display evidence for a radical change in the disk's radial structure. Using this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. These reduced accretion rates onto the star could be linked to forming planets. Future observations of transitional and pre-transitional disks will allow us to better quantify the signatures of planet formation in young disks.

We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

We present millimeter-wave observations of several molecular ions in the disk around the pre-main-sequence star DM Tau and use these to investigate the ionization fraction in different regions of the disk. New Submillimeter Array (SMA) observations of H2D+ J=1_10 - 1_11, N2H+ J=4-3 and CO J=3-2 are presented. H2D+ and N2H+ are not detected and using the CO 3-2 disk size the observations result in an upper limit of <0.47 K km s-1 for both lines, a factor of 2.5 below previous single-dish H2D+ observations. Assuming LTE, a disk midplane temperature of 10-20 K and estimates of the H2D+ o/p ratio, the observed limit corresponds to NH2D+ < 4 - 21 \times 1012 cm-2. We adopt a parametric model for the disk structure from the literature and use new IRAM 30 meter telescope observations of the H13CO+ J=3-2 line and previously published SMA observations of the N2H+ J=3-2, HCO+ J=3-2 and DCO+ J=3-2 lines to constrain the ionization fraction, xi, in three temperature regions in the disk where theoretical considerations suggest different ions should dominate: (1) a warm, upper layer with T>20 K where CO is in the gas-phase and HCO+ is most abundant, where we estimate xi \simeq 4 \times 10-10, (2) a cooler molecular layer with T = 16-20 K where N2H+ and DCO+ abundances are predicted to peak, with xi \simeq 3\times10-11, and (3) the cold, dense midplane with T<16 K where H3+ and its deuterated isotopologues are the main carriers of positive charge, with xi < 3\times10-10. While there are considerable uncertainties, these estimates are consistent with a decreasing ionization fraction into the deeper, colder, and denser disk layers. Stronger constraints on the ionization fraction in the disk midplane will require not only substantially more sensitive observations of the H2D+ 1_10 - 1_11 line, but also robust determinations of the o/p ratio, observations of D2H+ and stronger constraints on where N2 is present in the gas phase.

We present \sim1" resolution 1.3 millimeter dust continuum and spectral line (12CO and 13CO J=2-1) observations of the transitional disk system HD 135344B obtained with the Submillimeter Array. The disk shows a Keplerian rotation pattern with an inclination of \sim11{\deg}, based on the spatially and spectrally resolved 12CO and 13CO emission. The data show clear evidence for both dust and gas surface density reductions in the inner region of the disk (radius \lesssim 50 AU) from the continuum and 13CO J=2-1 data, respectively. The presence of this inner cavity in both the dust and gas is more consistent with clearing by giant planet formation than by photoevaporation or by grain growth. There is also an indication of global CO gas depletion in the disk, as the mass estimated from 13CO emission (\sim3.8 \times 10^(-4) M\odot) is about two orders of magnitude lower than that derived from the 1.3 mm continuum (\sim2.8 \times 10^(-2) M\odot).

We report Submillimeter Array (SMA) observations of CO (J=2--1, 3--2 and 6--5) and its isotopologues (13CO J=2--1, C18O J=2--1 and C17O J=3--2) in the disk around the Herbig Ae star HD 163296 at ~2" (250 AU) resolution, and interpret these data in the framework of a model that constrains the radial and vertical location of the line emission regions. First, we develop a physically self-consistent accretion disk model with an exponentially tapered edge that matches the spectral energy distribution and spatially resolved millimeter dust continuum emission. Then, we refine the vertical structure of the model using wide range of excitation conditions sampled by the CO lines, in particular the rarely observed J=6--5 transition. By fitting 13CO data in this structure, we further constrain the vertical distribution of CO to lie between a lower boundary below which CO freezes out onto dust grains (T ~ 19 K) and an upper boundary above which CO can be photodissociated (the hydrogen column density from the disk surface is ~ 10^{21} cm-2). The freeze-out at 19 K leads to a significant drop in the gas-phase CO column density beyond a radius of ~155 AU, a "CO snow line" that we directly resolve. By fitting the abundances of all CO isotopologues, we derive isotopic ratios of 12C/13C, 16O/18O and 18O/17O that are consistent with quiescent interstellar gas-phase values. This detailed model of the HD 163296 disk demonstrates the potential of a staged, parametric technique for constructing unified gas and dust structure models and constraining the distribution of molecular abundances using resolved multi-transition, multi-isotope observations.

We present 880 um Submillimeter Array observations of the debris disks around the young solar analogue HD 107146 and the multiple-planet host star HR 8799, at an angular resolution of 3" and 6", respectively. We spatially resolve the inner edge of the disk around HR 8799 for the first time. While the data are not sensitive enough (with rms noise of 1 mJy) to constrain the system geometry, we demonstrate that a model by Su et al. (2009) based on the spectral energy distribution (SED) with an inner radius of 150 AU predicts well the spatially resolved data. Furthermore, by modeling simultaneously the SED and visibilities, we demonstrate that the dust is distributed in a broad (of order 100 AU) annulus rather than a narrow ring. We also model the observed SED and visibilities for the HD 107146 debris disk and generate a model of the dust emission that extends in a broad band between 50 and 170 AU from the star. We perform an a posteriori comparison with existing 1.3 mm CARMA observations and demonstrate that a smooth, axisymmetric model reproduces well all of the available millimeter-wavelength data.

This is the second in a series of papers based on data from DISCS, a Submillimeter Array observing program aimed at spatially and spectrally resolving the chemical composition of 12 protoplanetary disks. We present data on six Southern sky sources - IM Lup, SAO 206462 (HD 135344b), HD 142527, AS 209, AS 205 and V4046 Sgr - which complement the six sources in the Taurus star forming region reported previously. CO 2-1 and HCO+ 3-2 emission are detected and resolved in all disks and show velocity patterns consistent with Keplerian rotation. Where detected, the emission from DCO+ 3-2, N2H+ 3-2, H2CO 3-2 and 4-3,HCN 3-2 and CN 2-1 are also generally spatially resolved. The detection rates are highest toward the M and K stars, while the F star SAO 206462 has only weak CN and HCN emission, and H2CO alone is detected toward HD 142527. These findings together with the statistics from the previous Taurus disks, support the hypothesis that high detection rates of many small molecules depend on the presence of a cold and protected disk midplane, which is less common around F and A stars compared to M and K stars. Disk-averaged variations in the proposed radiation tracer CN/HCN are found to be small, despite two orders of magnitude range of spectral types and accretion rates. In contrast, the resolved images suggest that the CN/HCN emission ratio varies with disk radius in at least two of the systems. There are no clear observational differences in the disk chemistry between the classical/full T Tauri disks and transitional disks. Furthermore, the observed line emission does not depend on measured accretion luminosities or the number of infrared lines detected, which suggests that the chemistry outside of 100 AU is not coupled to the physical processes that drive the chemistry in the innermost few AU.

Circumstellar disks are thought to experience a rapid "transition" phase in their evolution that can have a considerable impact on the formation and early development of planetary systems. We present new and archival high angular resolution (0.3" = 40-75 AU) Submillimeter Array (SMA) observations of the 880 micron dust continuum emission from 12 such transition disks in nearby star-forming regions. In each case, we directly resolve a dust-depleted disk cavity around the central star. Using radiative transfer calculations, we interpret these dust disk structures in a homogeneous, parametric model framework by reproducing their SMA visibilities and SEDs. The cavities in these disks are large (R_cav = 15-73 AU) and substantially depleted of small (~um-sized) dust grains, although their mass contents are still uncertain. The structures of the remnant material at larger radii are comparable to normal disks. We demonstrate that these large cavities are common among the millimeter-bright disk population, comprising at least 20% of the disks in the bright half of the millimeter luminosity (disk mass) distribution. Utilizing these results, we assess some of the physical mechanisms proposed to account for transition disk structures. As has been shown before, photoevaporation models do not produce the large cavity sizes, accretion rates, and disk masses representative of this sample. It would be difficult to achieve a sufficient decrease of the dust optical depths in these cavities by particle growth alone: substantial growth (to meter sizes or beyond) must occur in large (tens of AU) regions of low turbulence without also producing an abundance of small particles. Given those challenges, we suggest instead that the observations are most commensurate with dynamical clearing due to tidal interactions with low-mass companions --young brown dwarfs or giant planets on long-period orbits.

We present arcsecond-scale Submillimeter Array observations of the CO(3-2) line emission from the disks around the young stars HD 163296 and TW Hya at a spectral resolution of 44 m/s. These observations probe below the ~100 m/s turbulent linewidth inferred from lower-resolution observations, and allow us to place constraints on the turbulent linewidth in the disk atmospheres. We reproduce the observed CO(3-2) emission using two physical models of disk structure: (1) a power-law temperature distribution with a tapered density distribution following a simple functional form for an evolving accretion disk, and (2) the radiative transfer models developed by D'Alessio et al. that can reproduce the dust emission probed by the spectral energy distribution. Both types of models yield a low upper limit on the turbulent linewidth (Doppler b-parameter) in the TW Hya system (<40 m/s), and a tentative (3-sigma) detection of a ~300 m/s turbulent linewidth in the upper layers of the HD 163296 disk. These correspond to roughly <10% and 40% of the sound speed at size scales commensurate with the resolution of the data. The derived linewidths imply a turbulent viscosity coefficient, alpha, of order 0.01 and provide observational support for theoretical predictions of subsonic turbulence in protoplanetary accretion disks.

We present new results from a significant extension of our previous high angular resolution (0.3" = 40 AU) Submillimeter Array survey of the 880 um continuum emission from dusty circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. An expanded sample is constructed to probe disk structures that emit significantly lower millimeter luminosities (hence dust masses), down to the median value for T Tauri stars. Using a Monte Carlo radiative transfer code, the millimeter visibilities and broadband spectral energy distribution for each disk are simultaneously reproduced with a two-dimensional parametric model for a viscous accretion disk. We find wide ranges of characteristic radii (14-198 AU) and disk masses (0.004-0.143 M_sun), but a narrow distribution of surface density gradients (0.4-1.1) that is consistent with a uniform value $\gamma$ = 0.9 +/- 0.2 and independent of mass (or millimeter luminosity). In this sample, we find a correlation between the disk luminosity/mass and characteristic radius, such that fainter disks are both smaller and less massive. We suggest that this relationship is an imprint of the initial conditions inherited by the disks at their formation epoch, compare their angular momenta with those of molecular cloud cores, and speculate on how future observations can help constrain the distribution of viscous evolution timescales. No other correlations between disk and star properties are found. The inferred disk structures are briefly compared with theoretical models for giant planet formation, although resolution limitations do not permit us to directly comment on material inside R = 20 AU. However, there is some compelling evidence for dust evolution in the planet formation region: 4/17 disks in the sample show resolved regions of significantly reduced optical depths within ~20-40 AU of their central stars.

We have imaged the disk surrounding the nearby (D~73 pc), ~12 Myr, classical T Tauri binary system V4046 Sgr with the Submillimeter Array (SMA) at an angular resolution of ~2". We detect a rotating disk in 12CO(2-1) and 13CO(2-1) emission, and resolve the continuum emission at 1.3 mm. We infer disk gas and dust masses of ~110 and ~40 Earth masses, respectively. Fits to a power-law disk model indicate that the molecular disk extends to ~370 AU and is viewed at an inclination of between ~33 and ~39 degrees for dynamical stellar masses ranging from 1.8 $M_\odot$ down to 1.5 $M_\odot$ (the range of total mass previously determined for the central, 2.4 day spectroscopic binary). This range of disk inclination is consistent with that assumed in deducing the central binary mass (i.e., 35 degrees), suggesting that the V4046 Sgr binary system and its circumbinary, molecular disk are coplanar. In light of the system's age and binarity, the presence of an extensive molecular disk orbiting V4046 Sgr provides constraints on the timescales of processes related to Jovian planet formation, and demonstrates that circumbinary Jovian planets potentially could form around close binary systems.

The late stages of evolution of the primordial circumstellar disks surrounding young stars are poorly understood, yet vital to constrain theories of planet formation. We consider basic structural models for the disks around two ~10 Myr-old members of the nearby RCrA association, RX J1842.9-3532 and RX J1852.3-3700. We present new arcsecond-resolution maps of their 230 GHz continuum emission from the Submillimeter Array and unresolved CO(3-2) spectra from the Atacama Submillimeter Telescope Experiment. By combining these data with broadband fluxes from the literature and infrared fluxes and spectra from the catalog of the Formation and Evolution of Planetary Systems (FEPS) Legacy program on the Spitzer Space Telescope, we assemble a multiwavelength data set probing the gas and dust disks. Using the Monte Carlo radiative transfer code RADMC to model simultaneously the SED and millimeter continuum visibilities, we derive basic dust disk properties and identify an inner cavity of radius 16 AU in the disk around RX J1852.3-3700. We also identify an optically thin 5 AU cavity in the disk around RX J1842.9-3532, with a small amount of optically thick material close to the star. The molecular line observations suggest an intermediate disk inclination in RX J1842.9-3532, consistent with the continuum emission. In combination with the dust models, the molecular data allow us to derive a lower CO content than expected, suggesting that the process of gas clearing is likely underway in both systems, perhaps simultaneously with planet formation.

Chemistry plays an important role in the structure and evolution of protoplanetary disks, with implications for the composition of comets and planets. This is the first of a series of papers based on data from DISCS, a Submillimeter Array survey of the chemical composition of protoplanetary disks. The six Taurus sources in the program (DM Tau, AA Tau, LkCa 15, GM Aur, CQ Tau and MWC 480) range in stellar spectral type from M1 to A4 and offer an opportunity to test the effects of stellar luminosity on the disk chemistry. The disks were observed in 10 different lines at ~3" resolution and an rms of ~100 mJy beam-1 at ~0.5 km s-1. The four brightest lines are CO 2-1, HCO+ 3-2, CN 2_3-1_2 and HCN 3-2 and these are detected toward all sources (except for HCN toward CQ Tau). The weaker lines of CN 2_2-1_1, DCO+ 3-2, N2H+ 3-2, H2CO 3_03-2_02 and 4_14-3_13 are detected toward two to three disks each, and DCN 3-2 only toward LkCa 15. CH3OH 4_21-3_12 and c-C3H2 are not detected. There is no obvious difference between the T Tauri and Herbig Ae sources with regard to CN and HCN intensities. In contrast, DCO+, DCN, N2H+ and H2CO are detected only toward the T Tauri stars, suggesting that the disks around Herbig Ae stars lack cold regions for long enough timescales to allow for efficient deuterium chemistry, CO freeze-out, and grain chemistry.

Mid-infrared spectrophotometric observations have revealed a small sub-class of circumstellar disks with spectral energy distributions (SEDs) suggestive of large inner gaps with low dust content. However, such data provide only an indirect and model-dependent method of finding central holes. Imaging of protoplanetry disks provides an independent check of SED modeling. We present here the direct characterization of three 33-47 AU radii inner gaps, in the disks around LkHa 330, SR 21N and HD 135344B, via 340 GHz (880 micron) dust continuum aperture synthesis observations obtained with the Submillimeter Array (SMA). The large gaps are fully resolved at ~0\farcs3 by the SMA observations and mostly empty of dust, with less than 1 - 7.5 x 10^-6 Msolar of fine grained solids inside the holes. Gas (as traced by atomic accretion markers and CO 4.7 micron rovibrational emission) is still present in the inner regions of all three disks. For each, the inner hole exhibits a relatively steep rise in dust emission to the outer disk, a feature more likely to originate from the gravitational influence of a companion body than from a process expected to show a more shallow gradient like grain growth. Importantly, the good agreement of the spatially resolved data and spectrophotometry-based models lends confidence to current interpretations of SEDs, wherein the significant dust emission deficits arise from disks with inner gaps or holes. Further SED-based searches can therefore be expected to yield numerous additional candidates that can be examined at high spatial resolution.

We present a high angular resolution (0.3" = 40 AU) SMA survey of the 870 micron thermal continuum emission from 9 of the brightest, and therefore most massive, circumstellar disks in the ~1 Myr-old Ophiuchus star-forming region. Using 2-D radiative transfer calculations, we simultaneously fit the observed continuum visibilities and broadband spectral energy distribution for each disk with a parametric structure model. Compared to previous millimeter studies, this survey includes significant upgrades in modeling, data quality, and angular resolution that provide improved constraints on key structure parameters, particularly those that characterize the spatial distribution of mass in the disks. In the context of a surface density profile motivated by similarity solutions for viscous accretion disks, the best-fit models for the sample disks have characteristic radii R_c = 20-200 AU, high disk masses M_d = 0.005-0.14 M_sun, and a narrow range of radial surface density gradients around a median $\gamma$ = 0.9. These density structures are used in conjunction with accretion rate estimates from the literature to help characterize the viscous evolution of the disk material. Using the standard prescription for disk viscosities, those combined constraints indicate that $\alpha$ = 0.0005-0.08. Three of the sample disks show large (R = 20-40 AU) central cavities in their continuum emission morphologies, marking extensive zones where dust has been physically removed and/or has significantly diminished opacities. Based on the current requirements of planet formation models, these emission cavities and the structure constraints for the sample as a whole suggest that these young disks may eventually produce planetary systems, and have perhaps already started. (abridged)

We study the distribution and physical properties of molecular gas in the disc around the T Tauri star IM Lup on scales close to 200 AU. We investigate how well the gas and dust distributions compare and work towards a unified disc model that can explain both gas and dust emission. 12CO, 13CO, and C18O J=2-1 line emission, as well as the dust continuum at 1.3 mm, is observed at 1.8" resolution towards IM Lup using the Submillimeter Array. A detailed disc model based on the dust emission is tested against these observations with the aid of a molecular excitation and radiative transfer code. Apparent discrepancies between the gas and dust distribution are investigated by adopting simple modifications to the existing model. The disc is seen at an inclination of 54+/-3 degrees and is in Keplerian rotation around a 0.8-1.6 Msun star. The outer disc radius traced by molecular gas emission is 900 AU, while the dust continuum emission and scattered light images limit the amount of dust present beyond 400 AU and are consistent with the existing model that assumes a 400 AU radius. Our observations require a drastic density decrease close to 400 AU with the vertical gas column density at 900 AU in the range of 5.d20 - 1.d22 cm-2. We derive a gas-to-dust mass ratio of 100 or higher in disc regions beyond 400 AU. Within 400 AU from the star our observations are consistent with a gas-to-dust ratio of 100 but other values are not ruled out.

We present 0.3 arcsec resolution observations of the disk around GM Aurigae with the Submillimeter Array (SMA) at a wavelength of 860 um and with the Plateau de Bure Interferometer at a wavelength of 1.3 mm. These observations probe the distribution of disk material on spatial scales commensurate with the size of the inner hole predicted by models of the spectral energy distribution. The data clearly indicate a sharp decrease in millimeter optical depth at the disk center, consistent with a deficit of material at distances less than ~20 AU from the star. We refine the accretion disk model of Calvet et al. (2005) based on the unresolved spectral energy distribution (SED) and demonstrate that it reproduces well the spatially resolved millimeter continuum data at both available wavelengths. We also present complementary SMA observations of CO J=3-2 and J=2-1 emission from the disk at 2" resolution. The observed CO morphology is consistent with the continuum model prediction, with two significant deviations: (1) the emission displays a larger CO J=3-2/J=2-1 line ratio than predicted, which may indicate additional heating of gas in the upper disk layers; and (2) the position angle of the kinematic rotation pattern differs by 11 +/- 2 degrees from that measured at smaller scales from the dust continuum, which may indicate the presence of a warp. We note that photoevaporation, grain growth, and binarity are unlikely mechanisms for inducing the observed sharp decrease in opacity or surface density at the disk center. The inner hole plausibly results from the dynamical influence of a planet on the disk material. Warping induced by a planet could also potentially explain the difference in position angle between the continuum and CO data sets.

We investigate the physical structure of the gas component of the disk around the pre-main-sequence star HD169142. The 13CO and C18O J=2-1 line emission is observed from the disk with 1.4'' resolution using the Submillimeter Array. We adopt the disk physical structure derived from a model which fits the spectral energy distribution of HD169142. We obtain the full three-dimensional information on the CO emission with the aid of a molecular excitation and radiative transfer code. This information is used for the analysis of our observations and previous 12CO J=2-1 and 1.3 mm continuum data. The disk is in Keplerian rotation and seen at an inclination close to 13 deg from face-on. We conclude that the regions traced by different CO isotopologues are distinct in terms of their vertical location within the disk, their temperature and their column densities. With the given disk structure, we find that freeze-out is not efficient enough to remove a significant amount of CO from gas phase. Both observed lines match the model prediction both in flux and in the spatial structure of the emission. Therefore we use our data to derive the 13CO and C18O mass and consequently the 12CO mass using standard isotopic ratios. We constrain the total disk gas mass to (0.6-3.0)x10(-2) Msun. Adopting a maximum dust opacity of 2 cm2 per gram of dust we derive a minimum dust mass of 2.16x10(-4) Msun from the fit to the 1.3 mm data. Comparison of the derived gas and dust mass shows that the gas to dust mass ratio of 100 is only possible under the assumption of a dust opacity of 2 cm2/g and 12CO abundance of 10(-4) with respect to H2. However, our data are also compatible with a gas to dust ratio of 25, with a dust opacity of 1 cm2/g and 12CO abundance of 2x10(-4).

We present high spatial resolution (< 0.3" = 40$ AU) Submillimeter Array observations of the 865 micron continuum emission from the circumstellar disk around the young star DoAr 25. Despite its bright millimeter emission, this source exhibits only a comparatively small infrared excess and low accretion rate, suggesting that the material and structural properties of the inner disk may be in an advanced state of evolution. A simple model of the physical conditions in the disk is derived from the submillimeter visibilities and the complete spectral energy distribution using a Monte Carlo radiative transfer code. For the standard assumption of a homogeneous grain size distribution at all disk radii, the results indicate a shallow surface density profile, $\Sigma \propto r^{-p}$ with p = 0.34, significantly less steep than a steady-state accretion disk (p = 1) or the often adopted minimum mass solar nebula (p = 1.5). Even though the total mass of material is large (M_d = 0.10 M_sun), the densities inferred in the inner disk for such a model may be too low to facilitate any mode of planet formation. However, alternative models with steeper density gradients (p = 1) can explain the observations equally well if substantial grain growth in the planet formation region (r < 40 AU) has occurred. We discuss these data in the context of such models with dust properties that vary with radius and highlight their implications for understanding disk evolution and the early stages of planet formation.

[Abridged] We imaged a 2' x 2' region of the Orion Nebula cluster in 1.3 mm wavelength continuum emission with the recently commissioned Combined Array for Research in Millimeter Astronomy (CARMA) and with the Submillimeter Array (SMA). Our mosaics include >250 known near-IR cluster members, of which 36 are so-called "proplyds" that have been imaged previously with the Hubble Space Telescope. We detected 40 sources in 1 mm continuum emission, and several of them are spatially resolved with our observations. Dust masses inferred for detected sources range from 0.01 to 0.5 Msun, and the average disk mass for undetected sources is estimated to be ~0.001 Msun, approximately an order of magnitude smaller than the minimum mass solar nebula. Most stars in the ONC thus do not appear to currently possess sufficient mass in small dust grains to form Jupiter-mass (or larger) planets. Comparison with previous results for younger and older regions indicates that massive disks evolve significantly on ~Myr timescales. We also show that the percentage of stars in Orion surrounded by disks more massive than ~0.01 Msun is substantially lower than in Taurus, indicating that environment has an impact on the disk mass distribution. Finally, we explore potential correlations of disk mass with stellar mass and location within the cluster.

We present Submillimeter Array (SMA) observations of several deuterated species in the disk around the classical T Tauri star TW Hydrae at arcsecond scales, including detections of the DCN J=3-2 and DCO+ J=3-2 lines, and upper limits to the HDO 3(1,2)-2(2,1), ortho-H2D+ 1(1,0)-1(1,1) and para-D2H+ 1(1,0)-1(0,1) transitions. We also present observations of the HCN J=3-2, HCO+ J=3-2 and H13CO+ J=4-3 lines for comparison with their deuterated isotopologues. We constrain the radial and vertical distributions of various species in the disk by fitting the data using a model where the molecular emission from an irradiated accretion disk is sampled with a 2D Monte Carlo radiative transfer code. We find that the distribution of DCO+ differs markedly from that of HCO+. The D/H ratios inferred change by at least one order of magnitude (0.01 to 0.1) for radii <30 AU to >70 AU and there is a rapid falloff of the abundance of DCO+ at radii larger than 90 AU. Using a simple analytical chemical model, we constrain the degree of ionization, x(e-)=n(e-)/n(H2), to be ~10^-7 in the disk layer(s) where these molecules are present. Provided the distribution of DCN follows that of HCN, the ratio of DCN to HCN is determined to be 1.7\pm0.5 \times 10^-2; however, this ratio is very sensitive to the poorly constrained vertical distribution of HCN. The resolved radial distribution of DCO+ indicates that {\it in situ} deuterium fractionation remains active within the TW Hydrae disk and must be considered in the molecular evolution of circumstellar accretion disks.

We investigate the apparent discrepancy between gas and dust outer radii derived from millimeter observations of protoplanetary disks. Using 230 and 345 GHz continuum and CO J=3-2 data from the Submillimeter Array for four nearby disk systems (HD 163296, TW Hydrae, GM Aurigae, and MWC 480), we examine models of circumstellar disk structure and the effects of their treatment of the outer disk edge. We show that for these disks, models described by power laws in surface density and temperature that are truncated at an outer radius are incapable of reproducing both the gas and dust emission simultaneously: the outer radius derived from the dust continuum emission is always significantly smaller than the extent of the molecular gas disk traced by CO emission. However, a simple model motivated by similarity solutions of the time evolution of accretion disks that includes a tapered exponential edge in the surface density distribution (and the same number of free parameters) does much better at reproducing both the gas and dust emission. While this analysis does not rule out the disparate radii implied by the truncated power-law models, a realistic alternative disk model, grounded in the physics of accretion, provides a consistent picture for the extent of both the gas and dust.

Mid-infrared spectrophotometric observations have revealed a small sub-class of circumstellar disks with spectral energy distributions (SEDs) suggestive of large inner gaps with low dust content. However, such data provide only an indirect and model dependent method of finding central holes. We present here the direct characterization of a 40 AU radius inner gap in the disk around LkHa 330 through 340 GHz (880 micron) dust continuum imaging with the Submillimeter Array (SMA). This large gap is fully resolved by the SMA observations and mostly empty of dust with less than 1.3 x 10^-6 M_solar of solid particles inside of 40 AU. Gas (as traced by accretion markers and CO M-band emission) is still present in the inner disk and the outer edge of the gap rises steeply -- features in better agreement with the underlying cause being gravitational perturbation than a more gradual process such as grain growth. Importantly, the good agreement of the spatially resolved data and spectrophometry-based model lends confidence to current interpretations of SEDs with significant dust emission deficits as arising from disks with inner gaps or holes. Further SED-based searches can therefore be expected to yield numerous additional candidates that can be examined at high spatial resolution.

The recently recognized class of "transitional disk" systems consists of young starswith optically-thick outer disks but inner disks which are mostly devoid of small dust. Here we introduce a further class of "pre-transitional disks" with significant near-infrared excesses which indicate the presence of an optically thick inner disk separated from an optically thick outer disk; thus, the spectral energy distributions of pre-transitional disks suggest the incipient development of disk gaps rather than inner holes. In UX Tau A, our analysis of the Spitzer IRS spectrum finds that the near-infrared excess is produced by an inner optically thick disk and a gap of ~56 AU is present. The Spitzer IRS spectrum of LkCa 15 is suggestive of a gap of ~46 AU, confirming previous millimeter imaging. In addition, UX Tau A contains crystalline silicates in its disk at radii >~ 56 AU which poses a challenge to our understanding of the production of this crystalline material. In contrast, LkCa 15's silicates are amorphous and pristine. UX Tau A and LkCa 15 increase our knowledge of the diversity of dust clearing in low-mass star formation.

We present new multi-wavelength millimeter interferometric observations of the Herbig Ae star HD 163296 obtained with the IRAM/PBI, SMA and VLA arrays both in continuum and in the 12CO, 13CO and C18O emission lines. Gas and dust properties have been obtained comparing the observations with self-consistent disk models for the dust and CO emission. The circumstellar disk is resolved both in the continuum and in CO. We find strong evidence that the circumstellar material is in Keplerian rotation around a central star of 2.6 Msun. The disk inclination with respect to the line of sight is 46+-4 deg with a position angle of 128+-4 deg. The slope of the dust opacity measured between 0.87 and 7 mm (beta=1) confirms the presence of mm/cm-size grains in the disk midplane. The dust continuum emission is asymmetric and confined inside a radius of 200 AU while the CO emission extends up to 540 AU. The comparison between dust and CO temperature indicates that CO is present only in the disk interior. Finally, we obtain an increasing depletion of CO isotopomers from 12CO to 13CO and C18O. We argue that these results support the idea that the disk of HD 163296 is strongly evolved. In particular, we suggest that there is a strong depletion of dust relative to gas outside 200 AU; this may be due to the inward migration of large bodies that form in the outer disk or to clearing of a large gap in the dust distribution by a low mass companion.

We present Submillimeter Array observations of the Herbig Ae star HD169142 in 1.3 millimeter continuum emission and 12CO J=2-1 line emission at 1.5 arcsecond resolution that reveal a circumstellar disk. The continuum emission is centered on the star position and resolved, and provides a mass estimate of about 0.02 solar masses for the disk. The CO images show patterns in position and velocity that are well matched by a disk in Keplerian rotation with low inclination to the line-of-sight. We use radiative transfer calculations based on a flared, passive disk model to constrain the disk parameters by comparison to the spectral line emission. The derived disk radius is 235 AU, and the inclination is 13 degrees. The model also necessitates modest depletion of the CO molecules, similar to that found in Keplerian disks around T Tauri stars.

We present the first images of the 691.473 GHz CO J=6-5 line in a protoplanetary disk, obtained along with the 690 GHz dust continuum, toward the classical T Tauri star TW Hya using the Submillimeter Array. Imaging in the CO J=6-5 line reveals a rotating disk, consistent with previous observations of CO J=3-2 and 2-1 lines. Using an irradiated accretion disk model and 2D Monte Carlo radiative transfer, we find that additional surface heating is needed to fit simultaneously the absolute and relative intensities of the CO J=6-5, 3-2 and 2-1 lines. In particular, the vertical gas temperature gradient in the disk must be steeper than that of the dust, mostly likely because the CO emission lines probe nearer to the surface of the disk. We have used an idealized X-ray heating model to fit the line profiles of CO J=2-1 and 3-2 with Chi-square analysis, and the prediction of this model yields CO J=6-5 emission consistent with the observations.

We present high angular resolution (~ 1''), sensitive CO(2-1) line observations of the region OMC1 South in the Orion Nebula made using the Submillimeter Array (SMA). We detect the CO(2-1) high velocity outflow that was first found by Rodriguez-Franco et al. (1999a) with the IRAM 30 m. Our observations resolve the outflow, whose velocity-integrated emission has a deconvolved width of 0.89'' \pm 0.06'' (490 AU) and a projected length of ~ 48'' (21,000 AU) with very high redshifted and blueshifted gas with velocities of about \pm 80 km/s. This outflow is among the most collimated (~ 3 degrees) and youngest outflows (600 yr) that have been reported. The data show that this collimated outflow has been blowing in the same direction during the last 600 yr. At high velocities, the CO(2-1) outflow traces an extremely collimated jet, while at lower velocities the CO emission traces an envelope possibly produced by entrainment of ambient gas. Furthermore, we also detect for the first time a millimeter wavelength continuum source possibly associated with a class I protostar that we suggest could be the possible exciting source for this collimated outflow. However, the bolometric luminosity of this source appears to be far too low to account for the powerful molecular outflow.

Arcsecond-resolution spectral observations toward the protobinary system IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillimeter Array. Complex organic molecules such as CH3OH and HCOOCH3 were detected. Together with the rich organic inventory revealed, it clearly indicates the existence of two, rather than one, compact hot molecular cores (smaller than or equal to 400 AU in radius) associated with each of the protobinary components identified by their dust continuum emission in the inner star-forming core.

We present new observations of the FUV (1100-2200 Angstrom) radiation field and the near- to mid-IR (3--13.5 micron) spectral energy distribution (SED) of a sample of T Tauri stars selected on the basis of bright molecular disks (GM Aur, DM Tau, LkCa15). In each source we find evidence for Ly alpha induced H2 fluorescence and an additional source of FUV continuum emission below 1700 Angstroms. Comparison of the FUV spectra to a model of H2 excitation suggests that the strong continuum emission is due to electron impact excitation of H2. The ultimate source of this excitation is likely X-ray irradiation which creates hot photo-electrons mixed in the molecular layer. Analysis of the SED of each object finds the presence of inner disk gaps with sizes of a few AU in each of these young (~1 Myr) stellar systems. We propose that the presence of strong H2 continuum emission and inner disk clearing are related by the increased penetration power of high energy photons in gas rich regions with low grain opacity.

Using the partially completed Submillimeter Array with five antennas, we have observed the CO J=2-1 and 3-2 emission from the envelope surrounding the carbon star V Hya. The high-angular resolution (2"-4") maps show that V Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 km/s. The axis of this high velocity jet is perpendicular to the major axis of the flattened disk-like envelope, which is expanding with a velocity of ~16 km/s.There is a third kinematic component, a medium-velocity wind having a deprojected velocity of 40-120 km/s moving along the disk plane. Both the high velocity jet and the medium velocity wind have a dynamical time scale of a few hundred years. The flattened structure and the collimated jet observed in V Hya suggests that the formation of asymmetrical structure proceeds while the central star is still in the AGB phase

We have imaged the circumstellar envelope around the binary protostar L1551 IRS 5 in CS (J=7-6) and 343 GHz continuum emission at ~ 3 arcsec resolution using the Submillimeter Array. The continuum emission shows an elongated structure (~ 220 x 100 AU) around the binary perpendicular to the axis of the associated radio jet. The CS emission extends over ~ 400 AU, appears approximately circularly symmetric, and shows a velocity gradient from southeast (blueshifted) to northwest (redshifted). The direction of the velocity gradient is different from that observed in C18O (J=1-0). This may be because rotation is more dominant in the CS envelope than the C18O envelope, in which both infall and rotation exist. The CS emission may be divided into two velocity components: (1) a "high" velocity disk-like structure surrounding the protostar, +-1.0 - 1.5 km s-1 from the systemic velocity, and (2) a "low" velocity structure, located southwest of the protostar, < 1.0 km s -1 from the systemic velocity. The high-velocity component traces warm and dense gas with kinematics consistent with rotation around the protostar. The low-velocity component may arise from dense gas entrained in the outflow. Alternatively, this component may trace infalling and rotating gas in an envelope with a vertical structure.

We present ~2"-4" aperture synthesis observations of the circumstellar disk surrounding the nearby young star TW Hya in the CO J=2--1 and J=3--2 lines and associated dust continuum obtained with the partially completed Submillimeter Array. The extent and peak flux of the 230 and 345 GHz dust emission follow closely the predictions of the irradiated accretion disk model of Calvet et al. (2002). The resolved molecular line emission extends to a radius of at least 200 AU, the full extent of the disk visible in scattered light, and shows a clear pattern of Keplerian rotation. Comparison of the images with 2D Monte Carlo models constrains the disk inclination angle to 7+/-1 degrees. The CO emission is optically thick in both lines, and the kinetic temperature in the line formation region is ~20K. Substantial CO depletion, by an order of magnitude or more from canonical dark cloud values, is required to explain the characteristics of the line emission.

We present results from an observing campaign of the molecular content of the coma of comet C/1999 S4 (LINEAR) carried out jointly with the millimeter-arrays of the Berkeley-Illinois-Maryland Association (BIMA) and the Owens Valley Radio Observatory (OVRO). Using the BIMA array in autocorrelation (`single-dish') mode, we detected weak HCN J=1-0 emission from comet C/1999 S4 (LINEAR) at 14 +- 4 mK km/s averaged over the 143" beam. The three days over which emission was detected, 2000 July 21.9-24.2, immediately precede the reported full breakup of the nucleus of this comet. During this same period, we find an upper limit for HCN 1-0 of 144 mJy/beam km/s (203 mK km/s) in the 9"x12" synthesized beam of combined observations of BIMA and OVRO in cross-correlation (`imaging') mode. Together with reported values of HCN 1-0 emission in the 28" IRAM 30-meter beam, our data probe the spatial distribution of the HCN emission from radii of 1300 to 19,000 km. Using literature results of HCN excitation in cometary comae, we find that the relative line fluxes in the 12"x9", 28" and 143" beams are consistent with expectations for a nuclear source of HCN and expansion of the volatile gases and evaporating icy grains following a Haser model.

We present OVRO Millimeter Array lambda = 3.4 - 1.2 mm dust continuum and spectral line observations of the accretion disk encircling the T Tauri star LkCa 15. The 1.2 mm dust continuum emission is resolved, and gives a minimum diameter of 190 AU and an inclination angle of 57+/-5 degrees. There is a noticeable, but at present poorly constrained, decrease in the continuum spectral slope with frequency that may result from the coupled processes of grain growth and dust settling. Imaging of the fairly intense emission from the lowest rotational transitions of CO, 13CO and HCO+ reveals a rotating disk and emission extends to 750 AU and the characteristic radius of the disk is determined to be around 425 AU (HWHM) based on model fits to the CO velocity field. The disk mass derived from the CO isotopologues with ``typical'' dense cloud abundances is still nearly two orders of magnitude less than that inferred from the dust emission, which is probably due to extensive molecular depletion in the cold, dense disk midplane. N2H+ 1-0 emission has also been detected which, along with HCO+, sets a lower limit to the fractional ionization of 10^{-8} in the near-surface regions of protoplanetary disks. This first detection of N2H+ in circumstellar disks has also made possible a determination of the N2/CO ratio (~2) that is at least an order of magnitude larger than those in the envelopes of young stellar objects and dense clouds. The large N2/CO ratio indicates that our observations probe disk layers in which CO is depleted but some N2 remains in the gas phase. Such differential depletion can lead to large variations in the fractional ionization with height in the outer reaches of circumstellar disks, and may help to explain the relative nitrogen deficiency observed in comets.

Emission from the $2_{12}-1_{11}$ line of H$_2$CO has been detected and marginally resolved toward LkCa15 by the Nobeyama Millimeter Array. The column density of H$_2$CO is higher than that observed in DM Tau and than predicted by theoretical models of disk chemistry; also the line-intensity profile is less centrally peaked than that for CO. A similar behavior is observed in other organic gaseous molecules in the LkCa 15 disk.

We improve upon the radiative, hydrostatic equilibrium models of passive circumstellar disks constructed by Chiang & Goldreich (1997). New features include (1) account for a range of particle sizes, (2) employment of laboratory-based optical constants of representative grain materials, and (3) numerical solution of the equations of radiative and hydrostatic equilibrium within the original 2-layer (disk surface + disk interior) approximation. We explore how the spectral energy distribution (SED) of a face-on disk depends on grain size distributions, disk geometries and surface densities, and stellar photospheric temperatures. Observed SEDs of 3 Herbig Ae and 2 T Tauri stars, including spectra from the Long Wavelength Spectrometer (LWS) aboard the Infrared Space Observatory (ISO), are fitted with our models. Silicate emission bands from optically thin, superheated disk surface layers appear in nearly all systems. Water ice emission bands appear in LWS spectra of 2 of the coolest stars. Infrared excesses in several sources are consistent with vertical settling of photospheric grains. While this work furnishes further evidence that passive reprocessing of starlight by flared disks adequately explains the origin of infrared-to-millimeter wavelength excesses of young stars, we emphasize how the SED alone does not provide sufficient information to constrain particle sizes and disk masses uniquely.