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.

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 composition of a protoplanetary disk is set by a combination of interstellar inheritance and gas and grain surface chemical reactions within the disk. The survival of inherited molecules, as well as the disk in situ chemistry depends on the local temperature, density and irradiation environment, which can change over time due to stellar and disk evolution, as well as transport in the disk. We address one aspect of this coupling between the physical and chemical evolution in disks by following accretion streamlines of gas and small grains in the disk midplane, while simultaneously taking the evolving star into account. This approach is computationally efficient and enables us to take into account changing physical conditions without reducing the chemical network. We find that many species are enhanced in the inner disk midplane in the dynamic model due to inward transport of cosmic-ray driven chemical products, resulting in, e.g., orders-of magnitude hydrocarbon enhancements at 1 au, compared to a static disk. For several other chemical families, there is no difference between the static and dynamic models, indicative of a robust chemical reset, while yet others show differences between static and dynamic models that depend on complex interactions between physics and chemistry during the inward track. The importance of coupling dynamics and chemistry when modeling the chemical evolution of protoplanetary disks is thus depends on what chemistry is of interest.

Young stars emit strong flares of X-ray radiation that penetrate the surface layers of their associated protoplanetary disks. It is still an open question as to whether flares create significant changes in disk chemical composition. We present models of the time-evolving chemistry of gas-phase water during X-ray flaring events. The chemistry is modeled at point locations in the disk between 1 and 50 au at vertical heights ranging from the mid-plane to the surface. We find that strong, rare flares, i.e., those that increase the unattenuated X-ray ionization rate by a factor of 100 every few years, can temporarily increase the gas-phase water abundance relative to H can by more than a factor of $\sim3-5$ along the disk surface (Z/R $\ge$ 0.3). We report that a "typical" flare, i.e., those that increase the unattenuated X-ray ionization rate by a factor of a few every few weeks, will not lead to significant, observable changes. Dissociative recombination of H$_3$O$^+$, water adsorption and desorption onto dust grains, and ultraviolet photolysis of water and related species are found to be the three dominant processes regulating the gas-phase water abundance. While the changes are found to be significant, we find that the effect on gas phase water abundances throughout the disk is short-lived (days). Even though we do not see a substantial increase in long term water (gas and ice) production, the flares' large effects may be detectable as time varying inner disk water 'bursts' at radii between 5 and 30 au with future far infrared observations.

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.

CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemically reprocessed into less volatile species, which are frozen on dust grain surfaces as ice. This mechanism has the potential to change the primordial C/H ratio in the gas. However, current observations primarily probe the upper layers of the disk. It is not clear if the low volatile carbon abundances extend to the midplane, where planets form. We have run a grid of 198 chemical models, exploring how the chemical reprocessing of CO depends on disk mass, dust grain size distribution, temperature, cosmic ray and X-ray ionization rate, and initial water abundance. Building on our previous work focusing on the warm molecular layer, here we analyze the results for our grid of models in the disk midplane at 12 au. We find that either an ISM level cosmic ray ionization rate or the presence of UV photons due to a low dust surface density are needed to chemically reduce the midplane CO gas abundance by at least an order of magnitude within 1 Myr. In the majority of our models CO does not undergo substantial reprocessing by in situ chemistry and there is little change in the gas phase C/H and C/O ratios over the lifetime of the typical disk. However, in the small sub-set of disks where the disk midplane is subject to a source of ionization or photolysis, the gas phase C/O ratio increases by up to nearly 9 orders of magnitude due to conversion of CO into volatile hydrocarbons.

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

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

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.

Planets assemble in the midplanes of protoplanetary disks. The compositions of dust and gas in the disk midplane region determine the compositions of nascent planets, including their chemical hospitality to life. In this context, the distributions of volatile organic material across the planet and comet forming zones is of special interest. These are difficult to access in the disk midplane at IR and even millimeter wavelengths due to dust opacity, which can veil the midplane, low intrinsic molecular abundances due to efficient freeze-out, and, in the case of mid-sized organics, a mismatch between expected excitation temperatures and accessible line upper energy levels. At ngVLA wavelengths, the dust is optically thin, enabling observations into the planet forming disk midplane. ngVLA also has the requisite sensitivity. Using TW Hya as a case study, we show that ngVLA will be able to map out the distributions of diagnostic organics, such as CH3CN, in nearby protoplanetary disks.

