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.

Complex organic molecules (COMs) can be produced by energetic processing of interstellar ice mantles accreted on top of dust grains. Two COMs with proposed energetic ice formation pathways are formamide and acetaldehyde. Both have been detected in Solar System comets, and in different circumstellar and interstellar environments. In this work, we study the NH2CHO and CH3CHO formation upon UV photoprocessing of CO:NH3 and CO:CH4 ice samples. The conversion from NH2 radicals to NH2CHO is 2-16 times higher than the conversion from CH3 radicals to CH3CHO under the explored experimental conditions, likely because the formation of the latter competes with the formation of larger hydrocarbons. In addition, the conversion of NH2 into NH2CHO at 10 K increases with the NH3 abundance in the ice, and also with the temperature in CO-dominated CO:NH3 ices. This is consistent with the presence of a small NH2 and HCO reorientation barrier for the formation of NH2CHO, which is overcome with an increase in the ice temperature. The measured NH2CHO and CH3CHO formation efficiencies and rates are similar to those found during electron irradiation of the same ice samples under comparable conditions, suggesting that both UV photons and cosmic rays would have similar contributions to the solid-state formation of these species in space. Finally, the measured conversion yields (up to one order of magnitude higher for NH2CHO) suggest that in circumstellar environments, where the observed NH2CHO/CH3CHO abundance ratio is 0.1, there are likely additional ice and/or gas phase formation pathways for CH3CHO.

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

Phosphorus is a key ingredient in terrestrial biochemistry, but is rarely observed in the molecular ISM and therefore little is known about how it is inherited during the star and planet formation sequence. We present observations of the phosphorus-bearing molecules PO and PN towards the Class I low-mass protostar B1-a using the IRAM 30m telescope, representing the second detection of phosphorus carriers in a Solar-type star forming region. The P/H abundance contained in PO and PN is ~10$^{-10}$-10$^{-9}$ depending on the assumed source size, accounting for just 0.05-0.5% of the solar phosphorus abundance and implying significant sequestration of phosphorus in refractory material. Based on a comparison of the PO and PN line profiles with the shock tracers SiO, SO$_2$, and CH$_3$OH, the phosphorus molecule emission seems to originate from shocked gas and is likely associated with a protostellar outflow. We find a PO/PN column density ratio of ~1-3, which is consistent with the values measured in the shocked outflow of the low-mass protostar L1157, the massive star-forming regions W51 and W3(OH), and the galactic center GMC G+0.693-0.03. This narrow range of PO/PN ratios across sources with a range of environmental conditions is surprising, and likely encodes information on how phosphorus carriers are stored in grain mantles.

Jupiter's atmosphere is enriched in C, N, S, P, Ar, Kr and Xe with respect to solar abundances by a factor of ~3. Gas Giant envelopes are mainly enriched through the dissolution of solids in the atmosphere, and this constant enrichment factor is puzzling since several of the above elements are not expected to have been in the solid phase in Jupiter's feeding zone; most seriously, Ar and the main carrier of N, N2, only condense at the very low temperatures, 21-26 K, associated with the outer solar nebula. We propose that a plausible solution to the enigma of Jupiter's uniform enrichment pattern is that Jupiter's core formed exterior to the N2 and Ar snowlines, beyond 30 au, resulting in a Solar composition core in all volatiles heavier than Ne. During envelope accretion and planetesimal bombardment, some of the core mixed in with the envelope causing the observed enrichment pattern. We show that this scenario naturally produces the observed atmosphere composition, even with substantial pollution from N-poor pebble and planetesimal accretion in Jupiter's final feeding zone. We note that giant core formation at large nebular radii is consistent with recent models of gas giant core formation through pebble accretion, which requires the core to form exterior to Jupiter's current location to counter rapid inward migration during the core and envelope formation process. If this scenario is common, gas giant core formation may account for many of the gaps observed in protoplanetary disks between 10s and 100 au.

Planet atmosphere and hydrosphere compositions are fundamentally set by accretion of volatiles, and therefore by the division of volatiles between gas and solids in planet-forming disks. For hyper-volatiles such as CO, this division is regulated by a combination of binding energies, and by the ability of other ice components to entrap. Water ice is known for its ability to trap CO and other volatile species. In this study we explore whether another common interstellar and cometary ice component, CO2, is able to trap CO as well. We measure entrapment of CO molecules in CO2 ice through temperature programmed desorption (TPD) experiments on CO2:CO ice mixtures. We find that CO2 ice traps CO with a typical efficiency of 40-60% of the initially deposited CO molecules for a range of ice thicknesses between 7 and 50ML, and ice mixture ratios between 1:1 and 9:1. The entrapment efficiency increases with ice thickness and CO dilution. We also run analogous H2O:CO experiments and find that under comparable experimental conditions CO2 ice entraps CO more efficiently than H2O ice up to the onset of CO2 desorption at ~70K. We speculate that this may be due to different ice restructuring dynamics in H2O and CO2 ices around the CO desorption temperature. Importantly, the ability of CO2 to entrap CO may change the expected division between gas and solids for CO and other hyper-volatiles exterior to the CO2 snowline during planet formation.

We present ALMA observations of organic molecules towards five low-mass Class 0/I protostellar disk candidates in the Serpens cluster. Three sources (Ser-emb 1, Ser-emb 8, and Ser-emb 17) present emission of CH3OH as well as CH3OCH3, CH3OCHO, and CH2CO, while NH2CHO is detected in just Ser-emb 8 and Ser-emb 17. Detecting hot corino-type chemistry in three of five sources represents a high occurrence rate given the relative sparsity of these sources in the literature, and this suggests a possible link between protostellar disk formation and hot corino formation. For sources with CH3OH detections, we derive column densities of 10^{17}-10^{18} cm^{-2} and rotational temperatures of ~200-250 K. The CH3OH-normalized column density ratios of large, oxygen-bearing COMs in the Serpens sources and other hot corinos span two orders of magnitude, demonstrating a high degree of chemical diversity at the hot corino stage. Resolved observations of a larger sample of objects are needed to understand the origins of chemical diversity in hot corinos, and the relationship between different protostellar structural elements on disk-forming scales.

We have observed 29 transitions corresponding to 12 distinct species and 7 additional isotopologs toward the deeply embedded Class 0 young stellar object Ser-emb 1 with the Atacama Large Millimeter Array at 1 mm. The detected species include CH3OH and two complex organic molecules, CH3OCH3 and CH3OCHO. The emission of CH3OH and the two COMs is compact, and the CH3OH rotational temperature is 261 (46) K, implying that Ser-emb 1 hosts a hot corino. The derived CH3OH, CH3OCH3 and CH3OCHO column densities are at least 1.2 (0.4)E17 cm-2, 9.2 (3.8)E16 cm-2, and 9.1 (3.6)E16 cm-2, respectively, comparable to the values found for other Class 0 hot corinos. In addition, we observe evidence of rotation at compact scales: two of the more strongly detected lines, corresponding to C18O and H2CO present spatially resolved red- and blue-shifted compact emission orthogonal to the direction of a jet and outflow traced by CO, SiO, and several other molecules. The spatial coincidence of the hot corino emission and a possible disk in a compact region around the central protostar suggests that these structures may be physically and/or chemically related.

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.

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.

