BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Chunhua Qi

Chunhua Qi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAphysics.space-phastro-phastro-ph.IMphysics.ins-det

Catalog footprint

What is connected

50works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2026arXiv

Isotopic Ratios in the Disk of HD 163296

Isotopic abundance ratios in protoplanetary disks are critical for understanding volatile inheritance and chemical evolution in planet-forming environments. We present Atacama Large Millimeter/submillimeter Array observations of the rare isotopologue 13C18O(2-1) at approximately 0.3 arcsec resolution from the disk around the Herbig Ae star HD 163296, combined with archival observations of C17O(2-1), C18O(1-0), and C17O(1-0), to empirically constrain carbon and oxygen isotopic ratios without detailed disk modeling. Both the C17O/13C18O(2-1) and C18O/C17O(1-0) flux ratios rise sharply across the CO snowline and flatten beyond 1.5 arcsec (r >= 150 au), where the emission becomes optically thin. This transition, reflecting a steep drop in CO column density set by the disk's thermal structure, makes HD 163296 an optimal case for isotopic analysis. Using beam-averaged intensities of the four transitions measured in this optically thin region, we derive isotopic ratios of 12C/13C = 75.3 (+14.7/-11.4) and 18O/17O = 3.28 (+0.31/-0.26), both consistent with local interstellar medium values. The 16O/18O ratio remains weakly constrained due to moderate optical depth in the C18O(1-0) line and degeneracy with CO column density. These results demonstrate that rare CO isotopologues can provide robust, empirical constraints on isotopic ratios in disks when sharp structural transitions allow for the identification of optically thin regions, and establish HD 163296 as a benchmark for extending such studies to other systems with resolved snowline structures.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Tracing Pebble Drift History in Two Protoplanetary Disks with CO Enhancement

Pebble drift is an important mechanism for supplying the materials needed to build planets in the inner region of protoplanetary disks. Thus, constraining pebble drift's timescales and mass flux is essential to understanding planet formation history. Current pebble drift models suggest pebble fluxes can be constrained from the enhancement of gaseous volatile abundances when icy pebbles sublimate after drifting across key snowlines. In this work, we present ALMA observations of spatially resolved $^{13}$C$^{18}$O J=2-1 line emission inside the midplane CO snowline of the HD 163296 and MWC 480 protoplanetary disks. We use radiative transfer and thermochemical models to constrain the spatial distribution of CO gas column density. We find that both disks display centrally peaked CO abundance enhancement of up to ten times of ISM abundance levels. For HD 163296 and MWC 480, the inferred enhancements require 250-350 and 480-660 Earth Masses of pebbles to have drifted across their CO snowlines, respectively. These ranges fall within cumulative pebble mass flux ranges to grow gas giants in the interior to the CO snowline. The centrally peaked CO enhancement is unexpected in current pebble drift models, which predict CO enhancement peaks at the CO snowline or is uniform inside the snowline. We propose two hypotheses to explain the centrally-peaked CO enhancement, including a large CO desorption distance and CO trapped in water ice. By testing both hypotheses with the 1D gas and dust evolution code chemcomp, we find that volatile trapping (about 30\%) best reproduces the centrally peaked CO enhancement observed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Detection of HC18O+ in a protoplanetary disk: exploring oxygen isotope fractionation of CO

The oxygen isotope fractionation scenario, which has been developed to explain the oxygen isotope anomaly in the solar system materials, predicts that CO gas is depleted in 18O in protoplanetary disks, where segregation between solids and gas inside disks had already occurred. Based on ALMA observations, we report the first detection of HC18O+(4-3) in a Class II protoplanetary disk (TW Hya). This detection allows us to explore the oxygen isotope fractionation of CO in the TW Hya disk from optically thin HCO+ isotopologues as a proxy of optically thicker CO isotopologues. Using the H13CO+(4-3) data previously obtained with SMA, we find that the H13CO+/HC18O+ ratio in the central <100 au regions of the disk is 10.3 +- 3.2. We construct a chemical model of the TW Hya disk with carbon and oxygen isotope fractionation chemistry, and estimate the conversion factor from H13CO+/HC18O+ to 13CO/C18O. With the conversion factor (= 0.8), the 13CO/C18O ratio is estimated to be 8.3 +- 2.6, which is consistent with the elemental abundance ratio in the local ISM (8.1 +- 0.8) within error margin. Then there is no clear evidence of 18O depletion in CO gas of the disk, although we could not draw any robust conclusion due to large uncertainties. In conclusion, optically thin lines of HCO+ isotopologues are useful tracers of CO isotopic ratios, which are hardly constrained directly from optically thick lines of CO isotopologues. Future higher sensitivity observations of H13CO+ and HC18O+ would be able to allow us to better constrain the oxygen fractionation in the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS). A Circumplanetary Disk Candidate in Molecular Line Emission in the AS 209 Disk