One of the primary mechanisms for inferring the dynamical history of planets in our Solar System and in exoplanetary systems is through observation of elemental ratios (i.e. C/O). The ability to effectively use these observations relies critically on a robust understanding of the chemistry and evolutionary history of the observed abundances. Significant efforts have been devoted to this area from within astrochemistry circles, and these efforts should be supported going forward by the larger exoplanetary science community. In addition, the construction of a next-generation radio interferometer will be required to test many of these predictive models in situ, while simultaneously providing the resolution necessary to pinpoint the location of planets in formation.

We present new constraints on gas-phase C, N, and O abundances in the molecular layer of the IM Lup protoplanetary disk. Building on previous physical and chemical modeling of this disk, we use new ALMA observations of C$_2$H to constrain the C/O ratio in the molecular layer to be $\sim0.8$, i.e., higher than the solar value of $\sim0.54$. We use archival ALMA observations of HCN and H$^{13}$CN to show that no depletion of N is required (assuming an interstellar abundance of $7.5\times10^{-5}$ per H). These results suggest that an appreciable fraction of O is sequestered in water ice in large grains settled to the disk mid-plane. Similarly, a fraction of the available C is locked up in less volatile molecules. By contrast, N remains largely unprocessed, likely as N$_2$. This pattern of depletion suggests the presence of true abundance variations in this disk, and not a simple overall depletion of gas mass. If these results hold more generally, then combined CO, C$_2$H, and HCN observations of disks may provide a promising path for constraining gas-phase C/O and N/O during planet-formation. Together, these tracers offer the opportunity to link the volatile compositions of disks to the atmospheres of planets formed from them.

In this chapter we outline some of the basic understanding of the chemistry that accompanies planet formation. We discuss the basic physical environment which dictates the dominant chemical kinetic pathways for molecule formation. We focus on three zones from both observational and theoretical perspectives: (1) the planet forming midplane and ice/vapor transition zones (snow-lines), (2) the warm disk surface that is shielded from radiation, which can be readily accessed by todays observational facilities, and (3) the surface photodissociation layers where stellar radiation dominates. We end with a discussion of how chemistry influences planet formation along with how to probe the link between formation and ultimate atmospheric composition for gas giants and terrestrial worlds.

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.

Cometary studies suggest that the organic composition of the early Solar Nebula was rich in complex nitrile species such a CH$_3$CN. Recent ALMA detections in protoplanetary disks suggest that these species may be common during planet and comet formation, but connecting gas phase measurements to cometary abundances first requires constraints on formation chemistry and distributions of these species. We present here the detection of seven spatially resolved transitions of CH$_3$CN in the protoplanetary disk around the T-Tauri star TW Hya. Using a rotational diagram analysis we find a disk-averaged column density of N$_T$=1.45$^{+0.19}_{-0.15}\times10^{12}$ cm$^{-2}$ and a rotational temperature of T$_{rot}$=32.7$^{+3.9}_{-3.4}$ K. A radially resolved rotational diagram shows the rotational temperature to be constant across the disk, suggesting that the CH$_3$CN emission originates from a layer at z/r$\sim$0.3. Through comparison of the observations with predictions from a disk chemistry model, we find that grain-surface reactions likely dominate CH$_3$CN formation and that in situ disk chemistry is sufficient to explain the observed CH$_3$CN column density profile without invoking inheritance from the protostellar phase. However, the same model fails to reproduce a Solar System cometary abundance of CH$_3$CN relative to H$_2$O in the midplane, suggesting that either vigorous vertical mixing or some degree of inheritance from interstellar ices occurred in the Solar Nebula.

We present 870 $\mu$m ALMA observations of polarized dust emission toward the Class II protoplanetary disk IM Lup. We find that the orientation of the polarized emission is along the minor axis of the disk, and that the value of the polarization fraction increases steadily toward the center of the disk, reaching a peak value of ~1.1%. All of these characteristics are consistent with models of self-scattering of submillimeter-wave emission from an optically thin inclined disk. The distribution of the polarization position angles across the disk reveals that while the average orientation is along the minor axis, the polarization orientations show a significant spread in angles; this can also be explained by models of pure scattering. We compare the polarization with that of the Class I/II source HL Tau. A comparison of cuts of the polarization fraction across the major and minor axes of both sources reveals that IM Lup has a substantially higher polarization fraction than HL Tau toward the center of the disk. This enhanced polarization fraction could be due a number of factors, including higher optical depth in HL Tau, or scattering by larger dust grains in the more evolved IM Lup disk. However, models yield similar maximum grain sizes for both HL Tau (72 $\mu$m) and IM Lup (61 $\mu$m, this work). This reveals continued tension between grain-size estimates from scattering models and from models of the dust emission spectrum, which find that the bulk of the (unpolarized) emission in disks is most likely due to millimeter (or even centimeter) sized grains.