Oxygen atom addition and insertion reactions may provide a pathway to chemical complexity in ices that are too cold for radicals to diffuse and react. We have studied the ice-phase reactions of photo-produced oxygen atoms with C2 hydrocarbons under ISM-like conditions. The main products of oxygen atom reactions with ethane are ethanol and acetaldehyde; with ethylene are ethylene oxide and acetaldehyde; and with acetylene is ketene. The derived branching ratio from ethane to ethanol is ~0.74 and from ethylene to ethylene oxide is ~0.47. For all three hydrocarbons there is evidence of an effectively barrierless reaction with O(^1D) to form oxygen-bearing organic products; in the case of ethylene, there may be an additional barriered contribution of the ground-state O(^3P) atom. Thus, oxygen atom reactions with saturated and unsaturated hydrocarbons are a promising pathway to chemical complexity even at very low temperatures where the diffusion of radical species is thermally inaccessible.

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.

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.

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.

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.

Carbon chain molecules may be an important reservoir of reactive organics during star and planet formation. Carbon chains have been observed toward several low-mass young stellar objects (YSOs), but their typical abundances and chemical relationships in such sources are largely unconstrained. We present a carbon chain survey toward 16 deeply embedded (Class 0/I) low-mass protostars made with the IRAM 30 m telescope. Carbon chains are found to be common at this stage of protostellar evolution. We detect CCS, CCCS, HC$_3$N, HC$_5$N, l-C$_3$H, and C$_4$H toward 88%, 38%, 75%, 31%, 81%, and 88% of sources, respectively. Derived column densities for each molecule vary by one to two orders of magnitude across the sample. As derived from survival analysis, median column densities range between 1.2$\times 10^{11}$ cm$^{-2}$ (CCCS) and 1.5$\times 10^{13}$ cm$^{-2}$ (C$_4$H) and estimated fractional abundances with respect to hydrogen range between 2$\times 10^{-13}$ (CCCS) and 5$\times 10^{-11}$ (C$_4$H), which are low compared to cold cloud cores, warm carbon chain chemistry (WCCC) sources, and protostellar model predictions. We find significant correlations between molecules of the same carbon chain families, as well as between the cyanpolyynes (HC$_{\rm n}$N) and the pure hydrocarbon chains (C$_{\rm n}$H). This latter correlation is explained by a closely-related production chemistry of C$_{\rm{n}}$H and cyanpolyynes during low-mass star formation.

The organic content of protoplanetary disks sets the initial compositions of planets and comets, thereby influencing subsequent chemistry that is possible in nascent planetary systems. We present observations of the complex nitrile-bearing species CH3CN and HC3N towards the disks around the T Tauri stars AS 209, IM Lup, LkCa 15, and V4046 Sgr as well as the Herbig Ae stars MWC 480 and HD 163296. HC3N is detected towards all disks except IM Lup, and CH3CN is detected towards V4046 Sgr, MWC 480, and HD 163296. Rotational temperatures derived for disks with multiple detected lines range from 29-73K, indicating emission from the temperate molecular layer of the disk. V4046 Sgr and MWC 480 radial abundance profiles are constrained using a parametric model; the gas-phase CH3CN and HC3N abundances with respect to HCN are a few to tens of percent in the inner 100 AU of the disk, signifying a rich nitrile chemistry at planet- and comet-forming disk radii. We find consistent relative abundances of CH3CN, HC3N, and HCN between our disk sample, protostellar envelopes, and solar system comets; this is suggestive of a robust nitrile chemistry with similar outcomes under a wide range of physical conditions.

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

We 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 present experimental constraints on the insertion of oxygen atoms into methane to form methanol in astrophysical ice analogs. In gas-phase and theoretical studies this process has previously been demonstrated to have a very low or non-existent energy barrier, but the energetics and mechanisms have not yet been characterized in the solid state. We use a deuterium UV lamp filtered by a sapphire window to selectively dissociate O2 within a mixture of O2:CH4 and observe efficient production of CH3OH via O(1D) insertion. CH3OH growth curves are fit with a kinetic model, and we observe no temperature dependence of the reaction rate constant at temperatures below the oxygen desorption temperature of 25K. Through an analysis of side products we determine the branching ratio of ice-phase oxygen insertion into CH4: ~65% of insertions lead to CH3OH with the remainder leading instead to H2CO formation. There is no evidence for CH3 or OH radical formation, indicating that the fragmentation is not an important channel and that insertions typically lead to increased chemical complexity. CH3OH formation from O2 and CH4 diluted in a CO-dominated ice similarly shows no temperature dependence, consistent with expectations that insertion proceeds with a small or non-existent barrier. Oxygen insertion chemistry in ices should therefore be efficient under low-temperature ISM-like conditions, and could provide an important channel to complex organic molecule formation on grain surfaces in cold interstellar regions such as cloud cores and protoplanetary disk midplanes.

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.

Volatiles, especially CO, are important gas tracers of protoplanetary disks (PPDs). Freeze-out and sublimation processes determine their division between gas and solid phases, which affects both which disk regions can be traced by which volatiles, and the formation and composition of planets. Recently, multiple lines of evidence suggest that CO is substantially depleted from the gas in the outer regions of PPDs. In this letter, we show that the gas dynamics in the outer PPDs facilitates volatile depletion through a mechanism which we term "runaway freeze-out". Using a simple 1D model that incorporates dust settling, turbulent diffusion of dust and volatiles, as well as volatile freeze-out/sublimation processes, we show that as long as turbulence in the cold midplane is sufficiently weak to allow majority of the small grains to settle, CO in the warm surface layer can be turbulently mixed into the midplane region and depleted by freeze-out. The level of depletion sensitively depends on the level of disk turbulence. Based on recent disk simulations that suggest a layered turbulence profile with very weak midplane turbulence and strong turbulence at disk surface, CO and other volatiles can be efficiently depleted by a factor of a few over Myr timescales.

We aim to develop the Orion B Giant Molecular Cloud (GMC) as a local template for interpreting extra-galactic molecular line observations. We use the wide-band receiver at the IRAM-30m to spatially and spectrally resolve the Orion B GMC. The observations cover almost 1 square degree at 26" resolution with a bandwidth of 32 GHz from 84 to 116 GHz in only two tunings. Among the mapped spectral lines are the 12CO, 13CO, C18O, C17O, HCN, HNC, 12CN, CCH, HCO+, N2H+ (1-0), and 12CS, 32SO, SiO, c-C3H2, CH3OH (2-1) transitions. We introduce the molecular anatomy of the Orion B GMC, including relations between line intensities and gas column density or far-UV radiation fields, and correlations between selected line and line ratios. We also obtain a dust-traced gas mass that is less than about one third the CO-traced mass, using the standard Xco conversion factor. The presence of overluminous CO can be traced back to the dependence of the CO intensity on UV illumination. In fact, while most lines show some dependence on the UV radiation field, CN and CCH are the most sensitive. Moreover dense cloud cores are almost exclusively traced by N2H+. Other traditional high density tracers, such as HCN (1-0), are also easily detected in extended translucent regions at a typical density of about 500 H2 cm-3. In general, we find no straightforward relation between line critical density and the fraction of the line luminosity coming from dense gas regions. Our initial findings demonstrate that the relations between line (ratio) intensities and environment in GMCs are more complicated than often assumed. Sensitivity (i.e., the molecular column density), excitation, and above all chemistry contribute to the observed line intensity distributions. They must be considered together when developing the next generation of extra-galactic molecular line diagnostics of mass, density, temperature and radiation field.