We report the discovery of a circumplanetary disk (CPD) candidate embedded in the circumstellar disk of the T Tauri star AS 209 at a radial distance of about 200 au (on-sky separation of 1."4 from the star at a position angle of $161^\circ$), isolated via $^{13}$CO $J=2-1$ emission. This is the first instance of CPD detection via gaseous emission capable of tracing the overall CPD mass. The CPD is spatially unresolved with a $117\times82$ mas beam and manifests as a point source in $^{13}$CO, indicating that its diameter is $\lesssim14$ au. The CPD is embedded within an annular gap in the circumstellar disk previously identified using $^{12}$CO and near-infrared scattered light observations, and is associated with localized velocity perturbations in $^{12}$CO. The coincidence of these features suggests that they have a common origin: an embedded giant planet. We use the $^{13}$CO intensity to constrain the CPD gas temperature and mass. We find that the CPD temperature is $\gtrsim35$ K, higher than the circumstellar disk temperature at the radial location of the CPD, 22 K, suggesting that heating sources localized to the CPD must be present. The CPD gas mass is $\gtrsim 0.095 M_{\rm Jup} \simeq 30 M_{\rm Earth}$ adopting a standard $^{13}$CO abundance. From the non-detection of millimeter continuum emission at the location of the CPD ($3σ$ flux density $\lesssim26.4~μ$Jy), we infer that the CPD dust mass is $\lesssim 0.027 M_{\rm Earth} \simeq 2.2$ lunar masses, indicating a low dust-to-gas mass ratio of $\lesssim9\times10^{-4}$. We discuss the formation mechanism of the CPD-hosting giant planet on a wide orbit in the framework of gravitational instability and pebble accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Three Dimensional View of Gomez's Hamburger

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An ALMA Survey of H$_2$CO in Protoplanetary Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An evolutionary study of volatile chemistry in protoplanetary disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An unbiased ALMA spectral survey of the LkCa 15 and MWC 480 protoplanetary disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring turbulent motion in planet-forming disks with ALMA: A detection around DM Tau and non-detections around MWC 480 and V4046 Sgr

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The TW Hya Rosetta Stone Project I: Radial and vertical distributions of DCN and DCO+

Molecular D/H ratios are frequently used to probe the chemical past of Solar System volatiles. Yet it is unclear which parts of the Solar Nebula hosted an active deuterium fractionation chemistry. To address this question, we present 0".2-0".4 ALMA observations of DCO+ and DCN 2-1, 3-2 and 4-3 towards the nearby protoplanetary disk around TW Hya, taken as part of the TW Hya Rosetta Stone project, augmented with archival data. DCO+ is characterized by an excitation temperature of ~40 K across the 70 au radius pebble disk, indicative of emission from a warm, elevated molecular layer. Tentatively, DCN is present at even higher temperatures. Both DCO+ and DCN present substantial emission cavities in the inner disk, while in the outer disk the DCO+ and DCN morphologies diverge: most DCN emission originates from a narrow ring peaking around 30~au, with some additional diffuse DCN emission present at larger radii, while DCO+ is present in a broad structured ring that extends past the pebble disk. Based on parametric disk abundance models, these emission patterns can be explained by a near-constant DCN abundance exterior to the cavity, and an increasing DCO+ abundance with radius. There appears to be an active deuterium fractionation chemistry in multiple disk regions around TW Hya, but not in the cold planetesimal-forming midplane and in the inner disk. More observations are needed to explore whether deuterium fractionation is actually absent in these latter regions, and if its absence is a common feature, or something peculiar to the old TW Hya disk.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The TW Hya Rosetta Stone Project II: Spatially resolved emission of formaldehyde hints at low-temperature gas-phase formation

Formaldehyde (H$_2$CO) is an important precursor to organics like methanol (CH$_3$OH). It is important to understand the conditions that produce H$_2$CO and prebiotic molecules during star and planet formation. H$_2$CO possesses both gas-phase and solid-state formation pathways, involving either UV-produced radical precursors or CO ice and cold ($\lesssim 20$ K) dust grains. To understand which pathway dominates, gaseous H$_2$CO's ortho-to-para ratio (OPR) has been used as a probe, with a value of 3 indicating "warm" conditions and $<3$ linked to cold formation in the solid-state. We present spatially resolved ALMA observations of multiple ortho- and para-H$_2$CO transitions in the TW Hya protoplanetary disk to test H$_2$CO formation theories during planet formation. We find disk-averaged rotational temperatures and column densities of $33\pm2$ K, ($1.1\pm0.1)\times10^{12}$ cm$^{-2}$ and $25\pm2$ K, $(4.4\pm0.3)\times10^{11}$ cm$^{-2}$ for ortho- and para-H$_2$CO, respectively, and an OPR of $2.49\pm0.23$. A radially resolved analysis shows that the observed H$_2$CO emits mostly at rotational temperatures of 30-40 K, corresponding to a layer with $z/R\ge0.25$. The OPR is consistent with 3 within 60 au, the extent of the pebble disk, and decreases beyond 60 au to $2.0\pm0.5$. The latter corresponds to a spin temperature of 12 K, well below the rotational temperature. The combination of relatively uniform emitting conditions, a radial gradient in the OPR, and recent laboratory experiments and theory on OPR ratios after sublimation, lead us to speculate that gas-phase formation is responsible for the observed H$_2$CO across the TW Hya disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heavy Ion Induced SEU Sensitivity Evaluation of 3D Integrated SRAMs

Heavy ions induced single event upset (SEU) sensitivity of three-dimensional integrated SRAMs are evaluated by using Monte Carlo sumulation methods based on Geant4. The cross sections of SEUs and Multi Cell Upsets (MCUs) for 3D SRAM are simulated by using heavy ions with different energies and LETs. The results show that the sensitivity of different die of 3D SRAM has obvious discrepancies at low LET. Average percentage of MCUs of 3D SRAMs rises from 17.2% to 32.95% when LET increases from 42.19 MeV cm2/mg to 58.57MeV cm2/mg. As for a certain LET, the percentage of MCUs shows a notable distinction between face-to-face structure and back-to-face structure. For back-to-face structure, the percentage of MCUs increases with the deeper die. However, the face-to-face die presents the relatively low percentage of MCUs. The comparison of SEU cross sections for planar SRAMs and experiment data are conducted to indicate the effectiveness of our simulation method. Finally, we compare the upset cross sections of planar process and 3D integrated SRAMs. Results demonstrate that the sensitivity of 3D SRAMs is not more than that of planar SRAMs and the 3D structure can be become a great potential application for aerospace and military domain.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Resolved Gas Interior to the Dust Rings of the HD 141569 Disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Radial Distribution of H2 and CO in TW Hya as Revealed by Resolved ALMA Observations of CO Isotopologues

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Ring of C2H in the Molecular Disk Orbiting TW Hya

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analytical Formulas of Molecular Ion Abundances and N2H+ Ring in Protoplanetary Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chemical Imaging of the CO Snow Line in the HD 163296 Disk

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Double DCO+ rings reveal CO ice desorption in the outer disk around IM Lup

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for DCO+ as a probe of ionization in the warm disk surface

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HNC in Protoplanetary 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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Near-IR Polarized Scattered Light Imagery of the DoAr 28 Transitional Disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Next Generation Very Large Array Memo No. 6, Science Working Group 1: The Cradle of Life

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Steepening of the 820 micron continuum surface-brightness profile signals dust evolution in TW Hya's disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The cometary composition of a protoplanetary disk as revealed by complex cyanides

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the X-ray and Cosmic Ray Ionization Chemistry of the TW Hya Protoplanetary Disk: Evidence for a Sub-interstellar Cosmic Ray Rate

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, $ζ_{\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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Peculiar Near-Nucleus Outgassing of Comet 17P/Holmes During Its 2007 Outburst

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Physical and chemical structure of planet-forming disks probed by millimeter observations and modeling

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Significantly Low CO Abundance Toward the TW Hya Protoplanetary Disk: A Path to Active Carbon Chemistry?

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Spatially Resolved Vertical Temperature Gradient in the HD 163296 Disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First detection of c-C3H2 in a circumstellar disk

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

H2CO and N2H+ in Protoplanetary Disks: Evidence for a CO-ice Regulated Chemistry

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imaging of the CO Snow Line in a Solar Nebula Analog

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence for Multiple Pathways to Deuterium Enhancements in Protoplanetary Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Kinematics of the CO Gas in the Inner Regions of the TW Hya Disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Disk Imaging Survey of Chemistry with SMA: II. Southern Sky Protoplanetary Disk Data and Full Sample Statistics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Millimeter Observations of the Transition Disk around HD135344B (SAO 206462)

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolved Images of Large Cavities in Protoplanetary Transition Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolved Submillimeter Observations of the HR 8799 and HD 107146 Debris Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolving the CO Snow Line in the Disk around HD 163296

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Ionization Fraction in the DM Tau Protoplanetary Disk

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The TW Hya Disk at 870 microns: Comparison of CO and Dust Radial Structures

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Empirical Constraints on Turbulence in Protoplanetary Accretion Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Imaging the Molecular Disk Orbiting the Twin Young Suns of V4046 Sgr

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protoplanetary Disk Structures in Ophiuchus II: Extension to Fainter Sources

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A break in the gas and dust surface density of the disc around the T Tauri star IM Lup

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Spatially Resolved Inner Hole in the Disk around GM Aurigae

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

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Interferometric Observations of Formaldehyde in the Protoplanetary Disk around LkCa15

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

0 citations0 reviewsNo code linkedOpen access