CO is commonly used as a tracer of the total gas mass in both the interstellar medium and in protoplanetary disks. Recently there has been much debate about the utility of CO as a mass tracer in disks. Observations of CO in protoplanetary disks reveal a range of CO abundances, with measurements of low CO to dust mass ratios in numerous systems. One possibility is that carbon is removed from CO via chemistry. However, the full range of physical conditions conducive to this chemical reprocessing is not well understood. We perform a systematic survey of the time dependent chemistry in protoplanetary disks for 198 models with a range of physical conditions. We varying dust grain size distribution, temperature, comic ray and X-ray ionization rate, disk mass, and initial water abundance, detailing what physical conditions are necessary to activate the various CO depletion mechanisms in the warm molecular layer. We focus our analysis on the warm molecular layer in two regions: the outer disk (100 au) well outside the CO snowline and the inner disk (19 au) just inside the midplane CO snow line. After 1 Myr, we find that the majority of models have a CO abundance relative to H$_2$ less than $10^{-4}$ in the outer disk, while an abundance less than $10^{-5}$ requires the presence of cosmic rays. Inside the CO snow line, significant depletion of CO only occurs in models with a high cosmic ray rate. If cosmic rays are not present in young disks it is difficult to chemically remove carbon from CO. Additionally, removing water prior to CO depletion impedes the chemical processing of CO. Chemical processing alone cannot explain current observations of low CO abundances. Other mechanisms must also be involved.

We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the $^{12}$CO $J = 3-2$ line at a spatial resolution of 8 AU (0\farcs14). The CO emission exhibits a bright inner core, a shoulder at $r\approx70$ AU, and a prominent break in slope at $r\approx90$ AU. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a $^{12}$CO column density profile featuring a steep decrease at $r\approx15$ AU and a secondary bump peaking at $r\approx70$ AU. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disk's vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 $\mu$m continuum observations to produce a spectral index map at a spatial resolution of 2 AU. The spectral index rises sharply at the continuum emission gaps at radii of 25, 41, and 47 AU. This behavior suggests that the grains within the gaps are no larger than a few millimeters. Outside the continuum gaps, the low spectral index values of $\alpha\approx 2$ indicate either that grains up to centimeter size are present, or that the bright continuum rings are marginally optically thick at millimeter wavelengths.

We report the first detection of a substantial brightening event in an isotopologue of a key molecular ion, HCO$^+$, within a protoplanetary disk of a T Tauri star. The H$^{13}$CO$^+$ $J=3-2$ rotational transition was observed three times toward IM Lup between July 2014 and May 2015 with the Atacama Large Millimeter Array. The first two observations show similar spectrally integrated line and continuum fluxes, while the third observation shows a doubling in the disk integrated $J=3-2$ line flux compared to the continuum, which does not change between the three epochs. We explore models of an X-ray active star irradiating the disk via stellar flares, and find that the optically thin H$^{13}$CO$^+$ emission variation can potentially be explained via X-ray driven chemistry temporarily enhancing the HCO$^+$ abundance in the upper layers of the disk atmosphere during large or prolonged flaring events. If the HCO$^+$ enhancement is indeed caused by a X-ray flare, future observations should be able to spatially resolve these events and potentially enable us to watch the chemical aftermath of the high-energy stellar radiation propagating across the face of protoplanetary disks, providing a new pathway to explore ionization physics and chemistry, including electron density, in disks.

The initial mass distribution in the solar nebula is a critical input to planet formation models that seek to reproduce today's Solar System. Traditionally, constraints on the gas mass distribution are derived from observations of the dust emission from disks, but this approach suffers from large uncertainties in grain growth and gas-to-dust ratio. On the other hand, previous observations of gas tracers only probe surface layers above the bulk mass reservoir. Here we present the first partially spatially resolved observations of the $^{13}$C$^{18}$O J=3-2 line emission in the closest protoplanetary disk, TW Hya, a gas tracer that probes the bulk mass distribution. Combining it with the C$^{18}$O J=3-2 emission and the previously detected HD J=1-0 flux, we directly constrain the mid-plane temperature and optical depths of gas and dust emission. We report a gas mass distribution of 13$^{+8}_{-5}\times$(R/20.5AU)$^{-0.9^{+0.4}_{-0.3}}$ g cm$^{-2}$ in the expected formation zone of gas and ice giants (5-21AU). We find the total gas/millimeter-sized dust mass ratio is 140 in this region, suggesting that at least 2.4M_earth of dust aggregates have grown to >centimeter sizes (and perhaps much larger). The radial distribution of gas mass is consistent with a self-similar viscous disk profile but much flatter than the posterior extrapolation of mass distribution in our own and extrasolar planetary systems.