The elemental compositions of planets define their chemistry, and could potentially be used as beacons for their formation location if the elemental gas and grain ratios of planet birth environments, i.e. protoplanetary disks, are well understood. In disks, the ratios of volatile elements, such as C/O and N/O, are regulated by the abundance of the main C, N, O carriers, their ice binding environment, and the presence of snowlines of major volatiles at different distances from the central star. We explore the effects of disk dynamical processes, molecular compositions and abundances, and ice compositions on the snowline locations of the main C, O and N carriers, and the C/N/O ratios in gas and dust throughout the disk. The gas-phase N/O ratio enhancement in the outer disk (exterior to the H2O snowline) exceeds the C/O ratio enhancement for all reasonable volatile compositions. Ice compositions and disk dynamics individually change the snowline location of N2, the main nitrogen carrier, by a factor of 2-3, and when considered together the range of possible N2 snowline locations is ~11- ~79 AU in a standard disk model. Observations that anchor snowline locations at different stages of planet formation are therefore key to develop C/N/O ratios as a probe of planet formation zones.

The compositions of nascent planets depend on the compositions of their birth disks. In particular, the elemental compositions of Gas Giant gaseous envelopes depend on the elemental composition of the disk gas from which the envelope is accreted. Previous models demonstrated that sequential freeze-out of O and C-bearing volatiles in disks will result in an supersolar C/O ratios and subsolar C/H ratios in the gas between water and CO snowlines. This result does not take into account, however, the expected grain growth and radial drift of pebbles in disks, and the accompanying re-distribution of volatiles from the outer to the inner disk. Using a toy model we demonstrate that when drift is considered, CO is enhanced between the water and CO snowline, resulting in both supersolar C/O and C/H ratios in the disk gas in the Gas Giant formation zone. This result appears robust to the details of the disk model as long as there is substantial pebble drift across the CO snowline, and the efficiency of CO vapor diffusion is limited. Gas Giants that accrete their gaseous envelopes exterior to the water snowline and do not experience substantial core-envelope mixing, may thus present both superstellar C/O and C/H ratios in their atmospheres. Pebble drift will also affect the nitrogen and noble gas abundances in the planet forming zones, which may explain some of Jupiter's peculiar abundance patterns.

The chemistry of dense interstellar regions was analyzed using a time-dependent gas-grain astrochemical simulation and a new chemical network that incorporates deuterated chemistry taking into account nuclear spin-states for the hydrogen chemistry and its deuterated isotopologues. With this new network, the utility of the [HCO$^+$]/[DCO$^+$] abundance ratio as a probe of the cosmic ray ionization rate has been reexamined, with special attention paid to the effect of the initial value of the molecular hydrogen ortho-to-para ratio (OPR). After discussing the use of the probe for cold cores, we then compare our results with previous theoretical and observational results for a molecular cloud close to the supernova remnant W51C, which is thought to have an enhanced cosmic ray ionization rate $\zeta$ caused by the nearby $\gamma$-ray source. In addition, we attempt to use our approach to estimate the cosmic ray ionization rate for L1174, a dense core with an embedded star. Beyond the previously known sensitivity of [HCO$^+$]/[DCO$^+$] to $\zeta$, we demonstrate its additional dependence on the initial OPR and, secondarily, on the age of the source, its temperature, and its density. We conclude that the usefulness of the [HCO$^+$]/[DCO$^+$] abundance ratio to constrain the cosmic ray ionization rate in dense regions increases with source age and ionization rate as the ratio becomes far less sensitive to the initial value of the OPR.

The interstellar medium is characterized by a rich and diverse chemistry. Many of its complex organic molecules are proposed to form through radical chemistry in icy grain mantles. Radicals form readily when interstellar ices (composed of water and other volatiles) are exposed to UV photons and other sources of dissociative radiation, and, if sufficiently mobile, the radicals can react to form larger, more complex molecules. The resulting complex organic molecules (COMs) accompany star and planet formation, and may eventually seed the origins of life on nascent planets. Experiments of increasing sophistication have demonstrated that known interstellar COMs as well as the prebiotically interesting amino acids can form through ice photochemistry. We review these experiments and discuss the qualitative and quantitative kinetic and mechanistic constraints they have provided. We finally compare the effects of UV radiation with those of three other potential sources of radical production and chemistry in interstellar ices: electrons, ions and X-rays.

