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

L. Ilsedore Cleeves

L. Ilsedore Cleeves 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.GA

Catalog footprint

What is connected

33works
3topics
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

33 published item(s)

preprint2026arXiv

Protoplanetary disk cavities with JWST-MIRI: a dichotomy in molecular emission

The evolution of planet-forming regions in protoplanetary disks is of fundamental importance to understanding planet formation. Disks with a central deficit in dust emission, a "cavity", have long attracted interest as potential evidence for advanced disk clearing by protoplanets and/or winds. Before JWST, infrared spectra showed that these disks typically lack the strong molecular emission observed in full disks. In this work, we combine a sample of 12 disks with millimeter cavities of a range of sizes ($\sim2$-70 au) and different levels of millimeter and infrared continuum deficits. We analyze their molecular spectra as observed with MIRI on JWST, homogeneously reduced with the new JDISCS pipeline. This analysis demonstrates a stark dichotomy in molecular emission where "molecule-rich" (MR) cavities follow global trends between water, CO, and OH luminosity and accretion luminosity as in full disks, while "molecule-poor" (MP) cavities are significantly sub-luminous in all molecules except sometimes OH. Disk cavities generally show sub-luminous organic emission, higher OH/H$_2$O ratios, and suggest a lower water column density. The sub-thermal excitation of CO and water vibrational lines suggests a decreased gas density in the emitting layer in all cavities, supporting model expectations for C$_2$H$_2$ photodissociation. We discover a bifurcation in infrared index (lower in MR cavities) suggesting that the molecular dichotomy is linked to residual $μ$m-size dust within millimeter disk cavities. Put together, these results suggest a feedback process between dust depletion, gas density decrease, and molecule dissociation. Disk cavities may have a common evolutionary sequence where MR switch into MP over time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Classification of X-ray Flare Driven Chemical Variability in Protoplanetary Disks

Young stars are highly variable in the X-ray regime. In particular, bright X-ray flares can substantially enhance ionization in the surrounding protoplanetary disk. Since disk chemical evolution is impacted by ionization, X-ray flares have the potential to fundamentally alter the chemistry of planet forming regions. We present 2D disk chemical models that incorporate a stochastic X-ray flaring module, named \xgen, and examine the flares' overall chemical impact compared to models that assume a constant X-ray flux. We examine the impact of 500 years of flaring events and find global chemical changes on both short time scales (days) in response to discrete flaring events and long time-scales (centuries) in response to the cumulative impact of many flares. Individual X-ray flares most strongly affect small gas-phase cations, where a single flare can temporarily enhance the abundance of species such as H$_3^+$, HCO$^+$, CH$_3^+$, and C$^+$. We find that flares can also drive chemistry out of "steady state" over longer time periods, where the disk-integrated abundance of some species, such as O and O$_2$, changes by a few percent over the 500 year model. We also explore whether the specific history of X-ray flaring events (randomly drawn but from the same energy distribution) impacts the chemical evolution and find that it does not. Finally, we examine the impact of X-ray flares on the electron fraction. While most molecules modeled are not highly sensitive to flares, certain species, including observable molecules, are very reactive to the dynamic environment of a young star.

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
preprint2022arXiv

New Constraints on Protoplanetary Disk Gas Masses in Lupus

Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H$_2$ content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologues of CO. However, without additional constraints, the degeneracy between H$_2$ mass and the elemental composition of the gas leads to large uncertainties in such estimates. Here we explore the gas compositions of seven disks from the Lupus sample representing a range of CO-to-dust ratios. With Band 6 and 7 ALMA observations, we measure line emission for HCO$^+$, HCN, and N$_2$H$^+$. We find a tentative correlation among the line fluxes for these three molecular species across the sample, but no correlation with $^{13}$CO or sub-mm continuum fluxes. For the three disks where N$_2$H$^+$ is detected, we find that a combination of high disk gas masses and sub-interstellar C/H and O/H are needed to reproduce the observed values. We find increases of $\sim$10-100$\times$ previous mass estimates are required to match the observed line fluxes. This study highlights how multi-molecular studies are essential for constraining the physical and chemical properties of the gas in populations of protoplanetary disks and that CO isotopologues alone are not sufficient for determining the mass of many observed disks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The ODYSSEUS Survey. Motivation and First Results: Accretion, Ejection, and Disk Irradiation of CVSO 109