We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems.

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.

We model the radiatively driven flow from IM Lup -- a large protoplanetary disc expected to be irradiated by only a weak external radiation field (at least 10$^4$ times lower than the UV field irradiating the Orion Nebula Cluster proplyds). We find that material at large radii ($>400$AU) in this disc is sufficiently weakly gravitationally bound that significant mass loss can be induced. Given the estimated values of the disc mass and accretion rate, the viscous timescale is long ($\sim10$Myr) so the main evolutionary behaviour for the first Myr of the disc's lifetime is truncation of the disc by photoevaporation, with only modest changes effected by viscosity. We also produce approximate synthetic observations of our models, finding substantial emission from the flow which can explain the CO halo observed about IM Lup out to $\geq1000$AU. Solutions that are consistent with the extent of the observed CO emission generally imply that IM Lup is still in the process of having its disc outer radius truncated. We conclude that IM Lup is subject to substantial external photoevaporation, which raises the more general possibility that external irradiation of the largest discs can be of significant importance even in low mass star forming regions.

The spatial distribution of gas and solids in protoplanetary disks determines the composition and formation efficiency of planetary systems. A number of disks show starkly different distributions for the gas and small grains compared to millimeter-centimeter sized dust. We present new Atacama Large Millimeter/Submillimeter Array (ALMA) observations of the dust continuum, CO, $^{13}$CO, and C$^{18}$O in the IM Lup protoplanetary disk, one of the first systems where this dust-gas dichotomy was clearly seen. The $^{12}$CO is detected out to a radius of 970 AU, while the millimeter continuum emission is truncated at just 313 AU. Based upon this data, we have built a comprehensive physical and chemical model for the disk structure, which takes into account the complex, coupled nature of the gas and dust and the interplay between the local and external environment. We constrain the distributions of gas and dust, the gas temperatures, the CO abundances, the CO optical depths, and the incident external radiation field. We find that the reduction/removal of dust from the outer disk exposes this region to higher stellar and external radiation and decreases the rate of freeze-out, allowing CO to remain in the gas out to large radial distances. We estimate a gas-phase CO abundance of $5\%$ of the ISM value and a low external radiation field ($G_0\lesssim4$). The latter is consistent with that expected from the local stellar population. We additionally find tentative evidence for ring-like continuum substructure, suggestions of isotope-selective photodissociation, and a diffuse gas halo.

We report observations of resolved C2H emission rings within the gas-rich protoplanetary disks of TW Hya and DM Tau using the Atacama Large Millimeter Array (ALMA). In each case the emission ring is found to arise at the edge of the observable disk of mm-sized grains (pebbles) traced by (sub)mm-wave continuum emission. In addition, we detect a C3H2 emission ring with an identical spatial distribution to C2H in the TW Hya disk. This suggests that these are hydrocarbon rings (i.e. not limited to C2H). Using a detailed thermo-chemical model we show that reproducing the emission from C2H requires a strong UV field and C/O > 1 in the upper disk atmosphere and outer disk, beyond the edge of the pebble disk. This naturally arises in a disk where the ice-coated dust mass is spatially stratified due to the combined effects of coagulation, gravitational settling and drift. This stratification causes the disk surface and outer disk to have a greater permeability to UV photons. Furthermore the concentration of ices that transport key volatile carriers of oxygen and carbon in the midplane, along with photochemical erosion of CO, leads to an elemental C/O ratio that exceeds unity in the UV-dominated disk. Thus the motions of the grains, and not the gas, lead to a rich hydrocarbon chemistry in disk surface layers and in the outer disk midplane.