Interstellar complex organic molecules (COMs) are commonly observed during star formation, and are proposed to form through radical chemistry in icy grain mantles. Reactions between ions and neutral molecules in ices may provide an alternative cold channel to complexity, as ion-neutral reactions are thought to have low or even no energy barriers. Here we present a study of a the kinetics and mechanisms of a potential ion-generating acid-base reaction between NH$_{3}$ and HCOOH to form the salt NH$_{4}^{+}$HCOO$^{-}$. We observe salt growth at temperatures as low as 15K, indicating that this reaction is feasible in cold environments. The kinetics of salt growth are best fit by a two-step model involving a slow "pre-reaction" step followed by a fast reaction step. The reaction energy barrier is determined to be 70 $\pm$ 30K with a pre-exponential factor 1.4 $\pm$ 0.4 x 10$^{-3}$ s$^{-1}$. The pre-reaction rate varies under different experimental conditions and likely represents a combination of diffusion and orientation of reactant molecules. For a diffusion-limited case, the pre-reaction barrier is 770 $\pm$ 110K with a pre-exponential factor of $\sim$7.6 x 10$^{-3}$ s$^{-1}$. Acid-base chemistry of common ice constituents is thus a potential cold pathway to generating ions in interstellar ices.

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 \sigma$), reaching a peak signal-to-noise of $5.5\sigma$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation.

Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=2-1$ emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the C$^{18}$O emission consists of a central peak and a ring at a radius of $\sim1''$ (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.

We present long-baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870 micron continuum emission from the nearest gas-rich protoplanetary disk, around TW Hya, that trace millimeter-sized particles down to spatial scales as small as 1 AU (20 mas). These data reveal a series of concentric ring-shaped substructures in the form of bright zones and narrow dark annuli (1-6 AU) with modest contrasts (5-30%). We associate these features with concentrations of solids that have had their inward radial drift slowed or stopped, presumably at local gas pressure maxima. No significant non-axisymmetric structures are detected. Some of the observed features occur near temperatures that may be associated with the condensation fronts of major volatile species, but the relatively small brightness contrasts may also be a consequence of magnetized disk evolution (the so-called zonal flows). Other features, particularly a narrow dark annulus located only 1 AU from the star, could indicate interactions between the disk and young planets. These data signal that ordered substructures on ~AU scales can be common, fundamental factors in disk evolution, and that high resolution microwave imaging can help characterize them during the epoch of planet formation.

Large interstellar organic molecules are potential precursors of prebiotic molecules. Their formation pathways and chemical relationships with one another and simpler molecules are therefore of great interest. In this paper, we address the relationship between two classes of large organic molecules, carbon chains and saturated complex organic molecules (COMs), at the early stages of star formation through observations of C4H and CH3OH. We surveyed these molecules with the IRAM 30m telescope toward 16 deeply embedded low-mass protostars selected from the Spitzer c2d ice survey. We find that CH3OH and C4H are positively correlated indicating that these two classes of molecules can coexist during the embedded protostellar stage. The C4H/CH3OH gas abundance ratio tentatively correlates with the CH4/CH3OH ice abundance ratio in the same lines of sight. This relationship supports a scenario where carbon chain formation in protostellar envelopes begins with CH4 ice desorption.

The C/O ratio is a defining feature of both gas giant atmospheric and protoplanetary disk chemistry. In disks, the C/O ratio is regulated by the presence of snowlines of major volatiles at different distances from the central star. We explore the effect of radial drift of solids and viscous gas accretion onto the central star on the snowline locations of the main C and O carriers in a protoplanetary disk, H2O, CO2 and CO, and their consequences for the C/O ratio in gas and dust throughout the disk. We determine the snowline locations for a range of fixed initial particle sizes and disk types. For our fiducial disk model, we find that grains with sizes ~0.5 cm < s < 7 m for an irradiated disk, and ~0.001 cm < s < 7 m for an evolving and viscous disk, desorb at a size-dependent location in the disk, which is independent of the particle's initial position. The snowline radius decreases for larger particles, up to sizes of ~7 m. Compared to a static disk, we find that radial drift and gas accretion in a viscous disk move the H2O snowline inwards by up to 40%, the CO2 snowline by up to 60%, and the CO snowline by up to 50%. We thus determine an inner limit on the snowline locations when radial drift and gas accretion are accounted for.

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.

Observations of deuterium fractionation in the solar system, and in interstellar and circumstellar material, are commonly used to constrain the formation environment of volatiles. Toward protoplanetary disks, this approach has been limited by the small number of detected deuterated molecules, i.e. DCO$^+$ and DCN. Based on ALMA Cycle 2 observations toward the disk around the T Tauri star AS 209, we report the first detection of N$_2$D$^+$ (J=3-2) in a protoplanetary disk. These data are used together with previous Submillimeter Array observations of N$_2$H$^+$ (J=3-2) to estimate a disk-averaged D/H ratio of 0.3--0.5, an order of magnitude higher than disk-averaged ratios previously derived for DCN/HCN and DCO$^+$/HCO$^+$ around other young stars. The high fractionation in N$_2$H$^+$ is consistent with model predictions. The presence of abundant N$_2$D$^+$ toward AS 209 also suggests that N$_2$D$^+$ and the N$_2$D$^+$/N$_2$H$^+$ ratio can be developed into effective probes of deuterium chemistry, kinematics, and ionization processes outside the CO snowline of disks.

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.

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.

The mobility of atoms, molecules and radicals in icy grain mantles regulate ice restructuring, desorption, and chemistry in astrophysical environments. Interstellar ices are dominated by H2O, and diffusion on external and internal (pore) surfaces of H2O-rich ices is therefore a key process to constrain. This study aims to quantify the diffusion kinetics and barrier of the abundant ice constituent CO into H2O dominated ices at low temperatures (15-23 K), by measuring the mixing rate of initially layered H2O(:CO2)/CO ices. The mixed fraction of CO as a function of time is determined by monitoring the shape of the infrared CO stretching band. Mixing is observed at all investigated temperatures on minute time scales, and can be ascribed to CO diffusion in H2O ice pores. The diffusion coefficient and final mixed fraction depend on ice temperature, porosity, thickness and composition. The experiments are analyzed by applying Fick's diffusion equation under the assumption that mixing is due to CO diffusion into an immobile H2O ice. The extracted energy barrier for CO diffusion into amorphous H2O ice is ~160 K. This is effectively a surface diffusion barrier. The derived barrier is low compared to current surface diffusion barriers in use in astrochemical models. Its adoption may significantly change the expected timescales for different ice processes in interstellar environments.

Small hydrocarbons, such as C2H, C3H and C3H2 are more abundant in photo-dissociation regions (PDRs) than expected based on gas-phase chemical models. To explore the hydrocarbon chemistry further, we observed a key intermediate species, the hydrocarbon ion l-C3H+, in the Horsehead PDR with the Plateau de Bure Interferometer at high-angular resolution (6''). We compare with previous observations of C2H and c-C3H2 at similar angular resolution and new gas-phase chemical model predictions to constrain the dominant formation mechanisms of small hydrocarbons in low-UV flux PDRs. We find that, at the peak of the HCO emission (PDR position), the measured l-C3H+, C2H and c-C3H2 abundances are consistent with current gas-phase model predictions. However, in the first PDR layers, at the 7.7 mum PAH band emission peak, which are more exposed to the radiation field and where the density is lower, the C2H and c-C3H2 abundances are underestimated by an order of magnitude. At this position, the l-C3H+ abundance is also underpredicted by the model but only by a factor of a few. In addition, contrary to the model predictions, l-C3H+ peaks further out in the PDR than the other hydrocarbons, C2H and c-C3H2. This cannot be explained by an excitation effect. Current gas-phase photochemical models thus cannot explain the observed abundances of hydrocarbons, in particular in the first PDR layers. Our observations are consistent with a top-down hydrocarbon chemistry, in which large polyatomic molecules or small carbonaceous grains are photo-destroyed into smaller hydrocarbon molecules/precursors.

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.

Complex organic molecules (COMs) are known to be abundant toward some low-mass young stellar objects (YSOs), but how these detections relate to typical COM abundance are not yet understood. We aim to constrain the frequency distribution of COMs during low-mass star formation, beginning with this pilot survey of COM lines toward six embedded YSOs using the IRAM 30m telescope. The sample was selected from the Spitzer c2d ice sample and covers a range of ice abundances. We detect multiple COMs, including CH3CN, toward two of the YSOs, and tentatively toward a third. Abundances with respect to CH3OH vary between 0.7 and 10%. This sample is combined with previous COM observations and upper limits to obtain a frequency distributions of CH3CN, HCOOCH3, CH3OCH3 and CH3CHO. We find that for all molecules more than 50% of the sample have detections or upper limits of 1-10% with respect to CH3OH. Moderate abundances of COMs thus appear common during the early stages of low-mass star formation. A larger sample is required, however, to quantify the COM distributions, as well as to constrain the origins of observed variations across the sample.

HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system.

Volatiles are compounds with low sublimation temperatures, and they make up most of the condensible mass in typical planet-forming environments. They consist of relatively small, often hydrogenated, molecules based on the abundant elements carbon, nitrogen and oxygen. Volatiles are central to the process of planet formation, forming the backbone of a rich chemistry that sets the initial conditions for the formation of planetary atmospheres, and act as a solid mass reservoir catalyzing the formation of planets and planetesimals. This growth has been driven by rapid advances in observations and models of protoplanetary disks, and by a deepening understanding of the cosmochemistry of the solar system. Indeed, it is only in the past few years that representative samples of molecules have been discovered in great abundance throughout protoplanetary disks - enough to begin building a complete budget for the most abundant elements after hydrogen and helium. The spatial distributions of key volatiles are being mapped, snow lines are directly seen and quantified, and distinct chemical regions within protoplanetary disks are being identified, characterized and modeled. Theoretical processes invoked to explain the solar system record are now being observationally constrained in protoplanetary disks, including transport of icy bodies and concentration of bulk condensibles. The balance between chemical reset - processing of inner disk material strong enough to destroy its memory of past chemistry, and inheritance - the chemically gentle accretion of pristine material from the interstellar medium in the outer disk, ultimately determines the final composition of pre-planetary matter. This chapter focuses on making the first steps toward understanding whether the planet formation processes that led to our solar system are universal.

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.

Based on astrochemical theory, the complex molecular composition around high-mass YSOs should evolve from the outer envelope in toward the central hot region as a sequence of temperature dependent chemical pathways are activated in ices and in the gas-phase. The resulting complex molecules have been broadly classified into three generations dependent on the temperature (<25, >25, and >100 K) required for formation. We combine IRAM 30m and Submillimeter Array observations to explore the spatial distribution of organic molecules around the high-mass young stellar object NGC 7538 IRS9, whose weak complex molecule emission previously escaped detection, quantifying the emission and abundance profiles of key organic molecules as a function of distance from the central protostar. We find that emission from N-bearing organics and saturated O-bearing organics present large increases in emission around 8000 AU and R<3000 AU, while O-bearing molecules and hydrocarbons do not. The increase in flux from some complex molecules in the envelope, around 8000 AU or 25 K, is consistent with recent model predictions of an onset of complex ice chemistry at 20-30 K. The emission increase for some molecules at R<3000 AU suggests the presence of a weak hot core, where thermal ice evaporation and hot gas-phase chemistry drives the chemistry. Complex organics thus form at all radii and temperatures around this protostar, but the composition changes dramatically as the temperature increases, which is used to constrain the chemical generation(s) to which different classes of molecule belong.

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.

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.

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.

The C/O ratio is predicted to regulate the atmospheric chemistry in hot Jupiters. Recent observations suggest that some exo-planets, e.g. Wasp 12- b, have atmospheric C/O ratios substantially different from the solar value of 0.54. In this paper we present a mechanism that can produce such atmospheric deviations from the stellar C/O ratio. In protoplanetary disks, different snowlines of oxygen- and carbon-rich ices, especially water and carbon monoxide, will result in systematic variations in the C/O ratio both in the gas and in the condensed phase. In particular, between the H2O and CO snowlines most oxygen is present in icy grains - the building blocks of planetary cores in the core accretion model - while most carbon remains in the gas-phase. This region is coincidental with the giant-planet forming zone for a range of observed protoplanetary disks. Based on standard core accretion models of planet formation, gas giants that sweep up most of their atmospheres from disk gas outside of the water snowline will have C/O?1, while atmospheres significantly contaminated by evaporating planetesimals will have stellar or sub-stellar C/O when formed at the same disk radius. The overall metallicity will also depend on the atmosphere formation mechanism, and exoplanetary atmospheric compositions may therefore provide constraints on where and how a specific planet formed.

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.

UV-induced photodesorption of ice is a non-thermal evaporation process that can explain the presence of cold molecular gas in a range of interstellar regions. Information on the average UV photodesorption yield of astrophysically important ices exists for broadband UV lamp experiments. UV fields around low-mass pre-main sequence stars, around shocks and in many other astrophysical environments are however often dominated by discrete atomic and molecular emission lines. It is therefore crucial to consider the wavelength dependence of photodesorption yields and mechanisms. In this work, for the first time, the wavelength-dependent photodesorption of pure CO ice is explored between 90 and 170 nm. The experiments are performed under ultra high vacuum conditions using tunable synchrotron radiation. Ice photodesorption is simultaneously probed by infrared absorption spectroscopy in reflection mode of the ice and by quadrupole mass spectrometry of the gas phase. The experimental results for CO reveal a strong wavelength dependence directly linked to the vibronic transition strengths of CO ice, implying that photodesorption is induced by electronic transition (DIET). The observed dependence on the ice absorption spectra implies relatively low photodesorption yields at 121.6 nm (Ly-alpha), where CO barely absorbs, compared to the high yields found at wavelengths coinciding with transitions into the first electronic state of CO (singulet Pi at 150 nm); the CO photodesorption rates depend strongly on the UV profiles encountered in different star formation environments.

Ice desorption affects the evolution of the gas-phase chemistry during the protostellar stage, and also determines the chemical composition of comets forming in circumstellar disks. From observations, most volatile species are found in H2O-dominated ices. The aim of this study is first to experimentally determine how entrapment of volatiles in H2O ice depends on ice thickness, mixture ratio and heating rate, and second, to introduce an extended three-phase model (gas, ice surface and ice mantle) to describe ice mixture desorption with a minimum number of free parameters. Thermal H2O:CO2 ice desorption is investigated in temperature programmed desorption experiments of thin (10 - 40 ML) ice mixtures under ultra-high vacuum conditions. Desorption is simultaneously monitored by mass spectrometry and reflection-absorption infrared spectroscopy. The H2O:CO2 experiments are complemented with selected H2O:CO, and H2O:CO2:CO experiments. The results are modeled with rate equations that connect the gas, ice surface and ice mantle phases through surface desorption and mantle-surface diffusion. The fraction of trapped CO2 increases with ice thickness (10 - 32 ML) and H2O:CO2 mixing ratio (5:1 - 10:1), but not with one order of magnitude different heating rates. The fraction of trapped CO2 is 44 - 84 % with respect to the initial CO2 content for the investigated experimental conditions. This is reproduced quantitatively by the extended three-phase model that is introduced here. The H2O:CO and H2O:CO2:CO experiments are consistent with the H2O:CO2 desorption trends, suggesting that the model can be used for other ice species found in the interstellar medium to significantly improve the parameterization of ice desorption.

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

Ices regulate much of the chemistry during star formation and account for up to 80% of the available oxygen and carbon. In this paper, we use the Spitzer c2d ice survey, complimented with data sets on ices in cloud cores and high-mass protostars, to determine standard ice abundances and to present a coherent picture of the evolution of ices during low- and high-mass star formation. The median ice composition H2O:CO:CO2:CH3OH:NH3:CH4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 toward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. In the low-mass sample, the ice abundances with respect to H2O of CH4, NH3, and the component of CO2 mixed with H2O typically vary by <25%, indicative of co-formation with H2O. In contrast, some CO and CO2 ice components, XCN and CH3OH vary by factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, consistent with its OCN- identification. The origin(s) of the different levels of ice abundance variations are constrained by comparing ice inventories toward different types of protostars and background stars, through ice mapping, analysis of cloud-to-cloud variations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase is driven by hydrogenation of atoms, which results in a H2O-dominated ice. At later prestellar times, CO freezes out and variations in CO freeze-out levels and the subsequent CO-based chemistry can explain most of the observed ice abundance variations. The last important ice evolution stage is thermal and UV processing around protostars, resulting in CO desorption, ice segregation and formation of complex organic molecules. The distribution of cometary ice abundances are consistent with with the idea that most cometary ices have a protostellar origin.

Icy grain mantles are commonly observed through infrared spectroscopy toward dense clouds, cloud cores, protostellar envelopes and protoplanetary disks. Up to 80% of the available oxygen, carbon and nitrogen are found in such ices; the most common ice constituents - H2O, CO2 and CO - are second in abundance only to H2 in many star forming regions. In addition to being a molecular reservoir, ice chemistry is responsible for much of the chemical evolution from H2O to complex, prebiotic molecules. Combining the existing ISO, Spitzer, VLT and Keck ice data results in a large sample of ice sources (\sime80) that span all stages of star formation and a large range of protostellar luminosities (<0.1-105 L\odot). Here we summarize the different techniques that have been applied to mine this ice data set on information on typical ice compositions in different environments and what this implies about how ices form and evolve during star and planet formation. The focus is on how to maximize the use of empirical constraints from ice observations, followed by the application of information from experiments and models. This strategy is used to identify ice bands and to constrain which ices form early during cloud formation, which form later in the prestellar core and which require protostellar heat and/or UV radiation to form. The utility of statistical tests, survival analysis and ice maps is highlighted; the latter directly reveals that the prestellar ice formation takes place in two phases, associated with H2O and CO ice formation, respectively, and that most protostellar ice variation can be explained by differences in the prestellar CO ice formation stage. Finally, special attention is paid to the difficulty of observing complex ices directly and how gas observations, experiments and models help in constraining this ice chemistry stage.

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.

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.

The purpose of the 2010 NASA Laboratory Astrophysics Workshop (LAW) was, as given in the Charter from NASA, "to provide a forum within which the scientific community can review the current state of knowledge in the field of Laboratory Astrophysics, assess the critical data needs of NASA's current and future Space Astrophysics missions, and identify the challenges and opportunities facing the field as we begin a new decade". LAW 2010 was the fourth in a roughly quadrennial series of such workshops sponsored by the Astrophysics Division of the NASA Science Mission Directorate. In this White Paper, we report the findings of the workshop.

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

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.

UV irradiation of simple ices is proposed to efficiently produce complex organic species during star- and planet-formation. Through a series of laboratory experiments, we investigate the effects of the H2O concentration, the dominant ice constituent in space, on the photochemistry of more volatile species, especially CH4, in ice mixtures. In the experiments, thin (~40 ML) ice mixtures, kept at 20-60 K, are irradiated under ultra-high vacuum conditions with a broad-band UV hydrogen discharge lamp. Photodestruction cross sections of volatile species (CH4 and NH3) and production efficiencies of new species (C2H6, C2H4, CO, H2CO, CH3OH, CH3CHO and CH3CH2OH) in water-containing ice mixtures are determined using reflection-absorption infrared spectroscopy during irradiation and during a subsequent slow warm-up. The four major effects of increasing the H2O concentration are 1) an increase of the destruction efficiency of the volatile mixture constituent by up to an order of magnitude due to a reduction of back reactions following photodissociation, 2) a shift to products rich in oxygen e.g. CH3OH and H2CO, 3) trapping of up to a factor of five more of the formed radicals in the ice and 4) a disproportional increase in the diffusion barrier for the OH radical compared to the CH3 and HCO radicals. The radical diffusion temperature dependencies are consistent with calculated H2O-radical bond strengths. All the listed effects are potentially important for the production of complex organics in H2O-rich icy grain mantles around protostars and should thus be taken into account when modeling ice chemistry.

