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, encompass
ing 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 Hyas HCO$^+$ and N$_2$H$^+$ emission are best fit by a moderately hard X-ray spectra, as would be expected during the flaring state of the star, and a low CR ionization rate, $zeta_{rm CR}lesssim10^{-19}$ s$^{-1}$. This low CR rate is the first indication of the presence of CR exclusion by winds and/or magnetic fields in an actively accreting T Tauri disk system. With this new constraint, our best fit ionization structure predicts a low turbulence dead-zone extending from the inner edge of the disk out to $50-65$ AU. This region coincides with an observed concentration of millimeter grains, and we propose that the inner region of TW Hya is a dust (and possibly planet) growth factory as predicted by previous theoretical work.
We present high angular resolution Submillimeter Array observations ofthe outbursting Jupiter family comet 17P/Holmes on 2007 October 26-29, achieving a spatial resolution of 2.5, or ~3000 km at the comet distance. The observations resulted in detect
ions of the rotational lines CO 3-2, HCN 4-3, H$^{13}$CN 4-3, CS 7-6, H$_2$CO 3$_{1,2}$-2$_{1,1}$, H$_2$S 2$_{2,0}$-2$_{1,1}$, and multiple CH$_3$OH lines, along with the associated dust continuum at 221 and 349 GHz. The continuum has a spectral index of 2.7$pm$0.3, slightly steeper than blackbody emission from large dust particles. From the imaging data, we identify two components in the molecular emission. One component is characterized by a relatively broad line width (~1 km s$^{-1}$ FWHM) exhibiting a symmetric outgassing pattern with respect to the nucleus position. The second component has a narrower line width (<0.5 km s$^{-1}$ FWHM) with the line center red-shifted by 0.1-0.2 km s$^{-1}$ (cometocentric frame), and shows a velocity shift across the nucleus position with the position angle gradually changing from 66 to 30 degrees within the four days of observations. We determine distinctly different CO/HCN ratios for each of the components. For the broad-line component we find CO/HCN <7, while in the narrow-line component, CO/HCN = 40$pm$5. We hypothesize that the narrow-line component originates from the ice grain halo found in near-nucleus photometry, believed to be created by sublimating recently released ice grains around the nucleus during the outburst. In this interpretation, the high CO/HCN ratio of this component reflects the more pristine volatile composition of nucleus material released in the outburst.
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 aroun
d TW Hya, an analog of the solar nebula, using high spatial and spectral resolution Atacama Large Millimeter/Submillimeter Array (ALMA) observations of N2H+, a reactive ion present in large abundance only where CO is frozen out. The N2H+ emission is distributed in a large ring, with an inner radius that matches CO snow line model predictions. The extracted CO snow line radius of ~ 30 AU helps to assess models of the formation dynamics of the Solar System, when combined with measurements of the bulk composition of planets and comets.
We 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 1
63296 at 0.8 resolution. The emission is consistent with that arising from a Keplerian rotating disk. Two additional c-C3H2 transitions are also tentatively detected, bolstering the identification of this species, but with insufficient signal-to-noise ratio to constrain the spatial distribution. Using a previously developed model for the physical structure of this disk, we fit a radial power-law distribution model to the c-C3H2 6-5 emission and find that c-C3H2 is present in a ring structure from an inner radius of about 30 AU to an outer radius of about 165 AU. The column density is estimated to be 1e12-1e13 cm-2. The clear detection and intriguing ring structure suggest that c-C3H2 has the potential to become a useful probe of radiation penetration in disks.
We present Submillimeter Array observations of H2CO and N2H+ emission in the disks around the T Tauri star TW Hya and the Herbig Ae star HD 163296 at 2-6 resolution and discuss the distribution of these species with respect to CO freeze-out. The H2CO
and N2H+ emission toward HD 163296 does not peak at the continuum emission center that marks the stellar position but is instead significantly offset. Using a previously developed model for the physical structure of this disk, we show that the H2CO observations are reproduced if H2CO is present predominantly in the cold outer disk regions. A model where H2CO is present only beyond the CO snow line (estimated at a radius of 160 AU) matches the observations well. We also show that the average H2CO excitation temperature, calculated from two transitions of H2CO observed in these two disks and a larger sample of disks around T Tauri stars in the DISCS (the Disk Imaging Survey of Chemistry with SMA) program, is consistent with the CO freeze-out temperature of 20 K. In addition, we show that N2H+ and H2CO line fluxes in disks are strongly correlated, indicative of co-formation of these species across the sample. Taken together, these results imply that H2CO and N2H+ are generally present in disks only at low temperatures where CO depletes onto grains, consistent with fast destruction of N2H+ by gas-phase CO, and in situ formation of H2CO through hydrogenation of CO ice. In this scenario H2CO, CH3OH and N2H+ emission in disks should appear as rings with the inner edge at the CO midplane snow line. This prediction can be tested directly using observations from ALMA with higher resolution and better sensitivity.
We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lin
es exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.
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.
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 t
he 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.
