Source author record

Karin I. Oberg

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

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

preprint2022arXiv

First images of phosphorus molecules towards a proto-Solar analog

The chemistry of phosphorus in star- and planet-forming regions is poorly understood, despite the central role of phosphorus in terrestrial biochemistry. We present ALMA Band 3 and 4 observations of PO and PN towards the Class I protostar B1-a, representing the first spatially resolved observations of phosphorus carriers towards a Solar-type star forming region. The phosphorus molecules emit from two distinct clumps, which coincide with regions where the protostellar outflow (traced by SiO) interacts with a filament of dense gas (traced by CCS). Thus, the gas-phase phosphorus seems to originate from the shocking of dense interstellar clumps. Based on the observed emission patterns, PO and PN appear to be daughter products of a solid phosphorus carrier with an intermediate volatility between ices and silicate grains. Interstellar shocks may therefore play an important role in converting semi-refractory phosphorus to a more volatile form prior to incorporation into cometary ices. Indeed, the (PO+PN)/CH3OH ratio is similar in B1-a and comet 67P, implying a comparable reservoir of volatile phosphorus. The PO/PN ratio ranges from ~1-8 across B1-a. The northern emission clump exhibits a lower PO/PN ratio and weaker 13CH3OH emission than southern clump, indicating distinct shock physics and chemistry at the two positions. Resolved observations of P carriers towards additional sources are needed to better understand what regulates such variations in the PO/PN ratio in protostellar environments.

preprint2022arXiv

HCN snowlines in protoplanetary disks: constraints from ice desorption experiments

HCN is among the most commonly detected molecules in star- and planet-forming regions. It is of broad interest as a tracer of star-formation physics, a probe of nitrogen astrochemistry, and an ingredient in prebiotic chemical schemes. Despite this, one of the most fundamental astrochemical properties of HCN remains poorly characterized: its thermal desorption behavior. Here, we present a series of experiments to characterize the thermal desorption of HCN in astrophysically relevant conditions, with a focus on predicting the HCN sublimation fronts in protoplanetary disks. We derive HCN-HCN and HCN-H2O binding energies of 3207\pm197 K and 4192\pm68 K, which translate to disk midplane sublimation temperatures around 85 K and 103 K. For a typical midplane temperature profile, HCN should only begin to sublimate ~1-2 au exterior to the H2O snow line. Additionally, in H2O-dominated mixtures (20:1 H2O:HCN), we find that the majority of HCN remains trapped in the ice until H2O crystallizes. Thus, HCN may be retained in disk ices at almost all radii where H2O-rich planetesimals form. This implies that icy body impacts to planetary surfaces should commonly deliver this potential prebiotic ingredient. A remaining unknown is the extent to which HCN is pure or mixed with H2O in astrophysical ices, which impacts the HCN desorption behavior as well as the outcomes of ice-phase chemistry. Pure HCN and HCN:H2O mixtures exhibit distinct IR bands, raising the possibility that the James Webb Space Telescope will elucidate the mixing environment of HCN in star- and planet-forming regions and address these open questions.

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.

preprint2021arXiv

Hot corino chemistry in the Class I binary source Ser-emb 11

We report the detection of more than 120 emission lines corresponding to 8 complex organic molecules (CH3OH, CH3CH2OH, CH3OCH3, CH3OCHO, CH3COCH3, NH2CHO, CH2DCN, and CH3CH2CN) and 3 isotopologues (CH2DOH, 13CH3CN, and CH3C15N) toward the western component of the Ser-emb 11 binary young stellar object (YSO) using observations with the Atacama Large Millimeter/submillimeter Array at ~1 mm. The complex organic emission was unresolved with a ~0.5" beam (~220 au) in a compact region around the central protostar, and a population diagram analysis revealed excitation temperatures above 100 K for all COMs, indicating the presence of a hot corino. The estimated column densities were in the range of 10^17 - 10^18 cm^-2 for the O-bearing COMs, and three orders of magnitude lower for the N-bearing species. We also report the detection of H2CO and CH3OH emission in a nearby millimeter source that had not been previously catalogued. Ser-emb 11 is classified in the literature as a Class I source near the Class 0/I cutoff. The estimated COM relative abundances in Ser-emb 11 W and the other three Class I hot corino sources reported in the literature are consistent with those of Class 0 hot corinos, suggesting a continuity in the chemical composition of hot corinos during protostellar evolution.

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

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.