Nitrogen chemistry in protoplanetary disks and the freeze-out on dust particles is key to understand the formation of nitrogen bearing species in early solar system analogs. So far, ammonia has not been detected beyond the snowline in protoplanetary disks. We aim to find gas-phase ammonia in a protoplanetary disk and characterize its abundance with respect to water vapor. Using HIFI on the Herschel Space Observatory we detect, for the first time, the ground-state rotational emission of ortho-NH$_3$ in a protoplanetary disk, around TW Hya. We use detailed models of the disk's physical structure and the chemistry of ammonia and water to infer the amounts of gas-phase molecules of these species. We explore two radial distributions ( confined to $<$60 au like the millimeter-sized grains) and two vertical distributions (near the midplane where water is expected to photodesorb off icy grains) to describe the (unknown) location of the molecules. These distributions capture the effects of radial drift and vertical settling of ice-covered grains. We use physical-chemical models to reproduce the fluxes with assuming that water and ammonia are co-spatial. We infer ammonia gas-phase masses of 0.7-11.0 $\times$10$^{21}$ g. For water, we infer gas-phase masses of 0.2-16.0 $\times$10$^{22}$ g. This corresponds to NH$_3$/H$_2$O abundance ratios of 7\%-84\%, assuming that water and ammonia are co-located. Only in the most compact and settled adopted configuration is the inferred NH$_3$/H$_2$O consistent with interstellar ices and solar system bodies of $\sim$ 5\%-10\%. Volatile release in the midplane may occur via collisions between icy bodies if the available surface for subsequent freeze-out is significantly reduced, e.g., through growth of small grains into pebbles or larger.

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.

Deuterium-to-hydrogen (D/H) enrichments in molecular species provide clues about their original formation environment. The organic materials in primitive solar system bodies have generally higher D/H ratios and show greater D/H variation when compared to D/H in solar system water. We propose this difference arises at least in part due to 1) the availability of additional chemical fractionation pathways for organics beyond that for water, and 2) the higher volatility of key carbon reservoirs compared to oxygen. We test this hypothesis using detailed disk models, including a sophisticated, new disk ionization treatment with a low cosmic ray ionization rate, and find that disk chemistry leads to higher deuterium enrichment in organics compared to water, helped especially by fractionation via the precursors CH$_2$D$^+$/CH$_3^+$. We also find that the D/H ratio in individual species varies significantly depending on their particular formation pathways. For example, from $\sim20-40$ AU, CH$_4$ can reach $\rm{D/H\sim2\times10^{-3}}$, while D/H in CH$_3$OH remains locally unaltered. Finally, while the global organic D/H in our models can reproduce intermediately elevated D/H in the bulk hydrocarbon reservoir, our models are unable to reproduce the most deuterium-enriched organic materials in the solar system, and thus our model requires some inheritance from the cold interstellar medium from which the Sun formed.

An unsolved problem in step-wise core-accretion planet formation is that rapid radial drift in gas-rich protoplanetary disks should drive mm-/meter-sized particles inward to the central star before large bodies can form. One promising solution is to confine solids within small scale structures. Here we investigate dust structures in the (sub)mm continuum emission of four disks (TW Hya, HL Tau, HD 163296 and DM Tau), a sample of disks with the highest spatial resolution ALMA observations to date. We retrieve the surface brightness distributions using synthesized images and fitting visibilities with analytical functions. We find that the continuum emission of the four disks is ~axi-symmetric but rich in 10-30AU-sized radial structures, possibly due to physical gaps, surface density enhancements or localized dust opacity variations within the disks. These results suggest that small scale axi-symmetric dust structures are likely to be common, as a result of ubiquitous processes in disk evolution and planet formation. Compared with recent spatially resolved observations of CO snowlines in these same disks, all four systems show enhanced continuum emission from regions just beyond the CO condensation fronts, potentially suggesting a causal relationship between dust growth/trapping and snowlines.

Observations of protoplanetary disks suggest that the gas and dust follow significantly different radial distributions. This finding can be theoretically explained by a combination of radial drift and gas drag of intermediate-sized dust grains. Using a simple parametric model to approximate the different distributions of the gas and dust components, we calculate and examine the impact of radial drift on the global dust temperature structure. We find that the removal of large grains beyond the "truncation radius" allows this region to become significantly warmer from reprocessed stellar radiation shining down from the disk upper layers, increasing the outer disk temperature by $\sim10-30\%$. This change is sufficient to raise the local temperature to a value exceeding the CO desorption temperature. These findings imply that the disk density structures induced by radial drift are able to create multiple CO snow-lines. The inner disk CO is in the gas phase, freezing out near the classical snow-line at $R\sim20-40$ AU. Moving outward, the CO sublimates once again beyond the truncation radius (80 AU in our models) and subsequently re-freezes out at sufficiently large stellar distances, beyond $R\gtrsim130-200$ AU. We find that thermal desorption of CO in the outer disk becomes competitive with external UV photodesorption and that this additional transition from solid state CO to the gas-phase has significant implications for the C/O ratio in the outer disk.