NH3 and CH3OH are key molecules in astrochemical networks leading to the formation of more complex N- and O-bearing molecules, such as CH3CN and HCOOCH3. Despite a number of recent studies, little is known about their abundances in the solid state. (...) In this work, we investigate the ~ 8-10 micron region in the Spitzer IRS (InfraRed Spectrograph) spectra of 41 low-mass young stellar objects (YSOs). These data are part of a survey of interstellar ices in a sample of low-mass YSOs studied in earlier papers in this series. We used both an empirical and a local continuum method to correct for the contribution from the 10 micron silicate absorption in the recorded spectra. In addition, we conducted a systematic laboratory study of NH3- and CH3OH-containing ices to help interpret the astronomical spectra. We clearly detect a feature at ~9 micron in 24 low-mass YSOs. Within the uncertainty in continuum determination, we identify this feature with the NH3 nu_2 umbrella mode, and derive abundances with respect to water between ~2 and 15%. Simultaneously, we also revisited the case of CH3OH ice by studying the nu_4 C-O stretch mode of this molecule at ~9.7 micron in 16 objects, yielding abundances consistent with those derived by Boogert et al. 2008 (hereafter paper I) based on a simultaneous 9.75 and 3.53 micron data analysis. Our study indicates that NH3 is present primarily in H2O-rich ices, but that in some cases, such ices are insufficient to explain the observed narrow FWHM. The laboratory data point to CH3OH being in an almost pure methanol ice, or mixed mainly with CO or CO2, consistent with its formation through hydrogenation on grains. Finally, we use our derived NH3 abundances in combination with previously published abundances of other solid N-bearing species to find that up to 10-20 % of nitrogen is locked up in known ices.