preprint2020arXiv

Astrochemistry and compositions of planetary systems

Planets form and obtain their compositions in disks of gas and dust around young stars. The chemical compositions of these planet-forming disks regulate all aspects of planetary compositions from bulk elemental inventories to access to water and reactive organics, i.e. a planet's hospitality to life and its chemical origins. Disk chemical structures are in their turn governed by a combination of {\it in situ} chemical processes, and inheritance of molecules from the preceding evolutionary stages of the star formation process. In this review we present our current understanding of the chemical processes active in pre- and protostellar environments that set the initial conditions for disks, and the disk chemical processes that evolve the chemical conditions during the first million years of planet formation. We review recent observational, laboratory and theoretical discoveries that have led to the present view of the chemical environment within which planets form, and their effects on the compositions of nascent planetary systems. We also discuss the many unknowns that remain and outline some possible pathways to addressing them.

preprint2020arXiv

Formation of NH2CHO and CH3CHO upon UV photoprocessing of interstellar ice analogs

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

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.

preprint2016arXiv

Carbon Chains and Methanol toward Embedded Protostars

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

preprint2016arXiv

Evidence for a CO desorption front in the outer AS 209 disk

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

preprint2016arXiv

Excess C/O and C/H in outer protoplanetary disk gas

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

preprint2016arXiv

First detection of gas-phase methanol in a protoplanetary disk

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

preprint2016arXiv

On the Inference of the Cosmic Ray Ionization Rate $ζ$ from the HCO$^+$-to-DCO$^+$ Abundance Ratio: The Effect of Nuclear Spin

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

preprint2016arXiv

Photochemistry and astrochemistry: photochemical pathways to interstellar complex organic molecules

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

preprint2016arXiv

Ringed Substructure and a Gap at 1 AU in the Nearest Protoplanetary Disk

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

preprint2016arXiv

The Role of Ice Compositions for Snowlines and the C/N/O Ratios in Active Disks

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

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.

preprint2015arXiv

CO diffusion into amorphous H2O ices

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

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.

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.

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.

preprint2014arXiv

Complex organic molecules during low-mass star formation: Pilot survey results

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

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.

preprint2014arXiv

The HNC/HCN Ratio in Star-Forming Regions

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

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.

preprint2013arXiv

The Spatial Distribution of Organics toward the High-Mass YSO NGC 7538 IRS9

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

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.

preprint2011arXiv

CO ice photodesorption: A wavelength-dependent study

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

preprint2011arXiv

Complex molecules toward low-mass protostars: the Serpens core

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

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.

preprint2011arXiv

Ices in starless and starforming cores

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

preprint2011arXiv

Laboratory H2O:CO2 ice desorption data: entrapment dependencies and its parameterization with an extended three-phase model

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

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.

preprint2011arXiv

The effects of snowlines on C/O in planetary atmospheres

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

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.

preprint2011arXiv

The Spitzer ice legacy: Ice evolution from cores to protostars

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

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.

preprint2010arXiv

A cold complex chemistry toward the low-mass protostar B1-b: evidence for complex molecule production in ices

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

preprint2010arXiv

The c2d Spitzer Spectroscopic Survey of Ices Around Low-Mass Young Stellar Objects. IV. NH3 and CH3OH

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

preprint2010arXiv

The effect of H2O on ice photochemistry

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

preprint2009arXiv

Cold gas as an ice diagnostic toward low mass protostars

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

preprint2009arXiv

Photodesorption of ices I: CO, N2 and CO2

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

preprint2008arXiv

Photodesorption of Ices II: H2O and D2O

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

preprint2008arXiv

The c2d Spitzer spectroscopy survey of ices around low-mass young stellar objects, III: CH4

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

preprint2006arXiv

Effects of CO2 on H2O band profiles and band strengths in mixed H2O:CO2 ices

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