H$_2$CO ice on dust grains is an important precursor of complex organic molecules (COMs). H$_2$CO gas can be readily observed in protoplanetary disks and may be used to trace COM chemistry. However, its utility as a COM probe is currently limited by a lack of constraints on the relative contributions of two different formation pathways: on icy grain-surfaces and in the gas-phase. We use archival ALMA observations of the resolved distribution of H$_2$CO emission in the disk around the young low-mass star DM Tau to assess the relative importance of these formation routes. The observed H$_2$CO emission has a centrally peaked and radially broad brightness profile (extending out to 500 AU). We compare these observations with disk chemistry models with and without grain-surface formation reactions, and find that both gas and grain-surface chemistry are necessary to explain the spatial distribution of the emission. Gas-phase H$_2$CO production is responsible for the observed central peak, while grain-surface chemistry is required to reproduce the emission exterior to the CO snowline (where H$_2$CO mainly forms through the hydrogenation of CO ice before being non-thermally desorbed). These observations demonstrate that both gas and grain-surface pathways contribute to the observed H$_2$CO in disks, and that their relative contributions depend strongly on distance from the host star.

We examine changes in the molecular abundances resulting from increased heating due to a self-luminous planetary companion embedded within a narrow circumstellar disk gap. Using 3D models that include stellar and planetary irradiation, we find that luminous young planets locally heat up the parent circumstellar disk by many tens of Kelvin, resulting in efficient thermal desorption of molecular species that are otherwise locally frozen out. Furthermore, the heating is deposited over large regions of the disk, $\pm5$ AU radially and spanning $\lesssim60^\circ$ azimuthally. From the 3D chemical models, we compute rotational line emission models and full ALMA simulations, and find that the chemical signatures of the young planet are detectable as chemical asymmetries in $\sim10h$ observations. HCN and its isotopologues are particularly clear tracers of planetary heating for the models considered here, and emission from multiple transitions of the same species is detectable, which encodes temperature information in addition to possible velocity information from the spectra itself. We find submillimeter molecular emission will be a useful tool to study gas giant planet formation in situ, especially beyond $R\gtrsim10$ AU.

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.

The recent (apparent) passage of the Voyager 1 spacecraft into interstellar space provides us with front-row seats to the complex interplay between the solar wind and the protective surrounding bubble known as heliosphere. The heliosphere extends radially out to $\sim100$ AU from the sun, and within this sphere of influence, the solar wind modulates the incoming flux of galactic cosmic rays (CRs), especially those at low energies. Newly formed stars, which support both strong magnetic fields and winds, are expected to produce analogous regions of CR exclusion, perhaps at elevated levels. Such young stars are encircled by molecular gas-rich disks, and the net removal of CRs from the circumstellar environment significantly reduces the expected CR ionization rate in the disk gas, most likely by many orders-of-magnitude. The loss of ionization reduces disk turbulence, and thereby affects both planet-formation and active chemical processes in the disk. We present models of CR exclusion and explore the implications for turbulence and for predicted chemical abundances. We also discuss means by which ALMA can be used to search for extrasolar heliosphere-analogs around young stars.

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 Herschel Space Observatory PACS spectra of T Tauri stars, in which we detect amorphous and crystalline water ice features. Using irradiated accretion disk models, we determine the disk structure and ice abundance in each of the systems. Combining a model-independent comparison of the ice feature strength and disk size with a detailed analysis of the model ice location, we estimate that the ice emitting region is at disk radii >30AU, consistent with a proto-Kuiper belt. Vertically, the ice emits most below the photodesorption zone, consistent with Herschel observations of cold water vapor. The presence of crystallized water ice at a disk location a) colder than its crystallization temperature and b) where it should have been re-amorphized in ~1 Myr suggests that localized generation is occurring; the most likely cause appears to be micrometeorite impact or planetesimal collisions. Based on simple tests with UV models and different ice distributions, we suggest that the SED shape from 20 to 50 micron may probe the location of the water ice snow line in the disk upper layers. This project represents one of the first extra-solar probes of the spatial structure of the cometary ice reservoir thought to deliver water to terrestrial planets.