Gas-phase complex organic molecules have been detected toward a range of high- and low-mass star-forming regions at abundances which cannot be explained by any known gas-phase chemistry. Recent laboratory experiments show that UV irradiation of CH3OH-rich ices may be an important mechanism for producing complex molecules and releasing them into the gas-phase. To test this ice formation scenario we mapped the B1-b dust core and nearby protostar in CH3OH gas using the IRAM 30m telescope to identify locations of efficient non-thermal ice desorption. We find three CH3OH abundance peaks tracing two outflows and a quiescent region on the side of the core facing the protostar. The CH3OH gas has a rotational temperature of ~10 K at all locations. The quiescent CH3OH abundance peak and one outflow position were searched for complex molecules. Narrow, 0.6-0.8 km s-1 wide, HCOOCH3 and CH3CHO lines originating in cold gas are clearly detected, CH3OCH3 is tentatively detected and C2H5OH and HOCH2CHO are undetected toward the quiescent core, while no complex molecular lines were found toward the outflow. The core abundances with respect to CH3OH are ~2.3% and 1.1% for HCOOCH3 and CH3CHO, respectively, and the upper limits are 0.7-1.1%, which is similar to most other low-mass sources. The observed complex molecule characteristics toward B1-b and the pre-dominance of HCO-bearing species suggest a cold ice (below 25 K, the sublimation temperature of CO) formation pathway followed by non-thermal desorption through e.g. UV photons traveling through outflow cavities. The observed complex gas composition together with the lack of any evidence of warm gas-phase chemistry provide clear evidence of efficient complex molecule formation in cold interstellar ices.

(Abridged) Gas-phase complex organic molecules are commonly detected in the warm inner regions of protostellar envelopes. Recent models show that photochemistry in ices followed by desorption may explain the observed abundances. This study aims to experimentally quantify the broad-band UV-induced production rates of complex organics in CH3OH-rich ices at 20-70 K under ultra-high vacuum conditions. The reaction products are mainly identified by RAIRS and TPD experiments. Complex organics are readily formed in all experiments, both during irradiation and during a slow warm-up of the ices to 200 K after the UV lamp is turned off. The relative abundances of photoproducts depend on the UV fluence, the ice temperature, and whether pure CH3OH ice or CH3OH:CH4/CO ice mixtures are used. C2H6, CH3CHO, CH3CH2OH, CH3OCH3, HCOOCH3, HOCH2CHO and (CH2OH)2 are all detected in at least one experiment. The derived product-formation yields and their dependences on different experimental parameters, such as the initial ice composition, are used to estimate the CH3OH photodissociation branching ratios in ice and the relative diffusion barriers of the formed radicals. The experiments show that ice photochemistry in CH3OH ices is efficient enough to explain the observed abundances of complex organics around protostars and that ratios of complex molecules can be used to constrain their formation pathway.