The Hubble UV Legacy Library of Young Stars as Essential Standards (ULLYSES) Director's Discretionary Program of low-mass pre-main-sequence stars, coupled with forthcoming data from ALMA and JWST, will provide the foundation to revolutionize our understanding of the relationship between young stars and their protoplanetary disks. A comprehensive evaluation of the physics of disk evolution and planet formation requires understanding the intricate relationships between mass accretion, mass outflow, and disk structure. Here we describe the Outflows and Disks around Young Stars: Synergies for the Exploration of ULLYSES Spectra (ODYSSEUS) Survey and present initial results of the classical T Tauri Star CVSO 109 in Orion OB1b as a demonstration of the science that will result from the survey. ODYSSEUS will analyze the ULLYSES spectral database, ensuring a uniform and systematic approach in order to (1) measure how the accretion flow depends on the accretion rate and magnetic structures, (2) determine where winds and jets are launched and how mass-loss rates compare with accretion, and (3) establish the influence of FUV radiation on the chemistry of the warm inner regions of planet-forming disks. ODYSSEUS will also acquire and provide contemporaneous observations at X-ray, optical, NIR, and millimeter wavelengths to enhance the impact of the ULLYSES data. Our goal is to provide a consistent framework to accurately measure the level and evolution of mass accretion in protoplanetary disks, the properties and magnitudes of inner-disk mass loss, and the influence of UV radiation fields that determine ionization levels and drive disk chemistry.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamical Masses and Stellar Evolutionary Model Predictions of M-Stars

In this era of Gaia and ALMA, dynamical stellar mass measurements provide benchmarks that are independent of observations of stellar characteristics and their uncertainties. These benchmarks can then be used to validate and improve stellar evolutionary models, which can lead to both imprecise and inaccurate mass predictions for pre-main-sequence, low-mass stars. We present the dynamical stellar masses derived from disks around three M-stars (FP Tau, J0432+1827, and J1100-7619) using ALMA observations of $^{12}$CO (J=2--1) and $^{13}$CO (J=2--1) emission. These are the first dynamical stellar mass measurements for J0432+1827 and J1100-7619 and the most precise measurement for FP Tau. Fiducial stellar evolutionary model tracks, which do not include any treatment of magnetic activity, agree with the dynamical measurement of J0432+1827 but underpredict the mass by $\sim$60\% for FP Tau and $\sim$80\% for J1100-7619. Possible explanations for the underpredictions include inaccurate assumptions of stellar effective temperature, undetected binarity for J1100-7619, and that fiducial stellar evolutionary models are not complex enough to represent these stars. In the former case, the stellar effective temperatures would need to be increased by $\sim$40K to $\sim$340K to reconcile the fiducial model predictions with the dynamically-measured masses. In the latter case, we show that the dynamical masses can be reproduced using results from stellar evolutionary models with starspots, which incorporate fractional starspot coverage to represent the manifestation of magnetic activity. Folding in low-mass M-stars from the literature and assuming that the stellar effective temperatures are imprecise but accurate, we find tentative evidence of a relationship between fractional starspot coverage and observed effective temperature for these young, cool stars.

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

Chemistry Along Accretion Streams in a Viscously-Evolving Protoplanetary Disk

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

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
preprint2017arXiv

Composition of Early Planetary Atmospheres II: Coupled Dust and Chemical Evolution in Protoplanetary Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Exploring the Origins of Deuterium Enrichments in Solar Nebular Organics

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

First detection of gas-phase ammonia in a planet-forming disk

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hydrocarbon emission rings in protoplanetary disks induced by dust evolution

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Commonality of 10-30AU Sized Axisymmetric Dust Structures in Protoplanetary Disks

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

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

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

Indirect Detection of Forming Protoplanets via Chemical Asymmetries in Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiple Carbon monoxide Snow-lines in Disks Sculpted by Radial Drift

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Distribution and Chemistry of H$_2$CO in the DM Tau Protoplanetary Disk

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

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

Exploring the Origins of Carbon in Terrestrial Worlds

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The ancient heritage of water ice in the solar system

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wind-driven Exclusion of Cosmic Rays in the Protoplanetary Disk Environment

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

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

Exclusion of Cosmic Rays in Protoplanetary Disks: Stellar and Magnetic Effects

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Radionuclide Ionization in Protoplanetary Disks: Calculations of Decay Product Radiative Transfer

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Detection of the Water Reservoir in a Forming Planetary System

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transition disk chemistry and future prospects with ALMA

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

0 citations0 reviewsNo code linkedOpen access