Identifying the source of Earth's water is central to understanding the origins of life-fostering environments and to assessing the prevalence of such environments in space. Water throughout the solar system exhibits deuterium-to-hydrogen enrichments, a fossil relic of low-temperature, ion-derived chemistry within either (i) the parent molecular cloud or (ii) the solar nebula protoplanetary disk. Utilizing a comprehensive treatment of disk ionization, we find that ion-driven deuterium pathways are inefficient, curtailing the disk's deuterated water formation and its viability as the sole source for the solar system's water. This finding implies that if the solar system's formation was typical, abundant interstellar ices are available to all nascent planetary systems.

The chemical properties of protoplanetary disks are especially sensitive to their ionization environment. Sources of molecular gas ionization include cosmic rays, stellar X-rays and short-lived radionuclides, each of which varies with location in the disk. This behavior leads to a significant amount of chemical structure, especially in molecular ion abundances, which is imprinted in their submillimeter rotational line emission. Using an observationally motivated disk model, we make predictions for the dependence of chemical abundances on the assumed properties of the ionizing field. We calculate the emergent line intensity for abundant molecular ions and simulate sensitive observations with the Atacama Large Millimeter/Sub-millimeter Array (ALMA) for a disk at D=100 pc. The models readily distinguish between high ionization rates ($\zeta\gtrsim10^{-17}$ s$^{-1}$ per H$_2$) and below, but it becomes difficult to distinguish between low ionization models when $\zeta\lesssim10^{-19}$ s$^{-1}$. We find that \htdp\ emission is not detectable for sub-interstellar CR rates with ALMA (6h integration), and that \ntdp\ emission may be a more sensitive tracer of midplane ionization. HCO$^+$ traces X-rays and high CR rates ($\zeta_{\rm{CR}}\gtrsim10^{-17}$ s$^{-1}$), and provides a handle on the warm molecular ionization properties where CO is present in the gas. Furthermore, species like HCO$^+$, which emits from a wide radial region and samples a large gradient in temperature, can exhibit ring-like emission as a consequence of low-lying rotational level de-excitation near the star. This finding highlights a scenario where rings are not necessarily structural or chemical in nature, but simply a result of the underlying line excitation properties.

Given the central role of carbon in the chemistry of life, it is a fundamental question as to how carbon is supplied to the Earth, in what form and when. We provide an accounting of carbon found in solar system bodies, in particular a comparison between the organic content of meteorites and that in identified organics in the dense interstellar medium (ISM). Based on this accounting identified organics created by the chemistry of star formation could contain at most ~15% of the organic carbon content in primitive meteorites and significantly less for cometary organics, which represent the putative contributors to starting materials for the Earth. In the ISM ~30% of the elemental carbon is found in CO, either in the gas or ices, with a typical abundance of ~10^-4 (relative to H2). Recent observations of the TW Hya disk find that the gas phase abundance of CO is reduced by an order of magnitude compared to this value. We explore a solution where the volatile CO is destroyed via a gas phase processes, providing an additional source of carbon for organic material to be incorporated into planetesimals and cometesimals. This chemical processing mechanism requires warm grains (> 20 K), partially ionized gas, and sufficiently small <10 micron grains, i.e. a larger total grain surface area, such that freeze-out is efficient. Under these conditions static (non-turbulent) chemical models predict that a large fraction of the carbon nominally sequestered in CO can be the source of carbon for a wide variety of organics that are present as ice coatings on the surfaces of warm pre-planetesimal dust grains.

We present simple analytic solutions for the ionization rate $\zeta_{\rm{SLR}}$ arising from the decay of short-lived radionuclides (SLRs) within protoplanetary disks. We solve the radiative transfer problem for the decay products within the disk, and thereby allow for the loss of radiation at low disk surface densities; energy loss becomes important outside $R\gtrsim30$ for typical disk masses $M_g=0.04$ M$_\odot$. Previous studies of chemistry/physics in these disks have neglected the impact of ionization by SLRs, and often consider only cosmic rays (CRs), because of the high CR-rate present in the ISM. However, recent work suggests that the flux of CRs present in the circumstellar environment could be substantially reduced by relatively modest stellar winds, resulting in severely modulated CR ionization rates, $\zeta_{\rm{CR}}$, equal to or substantially below that of SLRs ($\zeta_{\rm{SLR}}\lesssim10^{-18}$ s$^{-1}$). We compute the net ionizing particle fluxes and corresponding ionization rates as a function of position within the disk for a variety of disk models. The resulting expressions are especially simple for the case of vertically gaussian disks (frequently assumed in the literature). Finally, we provide a power-law fit to the ionization rate in the midplane as a function of gas disk surface density and time. Depending on location in the disk, the ionization rates by SLRs are typically in the range $\zeta_{\rm{SLR}}\sim(1-10)\times10^{-19}$ s$^{-1}$.

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.