A longstanding problem in astrochemistry is how molecules can be maintained in the gas phase in dense inter- and circumstellar regions. Photodesorption is a non-thermal desorption mechanism, which may explain the small amounts of observed cold gas in cloud cores and disk mid-planes. This paper aims to determine the UV photodesorption yields and to constrain the photodesorption mechanisms of three astrochemically relevant ices: CO, N2 and CO2. In addition, the possibility of co-desorption in mixed and layered CO:N2 ices is explored. The ice photodesorption is studied experimentally under ultra high vacuum conditions and at 15-60 K using a hydrogen discharge lamp (7-10.5 eV). The ice desorption during irradiation is monitored by reflection absorption infrared spectroscopy of the ice and simultaneous mass spectrometry of the desorbed molecules. Both the UV photodesorption yields per incident photon and the photodesorption mechanisms are molecule specific. CO photodesorbs without dissociation from the surface layer of the ice. N2, which lacks an electronic transition in this wavelength range, has a photodesorption yield that is more than an order of magnitude lower. This yield increases significantly due to co-desorption when N2 is mixed in with or layered on top of CO ice. CO2 photodesorbs through dissociation and subsequent recombination from the top 10 layers of the ice. At low temperatures (15-18 K) the derived photodesorption yields are 2.7x10^-3 and <2x10^-4 molecules photon-1 for pure CO and N2, respectively. The CO2 photodesorption yield is 1.2x10^-3x(1-e^(-X/2.9)) + 1.1x10^-3x(1-e^(-X/4.6)) molecules photon-1, where X is the ice thickness in monolayers and the two parts of the expression represent a CO2 and CO photodesorption pathway.

Up to 90% of the chemical reactions during star formation occurs on ice surfaces, probably including the formation of complex organics. Only the most abundant ice species are however observed directly by infrared spectroscopy. This study aims to develop an indirect observational method of ices based on non-thermal ice desorption in the colder part of protostellar envelopes. For that purpose the IRAM 30m telescope was employed to observe two molecules that can be detected both in the gas and the ice, CH3 OH and HNCO, toward 4 low mass embedded protostars. Their respective gas-phase column densities are determined using rotational diagrams. The relationship between ice and gas phase abundances is subsequently determined. The observed gas and ice abundances span several orders of magnitude. Most of the CH3OH and HNCO gas along the lines of sight is inferred to be quiescent from the measured line widths and the derived excitation temperatures, and hence not affected by thermal desorption close to the protostar or in outflow shocks. The measured gas to ice ratio of ~10-4 agrees well with model predictions for non-thermal desorption under cold envelope conditions and there is a tentative correlation between ice and gas phase abundances. This indicates that non-thermal desorption products can serve as a signature of the ice composition. A larger sample is however necessary to provide a conclusive proof of concept.

Gaseous H2O has been detected in several cold astrophysical environments, where the observed abundances cannot be explained by thermal desorption of H2O ice or by H2O gas phase formation. These observations hence suggest an efficient non-thermal ice desorption mechanism. Here, we present experimentally determined UV photodesorption yields of H2O and D2O ice and deduce their photodesorption mechanism. The ice photodesorption is studied under ultra high vacuum conditions and at astrochemically relevant temperatures (18-100 K) using a hydrogen discharge lamp (7-10.5 eV), which simulates the interstellar UV field. The ice desorption during irradiation is monitored using reflection absorption infrared spectroscopy of the ice and simultaneous mass spectrometry of the desorbed species. The photodesorption yield per incident photon is identical for H2O and D2O and depends on both ice thickness and temperature. For ices thicker than 8 monolayers the photodesorption yield Y is linearly dependent on temperature due to increased diffusion of ice species such that Y(T) = 1E-3(1.3+0.032*T) UV photon-1, with a 60% uncertainty for the absolute yield. The increased diffusion also results in an increasing H2O:OH desorption product ratio with temperature. The yield does not depend on the substrate, the UV photon flux or the UV fluence. The yield is also independent on the initial ice structure since UV photons efficiently amorphize H2O ice. The results are consistent with theoretical predictions of H2O photodesorption and partly in agreement with a previous experimental study. Applying the experimentally determined yield to a Herbig Ae/Be star+disk model shows that UV photodesorption of ices increases the H2O content by orders of magnitude in the disk surface region compared to models where non-thermal desorption is ignored.

CH4 is proposed to be the starting point of a rich organic chemistry. Solid CH4 abundances have previously been determined mostly toward high mass star forming regions. Spitzer/IRS now provides a unique opportunity to probe solid CH4 toward low mass star forming regions as well. Infrared spectra from the Spitzer Space Telescope are presented to determine the solid CH4 abundance toward a large sample of low mass young stellar objects. 25 out of 52 ice sources in the $c2d$ (cores to disks) legacy have an absorption feature at 7.7 um, attributed to the bending mode of solid CH4. The solid CH4 / H2O abundances are 2-8%, except for three sources with abundances as high as 11-13%. These latter sources have relatively large uncertainties due to small total ice column densities. Toward sources with H2O column densities above 2E18 cm-2, the CH4 abundances (20 out of 25) are nearly constant at 4.7+/-1.6%. Correlation plots with solid H2O, CH3OH, CO2 and CO column densities and abundances relative to H2O reveal a closer relationship of solid CH4 with CO2 and H2O than with solid CO and CH3OH. The inferred solid CH4 abundances are consistent with models where CH4 is formed through sequential hydrogenation of C on grain surfaces. Finally the equal or higher abundances toward low mass young stellar objects compared with high mass objects and the correlation studies support this formation pathway as well, but not the two competing theories: formation from CH3OH and formation in gas phase with subsequent freeze-out.

With the goal to study the physical and chemical evolution of ices in solar-mass systems, a spectral survey is conducted of a sample of 41 low luminosity YSOs using 3-38 um Spitzer and ground-based spectra. The long-known 6.0 and 6.85 um bands are detected toward all sources, with the Class 0-type YSOs showing the deepest bands ever observed. In almost all sources the 6.0 um band is deeper than expected from the bending mode of pure solid H2O. The depth and shape variations of the remaining 5-7 um absorption indicate that it consists of 5 independent components, which, by comparison to laboratory studies, must be from at least 8 different carriers. Simple species are responsible for much of the absorption in the 5-7 um region, at abundances of 1-30% for CH3OH, 3-8% for NH3, 1-5% for HCOOH, ~6% for H2CO, and ~0.3% for HCOO- with respect to solid H2O. The 6.85 um band likely consists of one or two carriers, of which one is less volatile than H2O because its abundance relative to H2O is enhanced at lower H2O/tau_9.7 ratios. It does not survive in the diffuse interstellar medium (ISM), however. The similarity of the 6.85 um bands for YSOs and background stars indicates that its carrier(s) must be formed early in the molecular cloud evolution. If an NH4+ salt is the carrier its abundance with respect to solid H2O is typically 7%, and low temperature acid-base chemistry or cosmic ray induced reactions must have been involved in its formation. Possible origins are discussed for the carrier of an enigmatic, very broad absorption between 5 and 8 um. Finally, all the phenomena observed for ices toward massive YSOs are also observed toward low mass YSOs, indicating that processing of the ices by internal ultraviolet radiation fields is a minor factor in the early chemical evolution of the ices. [abridged]

A laboratory study on the band profiles and band strengths of H2O in CO ice, and vice versa, is presented and interpreted in terms of two models. The results show that a mutual interaction takes place between the two species in the solid, which alters the band positions and band strengths. It is found that the band strengths of the H2O bulk stretch, bending and libration vibrational bands decrease linearly by a factor of up to 2 when the CO concentration is increased from 0 to 80%. By contrast, the band strength of the free OH stretch increases linearly. The results are compared to a recently performed quantitative study on H2O:CO2 ice mixtures. It is shown that for mixing ratios of 1:0.5 H2O:X and higher, the H2O bending mode offers a good tracer to distinguish between CO2 or CO in H2O ice. Additionally, it is found that the band strength of the CO fundamental remains constant when the water concentration is increased in the ice. The integrated absorbance of the 2152 cm-1 CO feature, with respect to the total integrated CO absorption feature, is found to be a good indicator of the degree of mixing of CO in the H2O:CO laboratory ice system. From the change in the H2O absorption band strength in laboratory ices upon mixing we conclude that astronomical water ice column densities on various lines of sight can be underestimated by up to 25% if significant amounts of CO and CO2 are mixed in.

This paper presents Spitzer-IRS spectroscopy of the CO2 15.2 micron bending mode toward a sample of 50 embedded low-mass stars in nearby star-forming clouds, taken mostly from the ``Cores to Disks (c2d)'' Legacy program. The average abundance of solid CO2 relative to water in low-mass protostellar envelopes is 0.32 +/- 0.02, significantly higher than that found in quiescent molecular clouds and in massive star forming regions. It is found that a decomposition of all the observed CO2 bending mode profiles requires a minimum of five unique components. Roughly 2/3 of the CO2 ice is found in a water-rich environment, while most of the remaining 1/3 is found in a CO environment. Ground-based observations of solid CO toward a large subset of the c2d sample are used to further constrain the CO2:CO component and suggest a model in which low-density clouds form the CO2:H2O component and higher density clouds form the CO2:CO ice during and after the freeze-out of gas-phase CO. It is suggested that the subsequent evolution of the CO2 and CO profiles toward low-mass protostars, in particular the appearance of the splitting of the CO2 bending mode due to pure, crystalline CO2, is first caused by distillation of the CO2:CO component through evaporation of CO due to thermal processing to ~20-30 K in the inner regions of infalling envelopes. The formation of pure CO2 via segregation from the H2O rich mantle may contribute to the band splitting at higher levels of thermal processing (>50 K), but is harder to reconcile with the physical structure of protostellar envelopes around low-luminosity objects.

At the high densities and low temperatures found in star forming regions, all molecules other than H2 should stick on dust grains on timescales shorter than the cloud lifetimes. Yet these clouds are detected in the millimeter lines of gaseous CO. At these temperatures, thermal desorption is negligible and hence a non-thermal desorption mechanism is necessary to maintain molecules in the gas phase. Here, the first laboratory study of the photodesorption of pure CO ice under ultra high vacuum is presented, which gives a desorption rate of 3E-3 CO molecules per UV (7-10.5 eV) photon at 15 K. This rate is factors of 1E2-1E5 larger than previously estimated and is comparable to estimates of other non-thermal desorption rates. The experiments constrains the mechanism to a single photon desorption process of ice surface molecules. The measured efficiency of this process shows that the role of CO photodesorption in preventing total removal of molecules in the gas has been underestimated.

H2O is the most abundant component of astrophysical ices. In most lines of sight it is not possible to fit both the H2O 3 um stretching, the 6 um bending and the 13 um libration band intensities with a single pure H2O spectrum. Recent Spitzer observations have revealed CO2 ice in high abundances and it has been suggested that CO2 mixed into H2O ice can affect relative strengths of the 3 um and 6 um bands. We used laboratory infrared transmission spectroscopy of H2O:CO2 ice mixtures to investigate the effects of CO2 on H2O ice spectral features at 15-135 K. We find that the H2O peak profiles and band strengths are significantly different in H2O:CO2 ice mixtures compared to pure H2O ice. In all H2O:CO2 mixtures, a strong free-OH stretching band appears around 2.73 um, which can be used to put an upper limit on the CO2 concentration in the H2O ice. The H2O bending mode profile also changes drastically with CO2 concentration; the broad pure H2O band gives way to two narrow bands as the CO2 concentration is increased. This makes it crucial to constrain the environment of H2O ice to enable correct assignments of other species contributing to the interstellar 6 um absorption band. The amount of CO2 present in the H2O ice of B5:IRS1 is estimated by simultaneously comparing the H2O stretching and bending regions and the CO2 bending mode to laboratory spectra of H2O, CO2, H2O:CO2 and HCOOH.

We present Temperature Programmed Desorption (TPD) experiments of CO and N2 ices in pure, layered and mixed morphologies at various ice "thicknesses" and abundance ratios as well as simultaneously taken Reflection Absorption Infrared Spectra (RAIRS) of CO. A kinetic model has been developed to constrain the binding energies of CO and N2 in both pure and mixed environments and to derive the kinetics for desorption, mixing and segregation. For mixed ices N2 desorption occurs in a single step whereas for layered ices it proceeds in two steps, one corresponding to N2 desorption from a pure N2 ice environment and one corresponding to desorption from a mixed ice environment. The latter is dominant for astrophysically relevant ice "thicknesses". The ratio of the binding energies, Rbe, for pure N2 and CO is found to be 0.936 +/- 0.03, and to be close to 1 for mixed ice fractions. The model is applied to astrophysically relevant conditions for cold pre-stellar cores and for protostars which start to heat their surroundings. The importance of treating CO desorption with zeroth rather than first order kinetics is shown. The experiments also provide lower limits of 0.87 +/- 0.05 for the sticking probabilities of CO-CO, N2-CO and N2-N2 ices at 14 K. The combined results from the desorption experiments, the kinetic model, and the sticking probability data lead to the conclusion that these solid-state processes of CO and N2 are very similar under astrophysically relevant conditions. This conclusion affects the explanations for the observed anti-correlations of gaseous CO and N2H+ in pre-stellar and protostellar cores.

Absorption and emission lines originating from the nu3 C-H stretching manifold of gas phase CH4 were discovered in the high resolution (R=25,000) infrared L band spectrum along the line of sight toward NGC 7538 : IRS9. These observations provide a diagnostic of the complex dynamics and chemistry in a massive star forming region. The line shapes resemble P Cygni profiles with the absorption and emission components shifted by ~7 km/s with respect to the systemic velocity. Similar velocity components were observed in CO at 4.7 um, but in contrast to CH4, the CO shows deep absorption due to a high velocity outflow as well as absorption at the systemic velocity due to the cold outer envelope. It is concluded that the gas phase CH4 abundance varies by an order of magnitude in this line of sight: it is low in the envelope and the outflow (X[CH4]<0.4e-6), and at least a factor of 10 larger in the central core. The discovery of solid CH4 in independent ground and space based data sets shows that methane is nearly entirely frozen onto grains in the envelope. It thus appears that CH4 is formed by grain surface reactions, evaporates into the gas phase in the warm inner regions of protostellar cores and is efficiently destroyed in shocks related to outflows.