(Abridged) Cosmic rays (CRs) are thought to provide an important source of ionization in the outermost and densest regions of protoplanetary disks; however, it is unknown to what degree they are physically present. As is observed in the Solar System, stellar winds can inhibit the propagation of cosmic rays within the circumstellar environment and subsequently into the disk. In this work, we explore the hitherto neglected effects of cosmic ray modulation by both stellar winds and magnetic field structures and study how these processes act to reduce disk ionization rates. We construct a two-dimensional protoplanetary disk model of a T-Tauri star system, focusing on ionization from stellar and interstellar FUV, stellar X-ray photons, and cosmic rays. We show that stellar winds can power a Heliosphere-like analogue, i.e., a "T-Tauriosphere," diminishing cosmic ray ionization rates by several orders of magnitude at low to moderate CR energies (E_CR<1 GeV). We explore models of both the observed solar wind cosmic ray modulation and a highly simplified estimate for "elevated" cosmic ray modulation as would be expected from a young T-Tauri star. In the former (solar analogue) case, we estimate the ionization rate from galactic cosmic rays to be zeta_CR~(0.23-1.4) x 1e-18 s^-1. This range of values, which we consider to be the maximum CR ionization rate for the disk, is more than an order of magnitude lower than what is generally assumed in current models for disk chemistry and physics. In the latter case, i.e., for a "T-Tauriosphere," the ionization rate by cosmic rays is zeta_CR<1e-20 s^-1, which is 1000 times smaller than the interstellar value. Indeed, if winds are as efficient at cosmic ray modulation as predicted here, short-lived radionuclides (now extinct) would have provided the major source of ionization (zeta_RN~7.3e-19 s^-1) in the planet-forming zone of the young Solar Nebula.

Regions of massive star formation offer some of the best and most easily-observed examples of radiation hydrodynamics. Boundaries where fully-ionized H II regions transition to neutral/molecular photodissociation regions (PDRs) are of particular interest because marked temperature and density contrasts across the boundaries lead to evaporative flows and fluid dynamical instabilities that can evolve into spectacular pillar-like structures. When detached from their parent clouds, pillars become ionized globules that often harbor one or more young stars. H2 molecules at the interface between a PDR and an H II region absorb ultraviolet light from massive stars, and the resulting fluoresced infrared emission lines are an ideal way to trace this boundary independent of obscuring dust. This paper presents H2 images of four regions of massive star formation that illustrate different types of PDR boundaries. The Ara OB1 star formation region contains a striking long wall that has several wavy structures which are present in H2, but the emission is not particularly bright because the ambient UV fluxes are relatively low. In contrast, the Carina star formation region shows strong H2 fluorescence both along curved walls and at the edges of spectacular pillars that in some cases have become detached from their parent clouds. The less-spectacular but more well-known Eagle Nebula has two regions that have strong fluorescence in addition to its pillars. While somewhat older than the other regions, Cyg OB2 has the highest number of massive stars of the regions surveyed and contains many isolated, fluoresced globules that have head-tail morphologies which point towards the sources of ionizing radiation. These images provide a collection of potential astrophysical analogs that may relate to ablated interfaces observed in laser experiments of radiation hydrodynamics.

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 explore the chemical structure of a disk that contains a large central gap of R ~ 45 AU, as is commonly seen in transitional disk systems. In our chemical model of a disk with a cleared inner void, the midplane becomes revealed to the central star so that it is directly irradiated. The midplane material at the truncation radius is permissive to reprocessed optical heating radiation, but opaque to the photo-dissociating ultraviolet, creating an environment abundant in gas-phase molecules. Thus the disk midplane, which would otherwise for a full disk be dominated by near complete heavy element freeze-out, should become observable in molecular emission. If this prediction is correct this has exciting prospects for observations with the Atacama Large Millimeter/Submillimeter Array (ALMA), as the inner transition region should thus be readily detected and resolved, especially using high-J rotational transitions excited in the high density midplane gas. Therefore such observations will potentially provide us with a direct probe of the physics and chemistry at this actively evolving interface.

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extra-solar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on-board the Herschel Space Observatory has detected emission from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface hinting at a water ice reservoir equivalent to several thousand Earth Oceans in mass. The water's ortho-to-para ratio falls well below that of Solar System comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.