Source author record

Ewine F. van Dishoeck

Ewine F. van Dishoeck 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

124works
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

124 published item(s)

preprint2026arXiv

Circumplanetary Disk Candidate in the Disk of HD 163296 Traced by Localized Emission from Simple Organics

Atacama Large Millimeter/submillimeter Array observations suggest that the disc of HD 163296 is being actively shaped by embedded, yet unseen protoplanets, as indicated by numerous gas and dust substructures consistent with planet-disc interaction models. We report the first detection of simple organic molecules, HCN and C2H, tracing a candidate circumplanetary disc (CPD) in the HD 163296 system, located at an orbital radius of $R=88\pm7$ au and azimuth $ϕ=46\pm3^\circ$ (or $R=0.75''$, $\rm{PA}=350^\circ$ in projected sky coordinates), and originating near the midplane of the circumstellar disc. The signature is localised but spectrally resolved, and it overlaps with a previously reported planet candidate, P94, identified through kinematic perturbations traced by CO lines. We propose a scenario in which the observed chemical anomalies arise from increased heating driven by the forming planet and ongoing accretion through its CPD, facilitating the thermal desorption of species that would otherwise remain frozen out in the disc midplane, and potentially triggering the activation barriers of chemical reactions that lead to enhanced molecular production. Based on a first-order dynamical analysis of the HCN spectrum from the CPD--isolated with a 7$σ$ significance--we infer an upper limit on the planet mass of 1.8 $M_{\rm Jup}$, consistent with predictions from CO kinematics and constraints from direct imaging studies. By comparing the CPD sizes derived from our models with theoretical expectations where the CPD radius corresponds to roughly one-third of the planet's Hill radius, we favor CPD gas temperatures $T > 150$ K, planet masses $M_{\rm p} < 1.0$ $M_{\rm Jup}$, and CPD radii $R_{\rm CPD} < 2$ au.

preprint2025arXiv

MINDS. Young binary systems with JWST/MIRI: Variable water-rich primaries and extended emission

As part of the JWST GTO program MINDS, we analyze the mid-infrared emission of three Class II binary systems: VW Cha, WX Cha, and RW Aur, to investigate the impact of stellar multiplicity on the chemistry and physics of their inner disk. We analyze the 1D spectrum from JWST/MIRI-MRS for primary and secondary disks separately, extracted by combining forward modeling with a theoretical PSF and aperture photometry. We modeled the molecular lines with 0D slab models. We interpret the results by comparing our JWST spectra to VLT/CRIRES+, Spitzer/IRS, and ALMA. Primary and secondary disks are dramatically different in their mid-infrared emission, with primary disks showing H2O-rich spectra, and secondary disks being mostly line poor to the sensitivity of our spectra. When comparing MIRI-MRS to Spitzer/IRS, we observe large variability in the line emission of VW Cha A, as well as in the continuum of RW Aur A. The disks around VW Cha BC and RW Aur B show evidence of ionizing radiation, and a further comparison with ALMA at high angular resolution dust continuum suggest that the spectrum of RW Aur B is well explained by its ~4 au cavity. All the systems show [Ne II] jet emission, and three of them also show spatially resolved emission structures in H2, likely originated by outflows and dynamical interactions. Many of the observed features in the primary disks, such as enhanced water emission, could be linked to the increased accretion and radial drift produced by dynamical disk truncation. However, additional mechanisms are needed to explain the large differences between primary and secondary disks, potentially inner disk substructures. This work is an example of the need for combining multiple facilities to fully understand the observations from JWST.

preprint2022arXiv

A major asymmetric ice trap in a planet-forming disk: III. First detection of dimethyl ether

The complex organic molecules (COMs) detected in star-forming regions are the precursors of the prebiotic molecules that can lead to the emergence of life. By studying COMs in more evolved protoplanetary disks we can gain a better understanding of how they are incorporated into planets. This paper presents ALMA band 7 observations of the dust and ice trap in the protoplanetary disk around Oph IRS 48. We report the first detection of dimethyl ether (CH3OCH3) in a planet-forming disk and a tentative detection of methyl formate (CH3OCHO). We determined column densities for the detected molecules and upper limits on non-detected species using the CASSIS spectral analysis tool. The inferred column densities of CH3OCH3 and CH3OCHO with respect to methanol (CH3OH) are of order unity, indicating unusually high abundances of these species compared to other environments. Alternatively, the 12CH3OH emission is optically thick and beam diluted, implying a higher CH3OH column density and a smaller emitting area than originally thought. The presence of these complex molecules can be explained by thermal ice sublimation, where the dust cavity edge is heated by irradiation and the full volatile ice content is observable in the gas phase. This work confirms the presence of oxygen-bearing molecules more complex than CH3OH in protoplanetary disks for the first time. It also shows that it is indeed possible to trace the full interstellar journey of COMs across the different evolutionary stages of star, disk, and planet formation.

preprint2022arXiv

A new planet candidate detected in a dust gap of the disc around HD 163296 through localised kinematic signatures. An observational validation of the Discminer

We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at $R=94\pm6$ au, $ϕ=50\pm3^\circ$ (P94), is near the centre of one of the gaps in dust continuum emission, and is consistent with a planet mass of 1 $\rm M_{\rm Jup}$. The other planet, located at $R=261\pm4$ au, $ϕ=57\pm1^\circ$ (P261), is in the region where a velocity kink was previously observed in $^{12}$CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disc, and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles we detect gas gaps at $R=38$ au, $R=88$ au, and $R=136$ au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at $R=88$ au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the DISCMINER is capable of finding localised, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary discs.

preprint2022arXiv

Dependence of $X_{\rm CO}$ on metallicity, intensity, and spatial scale in a self-regulated interstellar medium

We study the CO(1-0)-to-H$_2$ conversion factor ($X_{\rm CO}$) and the line ratio of CO(2-1)-to-CO(1-0) ($R_{21}$) across a wide range of metallicity ($0.1 \leq Z/Z_\odot \leq 3$) in high-resolution (~0.2 pc) hydrodynamical simulations of a self-regulated multiphase interstellar medium. We construct synthetic CO emission maps via radiative transfer and systematically vary the "observational" beam size to quantify the scale dependence. We find that the kpc-scale $X_{\rm CO}$ can be over-estimated at low $Z$ if assuming steady-state chemistry or assuming that the star-forming gas is H$_2$-dominated. On parsec scales, $X_{\rm CO}$ varies by orders of magnitude from place to place, primarily driven by the transition from atomic carbon to CO. The pc-scale $X_{\rm CO}$ drops to the Milky Way value of $2\times 10^{20}\ {\rm cm^{-2}~(K~km~s^{-1})^{-1}}$ once dust shielding becomes effective, independent of $Z$. The CO lines become increasingly optically thin at lower $Z$, leading to a higher $R_{21}$. Most cloud area is filled by diffuse gas with high $X_{\rm CO}$ and low $R_{21}$, while most CO emission originates from dense gas with low $X_{\rm CO}$ and high $R_{21}$. Adopting a constant $X_{\rm CO}$ strongly over- (under-)estimates H$_2$ in dense (diffuse) gas. The line intensity negatively (positively) correlates with $X_{\rm CO}$ ($R_{21}$) as it is a proxy of column density (volume density). On large scales, $X_{\rm CO}$ and $R_{21}$ are dictated by beam averaging, and they are naturally biased towards values in dense gas. Our predicted $X_{\rm CO}$ is a multivariate function of $Z$, line intensity, and beam size, which can be used to more accurately infer the H$_2$ mass.

preprint2022arXiv

Sulphur monoxide emission tracing an embedded planet in the HD 100546 protoplanetary disk

Molecular line observations are powerful tracers of the physical and chemical conditions across the different evolutionary stages of star, disk and planet formation. Using the high angular resolution and unprecedented sensitivity of the Atacama Large Millimeter Array (ALMA) there is now a drive to detect small-scale gas structures in protoplanetary disks that can be attributed directly to forming planets. We report high angular resolution ALMA Band 7 observations of sulphur monoxide (SO) in the nearby planet-hosting disk around Herbig star HD 100546. SO is rarely detected in evolved protoplanetary disks but in other environments, it is most often utilised as a tracer of shocks. The SO emission from the HD 100546 disk is primarily originating from gas within the approx. 20 au mm-dust cavity and shows a clear azimuthal brightness asymmetry of a factor of 2. In addition, we see a significant difference in the line profile shape when comparing these new Cycle 7 data to Cycle 0 data of the same SO transitions. We discuss the different physical/chemical mechanisms that could be responsible for this asymmetry and time variability including disk winds, disk warps, and a shock triggered by a (forming) planet. We propose that the SO is enhanced in the cavity due to the presence of a giant planet. The SO asymmetry complements evidence for hot circumplanetary material around the giant planet HD 100546 c traced via CO ro-vibrational emission. This work sets the stage for further observational and modelling efforts to detect and understand the chemical imprint of a forming planet on its parent disk.

preprint2021arXiv

A cold accretion flow onto one component of a multiple protostellar system

Context: Gas accretion flows transport material from the cloud core onto the protostar. In multiple protostellar systems, it is not clear if the delivery mechanism is preferential or evenly distributed among the components. Aims: Gas accretion flows within IRAS16293 is explored out to 6000 AU. Methods: ALMA Band 3 observations of low-$J$ transitions of HNC, cyanopolyynes (HC$_3$N, HC$_5$N), and N$_2$H$^+$ are used to probe the cloud core structure at ~100 AU resolution. Additional Band 3 archival data provide low-$J$ HCN and SiO lines. These data are compared with the corresponding higher-$J$ lines from the PILS Band 7 data for excitation analysis. The HNC/HCN ratio is used as a temperature tracer. Results: The low-$J$ transitions of HC$_3$N, HC$_5$N, HNC and N$_2$H$^+$ trace extended and elongated structures from 6000 AU down to ~100 AU, without accompanying dust continuum emission. Two structures are identified: one traces a flow that is likely accreting toward the most luminous component of the system IRAS16293 A. Temperatures inferred from the HCN/HNC ratio suggest that the gas in this flow is cold, between 10 and 30 K. The other structure is part of an UV-irradiated cavity wall entrained by one of the outflows. The two outflows driven by IRAS16293 A present different molecular gas distributions. Conclusions: Accretion of cold gas is seen from 6000 AU scales onto IRAS16293 A, but not onto source B, indicates that cloud core material accretion is competitive due to feedback onto a dominant component in an embedded multiple protostellar system. The preferential delivery of material could explain the higher luminosity and multiplicity of source A compared to source B. The results of this work demonstrate that several different molecular species, and multiple transitions of each species, are needed to confirm and characterize accretion flows in protostellar cloud cores.

preprint2021arXiv

Hydrogenation of accreting C-atoms and CO molecules -- simulating ketene and acetaldehyde formation under dark and translucent cloud conditions

Simple and complex organic molecules (COMs) are observed along different phases of star and planet formation and have been successfully identified in prestellar environments such as dark and translucent clouds. Yet the picture of organic molecule formation at those earliest stages of star formation is not complete and an important reason is the lack of specific laboratory experiments that simulate carbon atom addition reactions on icy surfaces of interstellar grains. Here we present experiments in which CO molecules as well as C- and H-atoms are co-deposited with H$_2$O molecules on a 10 K surface mimicking the ongoing formation of an "H$_2$O-rich" ice mantle. To simulate the effect of impacting C-atoms and resulting surface reactions with ice components, a specialized C-atom beam source is used, implemented on SURFRESIDE$^3$, an UHV cryogenic setup. Formation of ketene (CH$_2$CO) in the solid state is observed "in situ" by means of reflection absorption IR spectroscopy. C$^1$$^8$O and D isotope labelled experiments are performed to further validate the formation of ketene. Data analysis supports that CH$_2$CO is formed through C-atom addition to a CO-molecule, followed by successive hydrogenation transferring the formed :CCO into ketene. Efficient formation of ketene is in line with the absence of an activation barrier in C+CO reaction reported in the literature. We also discuss and provide experimental evidence for the formation of acetaldehyde (CH$_3$CHO) and possible formation of ethanol (CH$_3$CH$_2$OH), two COM derivatives of CH$_2$CO hydrogenation. The underlying reaction network is presented and the astrochemical implications of the derived pathways are discussed.

preprint2021arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound.

preprint2020arXiv

Dust masses of young disks: constraining the initial solid reservoir for planet formation

In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $κ_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $κ_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks.

preprint2020arXiv

Protoplanetary disk masses in NGC 2024: Evidence for two populations

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.

preprint2020arXiv

Solving grain size inconsistency between ALMA polarization and VLA continuum in the Ophiuchus IRS 48 protoplanetary disk

The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $μ$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $μ$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $μ$m wavelength is optically thick ($τ_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $μ$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $μ$m thermal emission is optically thick.

preprint2020arXiv

Temperature structures of embedded disks: young disks in Taurus are warm

The chemical composition of gas and ice in disks around young stars set the bulk composition of planets. In contrast to protoplanetary disks (Class II), young disks that are still embedded in their natal envelope (Class 0 and I) are predicted to be too warm for CO to freeze out, as has been confirmed observationally for L1527 IRS. To establish whether young disks are generally warmer than their more evolved counterparts, we observed five young (Class 0/I and Class I) disks in Taurus with the Atacama Large Millimeter/submillimeter Array (ALMA), targeting C$^{17}$O $2-1$, H$_2$CO $3_{1,2}-2_{1,1}$, HDO $3_{1,2}-2_{2,1}$ and CH$_3$OH $5_K-4_K$ transitions at $0.48^{\prime\prime} \times 0.31^{\prime\prime}$ resolution. The different freeze-out temperatures of these species allow us to derive a global temperature structure. C$^{17}$O and H$_2$CO are detected in all disks, with no signs of CO freeze-out in the inner $\sim$100 au, and a CO abundance close to $\sim$10$^{-4}$. H$_2$CO emission originates in the surface layers of the two edge-on disks, as witnessed by the especially beautiful V-shaped emission pattern in IRAS~04302+2247. HDO and CH$_3$OH are not detected, with column density upper limits more than 100 times lower than for hot cores. Young disks are thus found to be warmer than more evolved protoplanetary disks around solar analogues, with no CO freeze-out (or only in the outermost part of $\gtrsim$100 au disks) or CO processing. However, they are not as warm as hot cores or disks around outbursting sources, and therefore do not have a large gas-phase reservoir of complex molecules.

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.

preprint2020arXiv

The TW Hya Rosetta Stone Project III: Resolving the Gaseous Thermal Profile of the Disk

The thermal structure of protoplanetary disks is a fundamental characteristic of the system that has wide reaching effects on disk evolution and planet formation. In this study, we constrain the 2D thermal structure of the protoplanetary disk TW Hya structure utilizing images of seven CO lines. This includes new ALMA observations of 12CO J=2-1 and C18O J=2-1 as well as archival ALMA observations of 12CO J=3-2, 13CO J=3-2, 6-5, C18O J= 3-2, 6-5. Additionally, we reproduce a Herschel observation of the HD J=1-0 line flux, the spectral energy distribution, and utilize a recent quantification of CO radial depletion in TW Hya. These observations were modeled using the thermochemical code RAC2D, and our best fit model reproduces all spatially resolved CO surface brightness profiles. The resulting thermal profile finds a disk mass of 0.025 Msun and a thin upper layer of gas depleted of small dust with a thickness of approx 1.2% of the corresponding radius. Using our final thermal structure, we find that CO alone is not a viable mass tracer as its abundance is degenerate with the total H2 surface density. Different mass models can readily match the spatially resolved CO line profiles with disparate abundance assumptions. Mass determination requires additional knowledge and, in this work, HD provides the additional constraint to derive the gas mass and supports the inference of CO depletion in the TW Hya disk. Our final thermal structure confirms the use of HD as a powerful probe of protoplanetary disk mass. Additionally, the method laid out in this paper is an employable strategy for extraction of disk temperatures and masses in the future.

preprint2020arXiv

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

preprint2020arXiv

Vertical gas accretion impacts the carbon-to-oxygen ratio of gas giant atmospheres

Recent theoretical, numerical, and observational work have suggested that when a growing planet opens a gap in its disk the flow of gas into the gap is dominated by gas falling vertically from a height of at least one gas scale height. Our primary objective is to include, for the first time, the chemical impact that accreting gas above the midplane will have on the resulting C/O. We compute the accretion of gas onto planetary cores beginning at different disk radii and track the chemical composition of the gas and small icy grains to predict the resulting carbon-to-oxygen ratio (C/O) in their atmospheres. In our model, all of the planets which began their evolution inward of 60 AU open a gap in the gas disk, and hence are chemically affected by the vertically accreting gas. Two important conclusions follow from this vertical flow: (1) more oxygen rich icy dust grains become available for accretion onto the planetary atmosphere. (2) The chemical composition of the gas dominates the final C/O of planets in the inner ($<$ 20 AU) part of the disk. This implies that with the launch of the James Webb Space Telescope we can trace the disk material that sets the chemical composition of exoplanetary atmospheres.

preprint2016arXiv

A primordial origin for molecular oxygen in comets: A chemical kinetics study of the formation and survival of O$_2$ ice from clouds to disks

Molecular oxygen has been confirmed as the fourth most abundant molecule in cometary material O$_2$/H$_2$O $\sim 4$ %) and is thought to have a primordial nature, i.e., coming from the interstellar cloud from which our solar system was formed. However, interstellar O$_2$ gas is notoriously difficult to detect and has only been observed in one potential precursor of a solar-like system. Here, the chemical and physical origin of O$_2$ in comets is investigated using sophisticated astrochemical models. Three origins are considered: i) in dark clouds, ii) during forming protostellar disks, and iii) during luminosity outbursts in disks. The dark cloud models show that reproduction of the observed abundance of O$_2$ and related species in comet 67P/C-G requires a low H/O ratio facilitated by a high total density ($\geq 10^5$ cm$^{-3}$), and a moderate cosmic ray ionisation rate ($\leq 10^{-16}$ s$^{-1}$) while a temperature of 20 K, slightly higher than the typical temperatures found in dark clouds, also enhances the production of O$_2$. Disk models show that O$_2$ can only be formed in the gas phase in intermediate disk layers, and cannot explain the strong correlation between O$_2$ and H$_2$O in comet 67P/C-G together with the weak correlation between other volatiles and H$_2$O. However, primordial O$_2$ ice can survive transport into the comet-forming regions of disks. Taken together, these models favour a dark cloud (or "primordial") origin for O$_2$ in comets, albeit for dark clouds which are warmer and denser than those usually considered as solar system progenitors.

preprint2016arXiv

Cometary ices in forming protoplanetary disc midplanes

Low-mass protostars are the extrasolar analogues of the natal Solar System. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few % of water ice) at large disc radii (R > 30 AU). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 AU and both are nearly solar in the case of pure infall, but both are sub-solar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically-significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the Solar Nebula.

preprint2016arXiv

Determining protoplanetary disk gas masses from CO isotopologues line observations

Despite intensive studies of protoplanetary disks, there is still no reliable way to determine their total mass and their surface density distribution, quantities that are crucial for describing both the structure and the evolution of disks up to the formation of planets. The goal of this work is to use less abundant CO isotopologues, whose detection is routine for ALMA, to infer the gas mass of disks. Isotope-selective effects need to be taken into account in the analysis, because they can significantly modify CO isotopologues line intensities. CO isotope-selective photodissociation has been implemented in the physical-chemical code DALI and 800 disk models have been run for a range of disk and stellar parameters. Dust and gas temperature structures have been computed self-consistently, together with a chemical calculation of the main species. Both disk structure and stellar parameters have been investigated. Total fluxes have been ray-traced for different CO isotopologues and for various transitions for different inclinations. A combination of 13CO and C18O total intensities allows inference of the total disk mass, although with non-negligible uncertainties. These can be overcome by employing spatially resolved observations, i.e. the disk's radial extent and inclination. Comparison with parametric models shows differences at the factor of a few level, especially for extremely low and high disk masses. Finally, total line intensities for different CO isotopologue and for various low-J transitions are provided and are fitted to simple formulae. The effects of a lower gas-phase carbon abundance and different gas/dust ratios are investigated as well, and comparison with other tracers is made. Disk masses can be determined within a factor of a few by comparing CO isotopologue lines observations with the simulated line fluxes, modulo the uncertainties in the volatile elemental abundances.

preprint2016arXiv

Do siblings always form and evolve simultaneously? Testing the coevality of multiple protostellar systems through SEDs

Multiplicity is common in field stars and among protostellar systems. Models suggest two paths of formation: turbulent fragmentation and protostellar disk fragmentation. We attempt to find whether or not the coevality frequency of multiple protostellar systems can help to better understand their formation mechanism. The coevality frequency is determined by constraining the relative evolutionary stages of the components in a multiple system. SEDs for known multiple protostars in Perseus were constructed from literature data. Herschel PACS photometric maps were used to sample the peak of the SED for systems with separations >7", a crucial aspect in determining the evolutionary stage of a protostellar system. Inclination effects and the surrounding envelope and outflows were considered to decouple source geometry from evolution. This together with the shape and derived properties from the SED was used to determine each system's coevality as accurately as possible. SED models were used to examine the frequency of non-coevality that is due to geometry. We find a non-coevality frequency of 33+/-10% from the comparison of SED shapes of resolved multiple systems. Other source parameters suggest a somewhat lower frequency of non-coevality. The frequency of apparent non-coevality that is due to random inclination angle pairings of model SEDs is 17+/-0.5%. Observations of the outflow of resolved multiple systems do not suggest significant misalignments within multiple systems. Effects of unresolved multiples on the SED shape are also investigated. We find that 1/3 of the multiple protostellar systems sampled here are non-coeval, which is more than expected from random geometric orientations. The other 2/3 are found to be coeval. Higher order multiples show a tendency to be non-coeval. The frequency of non-coevality found here is most likely due to formation and enhanced by dynamical evolution.

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.

preprint2016arXiv

Setting the volatile composition of (exo)planet-building material. Does chemical evolution in disk midplanes matter?

[Abridged] Chemical evolution in the protoplanetary disk midplane can modify the composition of ices and gases. We have investigated if and how chemical evolution affects the abundances and distributions of key volatile species in the midplane of a protoplanetary disk in the 0.2-30 AU range. A full chemical network including gas-phase, gas-grain interactions and grain-surface chemistry is used to evolve chemistry in time, for 1 Myr. Great diversity is observed in the relative abundance ratios of the main considered species: H2O, CO, CO2, CH4, O2, NH3 and N2. The choice of ionisation level, the choice of initial abundances, as well as the extent of chemical reaction types included are all factors that affect the chemical evolution. The only exception is the inheritance scenario with a low ionisation level, which results in negligible changes compared with the initial abundances, regardless of whether grain-surface chemistry is included. The chemical processing changes the C/O ratios for gas and ice significantly with radius and with model assumptions. For high ionisation levels, chemical processing becomes significant after a few times $10^{5}$~yrs. It is concluded that chemical evolution in the disk midplane needs to be considered in the determination of the volatile composition of planetesimals.

preprint2016arXiv

The (w)hole survey: an unbiased sample study of transition disk candidates based on Spitzer catalogs

Understanding disk evolution and dissipation is essential for studies of planet formation. Transition disks, i.e., disks with large dust cavities and gaps, are promising candidates of active evolution. About two dozen SED-selected candidates have been confirmed to have dust cavities through millimeter interferometric imaging, but this sample is biased towards the brightest disks. The Spitzer surveys of nearby low-mass star forming regions have resulted in more than 4000 Young Stellar Objects (YSOs). Using color criteria we have selected a sample of ~150 candidates, and an additional 40 candidates and known transition disks from the literature. The Spitzer data were complemented by new observations at longer wavelengths, including new JCMT and APEX submillimeter photometry, and WISE and Herschel-PACS mid and far-infrared photometry. Furthermore, optical spectroscopy was obtained and stellar types were derived for 85% of the sample, including information from the literature. The SEDs were fit to a grid of RADMC-3D disk models with a limited number of parameters: disk mass, inner disk mass, scale height and flaring, and disk cavity radius, where the latter is the main parameter of interest. A large fraction of the targets possibly have dust cavities based on the SED. The derived cavity sizes are consistent with imaging/modeling results in the literature, where available. Trends are found with Ldisk/Lstar and stellar mass and a possible connection with exoplanet orbital radii. A comparison with a previous study where color observables are used (Cieza et al. 2010) reveals large overlap between their category of planet-forming disks and our transition disks with cavities. A large number of the new transition disk candidates are suitable for follow-up observations with ALMA.

preprint2016arXiv

The CDF Archive: Herschel PACS and SPIRE Spectroscopic Data Pipeline and Products for Protostars and Young Stellar Objects

We present the COPS-DIGIT-FOOSH (CDF) Herschel spectroscopy data product archive, and related ancillary data products, along with data fidelity assessments, and a user-created archive in collaboration with the Herschel-PACS and SPIRE ICC groups. Our products include datacubes, contour maps, automated line fitting results, and best 1-D spectra products for all protostellar and disk sources observed with PACS in RangeScan mode for two observing programs: the DIGIT Open Time Key Program (KPOT_nevans_1 and SDP_nevans_1; PI: N. Evans), and the FOOSH Open Time Program (OT1_jgreen02_2; PI: J. Green). In addition, we provide our best SPIRE-FTS spectroscopic products for the COPS Open Time Program (OT2_jgreen02_6; PI: J. Green) and FOOSH sources. We include details of data processing, descriptions of output products, and tests of their reliability for user applications. We identify the parts of the dataset to be used with caution. The resulting absolute flux calibration has improved in almost all cases. Compared to previous reductions, the resulting rotational temperatures and numbers of CO molecules have changed substantially in some sources. On average, however, the rotational temperatures have not changed substantially (< 2%), but the number of warm (Trot ~ 300 K) CO molecules has increased by about 18%.

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.

preprint2015arXiv

A low-mass protostar's disk-envelope interface: disk-shadowing evidence from ALMA DCO+ observations of VLA1623

Due to instrumental limitations and a lack of disk detections, the structure between the envelope and the rotationally supported disk has been poorly studied. This is now possible with ALMA through observations of CO isotopologs and tracers of freezeout. Class 0 sources are ideal for such studies given their almost intact envelope and young disk. The structure of the disk-envelope interface of the prototypical Class 0 source, VLA1623A which has a confirmed Keplerian disk, is constrained from ALMA observations of DCO+ 3-2 and C18O 2-1. The physical structure of VLA1623 is obtained from the large-scale SED and continuum radiative transfer. An analytic model using a simple network coupled with radial density and temperature profiles is used as input for a 2D line radiative transfer calculation for comparison with the ALMA Cycle 0 12m array and Cycle 2 ACA observations of VLA1623. DCO+ emission shows a clumpy structure bordering VLA1623A's Keplerian disk, suggesting a cold ring-like structure at the disk-envelope interface. The radial position of the observed DCO+ peak is reproduced in our model only if the region's temperature is between 11-16K, lower than expected from models constrained by continuum and SED. Altering the density has little effect on the DCO+ position, but increased density is needed to reproduce the disk traced in C18O. The DCO+ emission around VLA1623A is the product of shadowing of the envelope by the disk. Disk-shadowing causes a drop in the gas temperature outside of the disk on >200AU scales, encouraging deuterated molecule production. This indicates that the physical structure of the disk-envelope interface differs from the rest of the envelope, highlighting the drastic impact that the disk has on the envelope and temperature structure. The results presented here show that DCO+ is an excellent cold temperature tracer.

preprint2015arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: IV. Mechanical and radiative feedback

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

preprint2015arXiv

Complex organic molecules in organic-poor massive young stellar objects

Massive young stellar objects (MYSOs) with hot cores are classic sources of complex organic molecules. The origins of these molecules in such sources, as well as the small- and large-scale differentiation between nitrogen- and oxygen-bearing complex species, are poorly understood. We aim to use complex molecule abundances toward a chemically less explored class of MYSOs with weak hot organic emission lines to constrain the impact of hot molecular cores and initial ice conditions on the chemical composition toward MYSOs. We use the IRAM 30m and the Submillimeter Array to search for complex organic molecules over 8-16 GHz in the 1~mm atmospheric window toward three MYSOs with known ice abundances, but without luminous molecular hot cores. Complex molecules are detected toward all three sources at comparable abundances with respect to CH$_3$OH to classical hot core sources. The relative importance of CH$_3$CHO, CH$_3$CCH, CH$_3$OCH$_3$, CH$_3$CN, and HNCO differ between the organic-poor MYSOs and hot cores, however. Furthermore, the N-bearing molecules are generally concentrated toward the source centers, while most O- and C-bearing molecules are present both in the center and in the colder envelope. Gas-phase HNCO/CH$_3$OH ratios are tentatively correlated with the ratios of NH$_3$ ice over CH$_3$OH ice in the same lines of sight, which is consistent with new gas-grain model predictions. Hot cores are not required to form complex organic molecules, and source temperature and initial ice composition both seem to affect complex organic distributions toward MYSOs. To quantify the relative impact of temperature and initial conditions requires, however, a larger spatially resolved survey of MYSOs with ice detections.

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

Molecular Oxygen in Oort Cloud Comet 1P/Halley

Recently the ROSINA mass spectrometer suite on board the European Space Agency's Rosetta spacecraft discovered an abundant amount of molecular oxygen, O2, in the coma of Jupiter family comet 67P/Churyumov-Gerasimenko of O2/H2O = 3.80+/-0.85%. It could be shown that O2 is indeed a parent species and that the derived abundances point to a primordial origin. One crucial question is whether the O2 abundance is peculiar to comet 67P/Churyumov-Gerasimenko or Jupiter family comets in general or whether also Oort cloud comets such as comet 1P/Halley contain similar amounts of molecular oxygen. We investigated mass spectra obtained by the Neutral Mass Spectrometer instrument obtained during the flyby by the European Space Agency's Giotto probe at comet 1P/Halley. Our investigation indicates that a production rate of O2 of 3.7+/-1.7% with respect to water is indeed compatible with the obtained Halley data and therefore that O2 might be a rather common and abundant parent species.

preprint2015arXiv

Molecular photodissociation

Photodissociation is the dominant removal process of molecules in any region exposed to intense ultraviolet (UV) radiation. This includes diffuse and translucent interstellar clouds, dense photon-dominated regions, high velocity shocks, the surface layers of protoplanetary disks, and cometary and exoplanetary atmospheres. The rate of photodissociation depends on the cross sections for absorption into a range of excited electronic states, as well as on the intensity and shape of the radiation field at each position into the region of interest. Thus, an acccurate determination of the photodissociation rate of even a simple molecule like water or carbon monoxide involves many detailed considerations ranging from its electronic structure to its dissociation dynamics and the specifics of the radiation field that the molecule is exposed to. In this review chapter, each of these steps in determining photodissociation rates is discussed systematically and examples are provided.

preprint2015arXiv

Photodesorption of H2O, HDO, and D2O ice and its impact on fractionation

The HDO/H2O ratio in interstellar gas is often used to draw conclusions on the origin of water in star-forming regions and on Earth. In cold cores and in the outer regions of protoplanetary disks, gas-phase water comes from photodesorption of water ice. We present fitting formulae for implementation in astrochemical models using photodesorption efficiencies for all water ice isotopologues obtained using classical molecular dynamics simulations. We investigate if the gas-phase HDO/H2O ratio reflects that present in the ice or whether fractionation can occur during photodesorption. Probabilities for the top four monolayers are presented for photodesorption of X (X=H,D) atoms, OX radicals, and X2O and HDO molecules following photodissociation of H2O, D2O, and HDO in H2O amorphous ice at temperatures from 10-100 K. Isotope effects are found for all products: (1) H atom photodesorption probabilities from H2O ice are larger than those for D atom photodesorption from D2O ice by a factor of 1.1; the ratio of H and D photodesorbed upon HDO photodissociation is a factor of 2. This process will enrich the ice in deuterium atoms over time; (2) the OD/OH photodesorption ratio upon D2O and H2O photodissociation is on average a factor of 2, but the ratio upon HDO photodissociation is almost constant at unity for all temperatures; (3) D atoms are more effective in kicking out neighbouring water molecules than H atoms. However, the ratio of the photodesorbed HDO and H2O molecules is equal to the HDO/H2O ratio in the ice, therefore, there is no isotope fractionation upon HDO and H2O photodesorption. Nevertheless, the enrichment of the ice in D atoms due to photodesorption can over time lead to an enhanced HDO/H2O ratio in the ice, and, when photodesorbed, also in the gas. The extent to which the ortho/para ratio of H2O can be modified by the photodesorption process is also discussed. (Abridged)

preprint2015arXiv

Resolved gas cavities in transitional disks inferred from CO isotopologues with ALMA

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

preprint2015arXiv

The complex chemistry of outflow cavity walls exposed: the case of low-mass protostars

Complex organic molecules are ubiquitous companions of young low-mass protostars. Recent observations suggest that their emission stems, not only from the traditional hot corino, but also from offset positions. In this work, 2D physicochemical modelling of an envelope-cavity system is carried out. Wavelength-dependent radiative transfer calculations are performed and a comprehensive gas-grain chemical network is used to simulate the physical and chemical structure. The morphology of the system delineates three distinct regions: the cavity wall layer with time-dependent and species-variant enhancements; a torus rich in complex organic ices, but not reflected in gas-phase abundances; and the remaining outer envelope abundant in simpler solid and gaseous molecules. Strongly irradiated regions, such as the cavity wall layer, are subject to frequent photodissociation in the solid phase. Subsequent recombination of the photoproducts leads to frequent reactive desorption, causing gas-phase enhancements of several orders of magnitude. This mechanism remains to be quantified with laboratory experiments. Direct photodesorption is found to be relatively inefficient. If radicals are not produced directly in the icy mantle, the formation of complex organics is impeded. For efficiency, a sufficient number of FUV photons needs to penetrate the envelope; and elevated cool dust temperatures need to enable grain-surface radical mobility. As a result, a high stellar luminosity and a sufficiently wide cavity favor chemical complexity. Furthermore within this paradigm, complex organics are demonstrated to have unique lifetimes and be grouped into early (formaldehyde, ketene, methanol, formic acid, methyl formate, acetic acid, glycolaldehyde) and late (acetaldehyde, dimethyl ether, ethanol) species.

preprint2015arXiv

The molecular composition of the planet-forming regions of protoplanetary disks across the luminosity regime

(Abridged) Near- to mid-IR observations of protoplanetary disks show that the inner regions (<10AU) are rich in small organic volatiles (e.g., C2H2 and HCN). Trends in the data suggest that disks around cooler stars (~3000K) are potentially more carbon- and molecule-rich than their hotter counterparts. Our aims are to explore the composition of the planet-forming region of disks around stars from M dwarf to Herbig Ae and compare with the observed trends. Models of the disk physical structure are coupled with a gas-grain chemical network to map the abundances in the planet-forming zone. N2 self shielding, X-ray-induced chemistry, and initial abundances, are investigated. The composition in the 'observable' atmosphere is compared with that in the midplane where the planet-building reservoir resides. M dwarf disk atmospheres are relatively more molecule rich than those for T Tauri or Herbig Ae disks. The weak far-UV flux helps retain this complexity which is enhanced by X-ray-induced ion-molecule chemistry. N2 self shielding has only a small effect and does not explain the higher C2H2/HCN ratios observed towards cooler stars. The models underproduce the OH/H2O column density ratios constrained in Herbig Ae disks, despite reproducing the absolute value for OH: H2O self shielding only increases this discrepency. The disk midplane content is sensitive to the initial main elemental reservoirs. The gas in the inner disk is generally more carbon rich than the midplane ices and is most significant for disks around cooler stars. The atmospheric C/O ratio appears larger than it actually is when calculated using observable tracers only because gas-phase O2 is predicted to be a significant oxygen reservoir. The models suggest that the gas in the inner regions of disks around cooler stars is more carbon rich; however, calculations of the molecular emission are necessary to confirm the observed trends.

preprint2014arXiv

A Substellar-Mass Protostar and its Outflow of IRAS 15398-3359 Revealed by Subarcsecond-Resolution Observations of H$_2$CO and CCH

Sub-arcsecond ($0.^{\prime\prime}5$) images of H$_2$CO and CCH line emission have been obtained in the $0.8$ mm band toward the low-mass protostar IRAS 15398-3359 in the Lupus 1 cloud as one of the Cycle 0 projects of the Atacama Large Millimeter/Submillimeter Array. We have detected a compact component concentrated in the vicinity of the protostar and a well-collimated outflow cavity extending along the northeast-southwest axis. The inclination angle of the outflow is found to be about $20^\circ$, or almost edge-on, based on the kinematic structure of the outflow cavity. This is in contrast to previous suggestions of a more pole-on geometry. The centrally concentrated component is interpreted by use of a model of the infalling rotating envelope with the estimated inclination angle, and the mass of the protostar is estimated to be less than $0.09\ M_\odot$. Higher spatial resolution data are needed to infer the presence of a rotationally supported disk for this source, hinted at by a weak high-velocity H$_2$CO emission associated with the protostar.

preprint2014arXiv

Effects of Reagent Rotation and Vibration on H + OH (v,j) $\to$ O + H2

The dynamics of the reaction H + OH $\to$ O (3P) + H2 have been studied in a series of quasi-classical trajectory (QCT) calculations and transition state theory (TST) methods using high quality 3A' and 3A'' potential energy surfaces (PESs). Accurate OH (v, j) state resolved cross sections and rate constants on both potential energy surfaces are presented and fitted for OH at (v = 0, j = 0-16) and (v = 1, j = 0-6). The cross sections were calculated for different collisional energies (Ec), ranging from the threshold energy at each specific rovibrational state up to 1.0 eV with step sizes of 0.1 eV or less. They increase steeply with collision energy when the barrier to reaction can be overcome, after which the cross sections stay nearly constant with energy. State resolved rate constants in the temperature range 200-2500 K are presented based on the cross sections. Total thermal rate constants were calculated by summing the rates for reaction on the 3A' and 3A'' potential energy surfaces weighted by 1/3 and taking into account the thermal populations of the rovibrational states of the OH molecules. It is shown that the improved canonical variational transition (CVT) treatments with the approximation of zero-curvature tunneling (ZCT) or small-curvature tunneling (SCT) produce results more in accord with the QCT results than the TST and CVT methods. The reactions are governed by the direct reaction mechanism. The rate constants for OH in excited vibrational and rotational states are orders of magnitude larger than the thermal rate constants, which needs to be taken into account in astrochemical models.

preprint2014arXiv

Herschel/PACS view of disks around low-mass stars and brown dwarfs in the TW Hya association

We conducted Herschel/PACS observations of five very low-mass stars or brown dwarfs located in the TW Hya association with the goal of characterizing the properties of disks in the low stellar mass regime. We detected all five targets at $70\,μ{\rm{m}}$ and $100\,μ{\rm{m}}$ and three targets at $160\,μ{\rm{m}}$. Our observations, combined with previous photometry from 2MASS, WISE, and SCUBA-2, enabled us to construct SEDs with extended wavelength coverage. Using sophisticated radiative transfer models, we analyzed the observed SEDs of the five detected objects with a hybrid fitting strategy that combines the model grids and the simulated annealing algorithm and evaluated the constraints on the disk properties via the Bayesian inference method. The modelling suggests that disks around low-mass stars and brown dwarfs are generally flatter than their higher mass counterparts, but the range of disk mass extends to well below the value found in T Tauri stars, and the disk scale heights are comparable in both groups. The inferred disk properties (i.e., disk mass, flaring, and scale height) in the low stellar mass regime are consistent with previous findings from large samples of brown dwarfs and very low-mass stars. We discuss the dependence of disk properties on their host stellar parameters and find a significant correlation between the Herschel far-IR fluxes and the stellar effective temperatures, probably indicating that the scaling between the stellar and disk masses (i.e., $M_{\rm{disk}} \propto M_{\star}$) observed mainly in low-mass stars may extend down to the brown dwarf regime.

preprint2014arXiv

Isotope selective photodissociation of N2 by the interstellar radiation field and cosmic rays

Photodissociation of 14N2 and 14N15N occurs in interstellar clouds, circumstellar envelopes, protoplanetary discs, and other environments due to UV radiation from stellar sources and the presence of cosmic rays. This source of N atoms initiates the formation of complex N-bearing species and influences their isotopic composition. To study the photodissociation rates of 14N15N by UV continuum radiation and both isotopologues in a field of cosmic ray induced photons. To determine the effect of these on the isotopic composition of more complex molecules. High-resolution photodissociation cross sections of N2 are used from an accurate and comprehensive quantum- mechanical model of the molecule based on laboratory experiments. A similarly high-resolution spectrum of H2 emission following interactions with cosmic rays has been constructed. The spectroscopic data are used to calculate dissociation rates which are input into isotopically differentiated chemical models, describing an interstellar cloud and a protoplanetary disc. The dissociation rate of 14N15N in a Draine field assuming 30K excitation is 1.73x10-10s-1 and the rate due to cosmic rays assuming an H2 ionisation rate of 10-16s-1 is about 10-15s-1, with up to a factor of 10 difference between isotopologues. Shielding functions for 14N15N by 14N2, H2, and H are presented. Incorporating these into an interstellar cloud model, an enhancement of the atomic 15N/14N ratio is obtained due to the self-shielding of external radiation at an extinction of about 1.5 mag. This effect is larger where grain growth has reduced the opacity of dust to ultraviolet radiation. The transfer of photolytic isotopic fractionation N2 to other molecules is significant in a disc model, and is species dependent with 15N enhancement approaching a factor of 10 for HCN.

preprint2014arXiv

Methanol Along the Path from Envelope to Protoplanetary Disc

Interstellar methanol is considered to be a parent species of larger, more complex organic molecules. A physicochemical simulation of infalling parcels of matter is performed for a low-mass star-forming system to trace the chemical evolution from cloud to disc. An axisymmetric 2D semi-analytic model generates the time-dependent density and velocity distributions, and full continuum radiative transfer is performed to calculate the dust temperature and the UV radiation field at each position as a function of time. A comprehensive gas-grain chemical network is employed to compute the chemical abundances along infall trajectories. Two physical scenarios are studied, one in which the dominant disc growth mechanism is viscous spreading, and another in which continuous infall of matter prevails. The results show that the infall path influences the abundance of methanol entering each type of disc, ranging from complete loss of methanol to an enhancement by a factor of > 1 relative to the prestellar phase. Critical chemical processes and parameters for the methanol chemistry under different physical conditions are identified. The exact abundance and distribution of methanol is important for the budget of complex organic molecules in discs, which will be incorporated into forming planetary system objects such as protoplanets and comets. These simulations show that the comet-forming zone contains less methanol than in the precollapse phase, which is dominantly of prestellar origin, but also with additional layers built up in the envelope during infall. Such intriguing links will soon be tested by upcoming data from the Rosetta mission.

preprint2014arXiv

Photodissociation and chemistry of N$_2$ in the circumstellar envelope of carbon-rich AGB stars

The envelopes of AGB stars are irradiated externally by ultraviolet photons; hence, the chemistry is sensitive to the photodissociation of N$_2$ and CO, which are major reservoirs of nitrogen and carbon, respectively. The photodissociation of N$_2$ has recently been quantified by laboratory and theoretical studies. Improvements have also been made for CO photodissociation. For the first time, we use accurate N$_2$ and CO photodissociation rates and shielding functions in a model of the circumstellar envelope of the carbon-rich AGB star, IRC +10216. We use a state-of-the-art chemical model of an AGB envelope, the latest CO and N$_2$ photodissociation data, and a new method for implementing molecular shielding functions in full spherical geometry with isotropic incident radiation. We compare computed column densities and radial distributions of molecules with observations. The transition of N$_2$ $\to$ N (also, CO $\to$ C $\to$ C$^+$) is shifted towards the outer envelope relative to previous models. This leads to different column densities and radial distributions of N-bearing species, especially those species whose formation/destruction processes largely depend on the availability of atomic or molecular nitrogen, for example, C$_n$N ($n$=1, 3, 5), C$_n$N$^-$ ($n$=1, 3, 5), HC$_n$N ($n$=1, 3, 5, 7, 9), H$_2$CN and CH$_2$CN. The chemistry of many species is directly or indirectly affected by the photodissociation of N$_2$ and CO, especially in the outer shell of AGB stars where photodissociation is important. Thus, it is important to include N$_2$ and CO shielding in astrochemical models of AGB envelopes and other irradiated environments. In general, while differences remain between our model of IRC +10216 and the observed molecular column densities, better agreement is found between the calculated and observed radii of peak abundance.

preprint2014arXiv

Protoplanetary disk masses from CO isotopologues line emission

One of the methods for deriving disk masses relies on direct observations of the gas, whose bulk mass is in the outer cold ($T\lesssim30$K) regions. This zone can be well traced by rotational lines of less abundant CO isotopologues, that probe the gas down to the midplane. The total CO gas mass is then obtained with the isotopologue ratios taken to be constant at the elemental isotope values found in the local ISM. This approach is however imprecise, because isotope selective processes are ignored. The aim of this work is an isotopologue selective treatment of CO isotopologues, in order to obtain a more accurate determination of disk masses. The isotope-selective photodissociation, the main process controlling the abundances of CO isotopologues in the CO-emissive layer, is properly treated for the first time in a full disk model (DALI, Bruderer et al. 2012; Bruderer 2013). The chemistry, thermal balance, line and continuum radiative transfer are all considered together with a chemical network that treats $^{13}$CO, C$^{18}$O, C$^{17}$O, isotopes of all included atoms, and molecules, as independent species. Isotope selective processes lead to regions in the disk where the isotopologues abundance ratios are considerably different from the elemental ratio. The results of this work show that considering CO isotopologue ratios as constants can lead to an underestimate of disk masses by up to almost two orders of magnitude if grains have grown to larger sizes. This may explain observed discrepancies in mass determinations from different tracers. The dependence of the various isotopologues emission on stellar and disk parameters is investigated. Including CO isotope selective processes is crucial to determine the gas mass of the disk accurately (through ALMA observations) and thus to provide the amount of gas which may eventually form planets or change the dynamics of forming planetary systems.

preprint2014arXiv

Ro-vibrational excitation of an organic molecule (HCN) in protoplanetary disks

(Abridged) Organic molecules are important constituents of protoplanetary disks. Their ro-vibrational lines observed in the near- and mid-infrared are commonly detected toward T Tauri disks. These lines are the only way to probe the chemistry in the inner few au where terrestrial planets form. To understand this chemistry, accurate molecular abundances have to be determined. This is complicated by excitation effects. Most analyses so far have made the assumption of local thermal equilibrium (LTE). Starting from estimates for the collisional rate coefficients of HCN, non-LTE slab models of the HCN emission were calculated to study the importance of different excitation mechanisms. Using a new radiative transfer model, the HCN emission from a full two-dimensional disk was then modeled to study the effect of the non-LTE excitation, together with the line formation. We ran models tailored to the T Tauri disk AS 205 (N) where HCN lines in both the 3 μm and 14 μm bands have been observed by VLT-CRIRES and the Spitzer Space Telescope. Reproducing the observed 3 μm / 14 μm flux ratios requires very high densities and kinetic temperatures ($n > 10^{14}$ cm$^{-3}$ and $T > 750$ K), if only collisional excitation is accounted for. Radiative pumping can, however, excite the lines easily out to considerable radii $\sim$ 10 au. Consequently, abundances derived from LTE and non-LTE models do not differ by more than a factor of about 3. Models with both a strongly enhanced abundance within $\sim$ 1 au (jump abundance) and constant abundance can reproduce the current observations, but future observations with the MIRI instrument on JWST and METIS on the E-ELT can easily distinguish between the scenarios and test chemical models. Depending on the scenario, ALMA can detect rotational lines within vibrationally excited levels.

preprint2014arXiv

Warm formaldehyde in the Oph IRS 48 transitional disk

Simple molecules like H2CO and CH3OH in protoplanetary disks are the starting point for the production of more complex organic molecules. So far, the observed chemical complexity in disks has been limited due to freeze out of molecules onto grains in the bulk of the cold outer disk. Complex molecules can be studied more directly in transitional disks with large inner holes, as these have a higher potential of detection, through UV heating of the outer disk and the directly exposed midplane at the wall. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 9 (~680 GHz) line data of the transitional disk Oph IRS 48, previously shown to have a large dust trap, to search for complex molecules in regions where planetesimals are forming. We report the detection of the H2CO 9(1,8)-8(1,7) line at 674 GHz, which is spatially resolved as a semi-ring at ~60 AU radius centered south from the star. The inferred H2CO abundance is ~10^{-8} derived by combining a physical disk model of the source with a non-LTE excitation calculation. Upper limits for CH3OH lines in the same disk give an abundance ratio H2CO/CH3OH>0.3, which points to both ice formation and gas-phase routes playing a role in the H2CO production. Upper limits on the abundances of H13CO+, CN and several other molecules in the disk are also derived and found to be consistent with full chemical models. The detection of the H2CO line demonstrates the start of complex organic molecules in a planet-forming disk. Future ALMA observations should be able to push down the abundance detection limits of other molecules by 1-2 orders of magnitude and test chemical models of organic molecules in (transitional) disks.

preprint2014arXiv

Water and methanol in low-mass protostellar outflows: gas-phase synthesis, ice sputtering and destruction

Water in outflows from protostars originates either as a result of gas-phase synthesis from atomic oxygen at T > 200 K, or from sputtered ice mantles containing water ice. We aim to quantify the contribution of the two mechanisms that lead to water in outflows, by comparing observations of gas-phase water to methanol (a grain surface product) towards three low-mass protostars in NGC1333. In doing so, we also quantify the amount of methanol destroyed in outflows. To do this, we make use of JCMT and Herschel-HIFI data of H2O, CH3OH and CO emission lines and compare them to RADEX non-LTE excitation simulations. We find up to one order of magnitude decrease in the column density ratio of CH3OH over H2O as the velocity increases in the line wings up to ~15 km/s. An independent decrease in X(CH3OH) with respect to CO of up to one order of magnitude is also found in these objects. We conclude that gas-phase formation of H2O must be active at high velocities (above 10 km/s, relative to the source velocity) to re-form the water destroyed during sputtering. In addition, the transition from sputtered water at low velocities to formed water at high velocities must be gradual. We place an upper limit of two orders of magnitude on the destruction of methanol by sputtering effects.

preprint2014arXiv

Water: from clouds to planets

Results from recent space missions, in particular Spitzer and Herschel, have lead to significant progress in our understanding of the formation and transport of water from clouds to disks, planetesimals, and planets. In this review, we provide the underpinnings for the basic molecular physics and chemistry of water and outline these advances in the context of water formation in space, its transport to a forming disk, its evolution in the disk, and finally the delivery to forming terrestrial worlds and accretion by gas giants. Throughout, we pay close attention to the disposition of water as vapor or solid and whether it might be subject to processing at any stage. The context of the water in the solar system and the isotopic ratios (D/H) in various bodies are discussed as grounding data point for this evolution. Additional advances include growing knowledge of the composition of atmospheres of extra-solar gas giants, which may be influenced by the variable phases of water in the protoplanetary disk. Further, the architecture of extra-solar systems leaves strong hints of dynamical interactions, which are important for the delivery of water and subsequent evolution of planetary systems. We conclude with an exploration of water on Earth and note that all of the processes and key parameters identified here should also hold for exoplanetary systems.

preprint2013arXiv

A Keplerian disk around a Class 0 source: ALMA observations of VLA1623A

Context: Rotationally supported disks are critical in the star formation process. The questions of when do they form and what factors influence or hinder their formation have been studied but are largely unanswered. Observations of early stage YSOs are needed to probe disk formation. Aims: VLA1623 is a triple non-coeval protostellar system, with a weak magnetic field perpendicular to the outflow, whose Class 0 component, VLA1623A, shows a disk-like structure in continuum with signatures of rotation in line emission. We aim to determine whether this structure is in part or in whole a rotationally supported disk, i.e. a Keplerian disk, and what are its characteristics. Methods: ALMA Cycle 0 Early Science 1.3 mm continuum and C$^{18}$O (2-1) observations in the extended configuration are presented here and used to perform an analysis of the disk-like structure using PV diagrams and thin disk modelling with the addition of foreground absorption. Results: The PV diagrams of the C$^{18}$O line emission suggest the presence of a rotationally supported component with a radius of at least 50 AU. Kinematical modelling of the line emission shows that the disk out to 180 AU is actually rotationally supported, with the rotation being well described by Keplerian rotation out to at least 150 AU, and the central source mass to be $\sim$0.2 M$_{sun}$ for an inclination of 55$^{\circ}$. Pure infall and conserved angular momentum rotation models are excluded. Conclusions: VLA1623A, a very young Class 0 source, presents a disk with an outer radius $R_{\rm out}$ = 180 AU with a Keplerian velocity structure out to at least 150 AU. The weak magnetic fields and recent fragmentation in this region of rho Ophiuchus may have played a lead role in the formation of the disk.

preprint2013arXiv

A Major Asymmetric Dust Trap in a Transition Disk

The statistics of discovered exoplanets suggest that planets form efficiently. However, there are fundamental unsolved problems, such as excessive inward drift of particles in protoplanetary disks during planet formation. Recent theories invoke dust traps to overcome this problem. We report the detection of a dust trap in the disk around the star Oph IRS 48 using observations from the Atacama Large Millimeter/submillimeter Array (ALMA). The 0.44-millimeter-wavelength continuum map shows high-contrast crescent-shaped emission on one side of the star originating from millimeter-sized grains, whereas both the mid-infrared image (micrometer-sized dust) and the gas traced by the carbon monoxide 6-5 rotational line suggest rings centered on the star. The difference in distribution of big grains versus small grains/gas can be modeled with a vortex-shaped dust trap triggered by a companion.

preprint2013arXiv

A recent accretion burst in the low-mass protostar IRAS 15398-3359: ALMA imaging of its related chemistry

Low-mass protostars have been suggested to show highly variable accretion rates through-out their evolution. Such changes in accretion, and related heating of their ambient envelopes, may trigger significant chemical variations on different spatial scales and from source-to-source. We present images of emission from C17O, H13CO+, CH3OH, C34S and C2H toward the low-mass protostar IRAS 15398-3359 on 0.5" (75 AU diameter) scales with the Atacama Large Millimeter/submillimeter Array (ALMA) at 340 GHz. The resolved images show that the emission from H13CO+ is only present in a ring-like structure with a radius of about 1-1.5" (150-200 AU) whereas the CO and other high dipole moment molecules are centrally condensed toward the location of the central protostar. We propose that HCO+ is destroyed by water vapor present on small scales. The origin of this water vapor is likely an accretion burst during the last 100-1000 years increasing the luminosity of IRAS 15398-3359 by a factor of 100 above its current luminosity. Such a burst in luminosity can also explain the centrally condensed CH3OH and extended warm carbon-chain chemistry observed in this source and furthermore be reflected in the relative faintness of its compact continuum emission compared to other protostars.

preprint2013arXiv

Deep observations of O2 toward a low-mass protostar with Herschel-HIFI

According to traditional gas-phase chemical models, O2 should be abundant in molecular clouds, but until recently, attempts to detect interstellar O2 line emission with ground- and space-based observatories have failed. Following the multi-line detections of O2 with low abundances in the Orion and rho Oph A molecular clouds with Herschel, it is important to investigate other environments, and we here quantify the O2 abundance near a solar-mass protostar. Observations of O2, at 487 GHz toward a deeply embedded low-mass Class 0 protostar, NGC 1333-IRAS 4A, are presented, using the HIFI instrument on the Herschel Space Observatory. Complementary data of the chemically related NO and CO molecules are obtained as well. The high spectral resolution data are analysed using radiative transfer models to infer column densities and abundances, and are tested directly against full gas-grain chemical models. The deep HIFI spectrum fails to show O2 at the velocity of the dense protostellar envelope, implying one of the lowest abundance upper limits of O2/H2 at <6x10^-9 (3 sigma). However, a tentative (4.5 sigma) detection of O2 is seen at the velocity of the surrounding NGC 1333 molecular cloud, shifted by 1 km/s relative to the protostar. For the protostellar envelope, pure gas-phase models and gas-grain chemical models require a long pre-collapse phase (~0.7-1x10^6 years), during which atomic and molecular oxygen are frozen out onto dust grains and fully converted to H2O, to avoid overproduction of O2 in the dense envelope. The same model also reproduces the limits on the chemically related NO molecule. The tentative detection of O2 in the surrounding cloud is consistent with a low-density PDR model with small changes in reaction rates. The low O2 abundance in the collapsing envelope around a low-mass protostar suggests that the gas and ice entering protoplanetary disks is very poor in O2.

preprint2013arXiv

DIGIT survey of far-infrared lines from protoplanetary discs: II. CO

CO is an important component of a protoplanetary disc as it is one of the most abundant gas phase species. Furthermore, observations of the CO transitions can be used as a diagnostic of the gas, tracing conditions in both the inner and outer disc. We present Herschel/PACS spectroscopy of a sample of 22 Herbig Ae/Be (HAEBEs) and 8 T Tauri stars (TTS), covering the pure rotational CO transitions from J = 14-13 up to J = 49-48. CO is detected in only 5 HAEBEs, AB Aur, HD 36112, HD 97048, HD 100546 and IRS 48, and in 4 TTS, AS 205, S CrA, RU Lup, and DG Tau. The highest transition seen is J = 36-35, with Eup of 3669 K, detected in HD 100546 and DG Tau. We construct rotational diagrams for the discs with at least 3 CO detections to derive Trot, and find average temperatures of 270 K for the HAEBEs and 485 K for the TTS. HD 100546 requires an extra temperature component with Trot ~ 900-1000 K, suggesting a range of temperatures in its disc atmosphere, consistent with thermo-chemical disc models. In HAEBEs, the objects with CO detections all have flared discs in which the gas and dust are thermally decoupled. We use a small model grid to analyse our observations and find that an increased amount of flaring means higher line flux, as it increases the mass in warm gas. CO is not detected in our flat discs as the emission is below the detection limit. We find that HAEBE sources with CO detections have high Luv and strong PAH emission, again connected to the heating of the gas. In TTS, the objects with CO detections are all sources with evidence for a disc wind or outflow. For both groups of objects, sources with CO detections generally have high UV luminosity (either stellar in HAEBEs or due to accretion in TTS), but this is not a sufficient condition for the detection of the far-IR CO lines.

preprint2013arXiv

Embedded Protostars in the Dust, Ice, and Gas In Time (DIGIT) Key Program: Continuum SEDs, and an Inventory of Characteristic Far-Infrared Lines from PACS Spectroscopy

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

preprint2013arXiv

Evolution of CO lines in time-dependent models of protostellar disk formation

(Abridged) Star and planet formation theories predict an evolution in the density, temperature, and velocity structure as the envelope collapses and forms an accretion disk. The aim of this work is to model the evolution of the molecular excitation, line profiles, and related observables during low-mass star formation. Specifically, the signatures of disks during the deeply embedded stage are investigated. Semi-analytic 2D axisymmetric models have been used to describe the evolution of the density, stellar mass, and luminosity from the pre-stellar to the T-Tauri phase. A full radiative transfer calculation is carried out to accurately determine the time-dependent dust temperatures and CO abundance structure. We present non-LTE near-IR, FIR, and submm lines of CO have been simulated at a number of time steps. In contrast to the dust temperature, the CO excitation temperature derived from submm/FIR lines does not vary during the protostellar evolution, consistent with C18O observations obtained with Herschel and from ground-based telescopes. The near-IR spectra provide complementary information to the submm lines by probing not only the cold outer envelope but also the warm inner region. The near-IR high-J (>8) absorption lines are particularly sensitive to the physical structure of the inner few AU, which does show evolution. High signal-to-noise ratio subarcsec resolution data with ALMA are needed to detect the presence of small rotationally supported disks during the Stage 0 phase and various diagnostics are discussed.

preprint2013arXiv

Gas structure inside dust cavities of transition disks: Oph IRS 48 observed by ALMA

(Abridged) Transition disks are recognized by the absence of emission of small dust grains inside a radius of up to several 10s of AUs. Due to the lack of angular resolution and sensitivity, the gas content of such dust holes has not yet been determined, but is of importance to constrain the mechanism leading to the dust holes. Transition disks are thought to currently undergo the process of dispersal, setting an end to the giant planet formation process. We present new high-resolution observations with the Atacama Large Millimeter/ submillimeter Array (ALMA) of gas lines towards the transition disk Oph IRS 48 previously shown to host a large dust trap. ALMA has detected the $J=6-5$ line of $^{12}$CO and C$^{17}$O around 690 GHz (434 $μ$m) at a resolution of $\sim$0.25$''$ corresponding to $\sim$30 AU (FWHM). The observed gas lines are used to set constraints on the gas surface density profile. New models of the physical-chemical structure of gas and dust in Oph IRS 48 are developed to reproduce the CO line emission together with the spectral energy distribution (SED) and the VLT-VISIR 18.7 $μ$m dust continuum images. Integrated intensity cuts and the total spectrum from models having different trial gas surface density profiles are compared to observations. Using the derived surface density profiles, predictions for other CO isotopologues are made, which can be tested by future ALMA observations of the object. The derived gas surface density profile points to the clearing of the cavity by one or more massive planet/companion rather than just photoevaporation or grain-growth.

preprint2013arXiv

High-J CO survey of low-mass protostars observed with Herschel-HIFI

The evolution of deeply embedded young stellar objects (YSOs) has traditionally been traced through dust continuum spectral energy distributions (SEDs), but the use of CO excitation as an evolutionary probe has not yet been explored due to lack of high-J CO observations. The aim is to constrain the physical characteristics (excitation, kinematics, column density) of the warm gas toward 26 low-mass Class 0 and I YSOs using spectrally-resolved Herschel Space Observatory data of high-J lines of CO. Data are complemented by ground-based observations from APEX and the JCMT to compare those with the colder gas traced by lower-J CO lines. This is the first large spectrally resolved high-J CO survey conducted for these types of sources. The median excitation temperatures for 12CO, 13CO and C18O derived from single-temperature fits to the J_u=2-10 integrated intensities are ~70K, 48K and 37K, respectively, with no significant difference between Class 0 and I sources and no trend with M_env or L_bol. Thus, in contrast with the continuum SEDs, the spectral line energy distributions (SLEDs) do not show an evolution during the embedded stage. In contrast, the integrated line intensities of all CO isotopologs show a clear decrease with evolutionary stage as the envelope is dispersed. The H2O 110-101/CO 10-9 intensity ratio does not change significantly with velocity, in contrast with the H2O/CO 3-2 ratio, indicating that CO 10-9 is the lowest transition for which the line wings probe the same warm shocked gas as H2O. Modeling of the full suite of C18O lines indicates an abundance profile for Class 0 sources that is consistent with a freeze-out zone below 25 K and evaporation at higher temperatures, with but some fraction of the CO transformed into other species in the cold phase. In contrast, observations for two Class I sources in Ophiuchus are consistent with a constant high CO abundance profile.

preprint2013arXiv

Hot water in the inner 100 AU of the Class 0 protostar NGC1333 IRAS2A

Evaporation of water ice above 100 K in the inner few 100 AU of low-mass embedded protostars (the so-called hot core) should produce quiescent water vapor abundances of ~10^-4 relative to H2. Observational evidence so far points at abundances of only a few 10^-6. However, these values are based on spherical models, which are known from interferometric studies to be inaccurate on the relevant spatial scales. Are hot cores really that much drier than expected, or are the low abundances an artifact of the inaccurate physical models? We present deep velocity-resolved Herschel-HIFI spectra of the 3(12)-3(03) lines of H2-16O and H2-18O (1097 GHz, Eup/k = 249 K) in the low-mass Class 0 protostar NGC1333 IRAS2A. A spherical radiative transfer model with a power-law density profile is unable to reproduce both the HIFI data and existing interferometric data on the H2-18O 3(13)-2(20) line (203 GHz, Eup/k = 204 K). Instead, the HIFI spectra likely show optically thick emission from a hot core with a radius of about 100 AU. The mass of the hot core is estimated from the C18O J=9-8 and 10-9 lines. We derive a lower limit to the hot water abundance of 2x10^-5, consistent with the theoretical predictions of ~10^-4. The revised HDO/H2O abundance ratio is 1x10^-3, an order of magnitude lower than previously estimated.

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

Interstellar water chemistry: from laboratory to observations

Water is observed throughout the universe, from diffuse interstellar clouds to protoplanetary disks around young stars, and from comets in our own solar system and exoplanetary atmospheres to galaxies at high redshifts. This review summarizes the spectroscopy and excitation of water in interstellar space as well as the basic chemical processes that form and destroy water under interstellar conditions. Three major routes to water formation are identified: low temperature ion-molecule chemistry, high-temperature neutral-neutral chemistry and gas-ice chemistry. The rate coefficients of several important processes entering the networks are discussed in detail; several of them have been determined only in the last decade through laboratory experiments and theoretical calculations. Astronomical examples of each of the different chemical routes are presented using data from powerful new telescopes, in particular the Herschel Space Observatory. Basic chemical physics studies remain critically important to analyze astronomical data.

preprint2013arXiv

Measuring Protoplanetary Disk Accretion with H I Pf-beta

In this work, we introduce the use of H I Pfund-beta (4.6538 micron) as a tracer of mass accretion from protoplanetary disks onto young stars. Pf-beta was serendipitously observed in NIRSPEC and CRIRES surveys of CO fundamental emission, amounting to a sample size of 120 young stars with detected Pf-beta emission. Using a subsample of disks with previously measured accretion luminosities, we show that Pf-beta line luminosity is well correlated with accretion luminosity over a range of at least three orders of magnitude. We use this correlation to derive accretion luminosities for all 120 targets, 65 of which are previously unreported in the literature. We also report accretion rates for 67 targets, 16 previously unmeasured. Our large sample size and our ability to probe high extinction values allow for relatively unbiased comparisons between different types of disks. We find that the transitional disks in our sample have lower than average Pf-beta line luminosities, and thus accretion luminosities, at a marginally significant level. We also show that high Pf-beta equivalent width is a signature of transitional disks with high inner disk gas/dust ratios. In contrast, we find that disks with signatures of slow disk winds have Pf-beta luminosities comparable to those of other disks in our sample. Finally, we investigate accretion rates for stage I disks, including significantly embedded targets. We find that stage I and stage II disks have statistically indistinguishable Pf-beta line luminosities, implying similar accretion rates, and that the accretion rates of stage I disks are too low to be consistent with quiescent accretion. Our results are instead consistent with both observational and theoretical evidence that stage I objects experience episodic, rather than quiescent, accretion.

preprint2013arXiv

Molecular Dynamics Simulations of CO2 Formation in Interstellar Ices

CO2 ice is one of the most abundant components in ice-coated interstellar ices besides H2O and CO, but the most favorable path to CO2 ice is still unclear. Molecular dynamics calculations on the ultraviolet photodissociation of different kinds of CO-H2O ice systems have been performed at 10 K in order to demonstrate that the reaction between CO and an OH molecule resulting from H2O photodissociation through the first excited state is a possible route to form CO2 ice. However, our calculations, which take into account different ice surface models, suggest that there is another product with a higher formation probability ((3.00+-0.07)x10-2), which is the HOCO complex, whereas the formation of CO2 has a probability of only (3.6+-0.7)x10-4. The initial location of the CO is key to obtain reaction and form CO2: the CO needs to be located deep into the ice. The HOCO complex becomes trapped in the cold ice surface in the trans-HOCO minimum because it quickly loses its internal energy to the surrounding ice, preventing further reaction to H + CO2. Several laboratory experiments have been carried out recently and they confirm that CO2 can also be formed through other, different routes. Here we compare our theoretical results with the data available from experi- ments studying the formation of CO2 through a similar pathway as ours, even though the initial conditions were not exactly the same. Our results also show that the HCO van der Waals complex can be formed through the interaction of CO with the H atom that is formed as a product of H2O photodissociation. Thus, the reaction of the H atom photofragment following H2O photodissociation with CO can be a possible route to form HCO ice.

preprint2013arXiv

Neutral and Ionized Hydrides in Star-forming Regions -- Observations with Herschel/HIFI

The cosmic abundance of hydrides depends critically on high-energy UV, X-ray, and particle irradiation. Here we study hydrides in star-forming regions where irradiation by the young stellar object can be substantial, and density and temperature can be much enhanced over interstellar values. Lines of OH, CH, NH, SH and their ions OH+, CH+, NH+, SH+, H2O+, and H3O+ were observed in star-forming regions by the HIFI spectrometer onboard the Herschel Space Observatory. Molecular column densities are derived from observed ground-state lines, models, or rotational diagrams. We report here on two prototypical high-mass regions, AFGL 2591 and W3 IRS5, and compare them to chemical calculations making assumptions on the high-energy irradiation. A model assuming no ionizing protostellar emission is compared with (i) a model assuming strong protostellar X-ray emission and (ii) a two-dimensional (2D) model including emission in the far UV (FUV, 6 -- 13.6 eV) irradiating the outflow walls that separate the outflowing gas and infalling envelope material. We confirm that the effect of FUV in two dimensional models with enlarged irradiated surfaces is clearly noticeable. A molecule that is very sensitive to FUV irradiation is CH+, enhanced in abundance by more than 5 orders of magnitude. The HIFI observations of CH+ lines agree with the two-dimensional FUV model by Bruderer et al. which computes abundances, non-LTE excitation and line radiative transfer.{Ref 20} It is concluded that CH+ is a good FUV tracer in star-forming regions. The effect of potential X-ray irradiation is not excluded, but cannot be demonstrated by the present data.

preprint2013arXiv

Outflow forces of low mass embedded objects in Ophiuchus: a quantitative comparison of analysis methods

The outflow force of molecular bipolar outflows is a key parameter in theories of young stellar feedback on their surroundings. The focus of many outflow studies is the correlation between the outflow force, bolometric luminosity and envelope mass. However, it is difficult to combine the results of different studies in large evolutionary plots over many orders of magnitude due to the range of data quality, analysis methods and corrections for observational effects such as opacity and inclination. We aim to determine the outflow force for a sample of low luminosity embedded sources. We will quantify the influence of the analysis method and the assumptions entering the calculation of the outflow force. We use the James Clerk Maxwell Telescope to map 12CO J=3-2 over 2'x2' regions around 16 Class I sources of a well-defined sample in Ophiuchus at 15" resolution. The outflow force is then calculated using seven different methods differing e.g. in the use of intensity-weighted emission and correction factors for inclination. The results from the analysis methods differ from each other by up to a factor of 6, whereas observational properties and choices in the analysis procedure affect the outflow force by up to a factor of 4. For the sample of Class I objects, bipolar outflows are detected around 13 sources including 5 new detections, where the three non-detections are confused by nearby outflows from other sources. When combining outflow forces from different studies, a scatter by up to a factor of 5 can be expected. Although the true outflow force remains unknown, the separation method (separate calculation of dynamical time and momentum) is least affected by the uncertain observational parameters. The correlations between outflow force, bolometric luminosity and envelope mass are further confirmed down to low luminosity sources.

preprint2013arXiv

Photodissociation of interstellar N2

Molecular nitrogen is one of the key species in the chemistry of interstellar clouds and protoplanetary disks and the partitioning of nitrogen between N and N2 controls the formation of more complex prebiotic nitrogen-containing species. The aim of this work is to gain a better understanding of the interstellar N2 photodissociation processes based on recent detailed theoretical and experimental work and to provide accurate rates for use in chemical models. We simulated the full high-resolution line-by-line absorption + dissociation spectrum of N2 over the relevant 912-1000 Å wavelength range, by using a quantum-mechanical model which solves the coupled-channels Schrödinger equation. The simulated N2 spectra were compared with the absorption spectra of H2, H, CO, and dust to compute photodissociation rates in various radiation fields and shielding functions. The effects of the new rates in interstellar cloud models were illustrated for diffuse and translucent clouds, a dense photon dominated region and a protoplanetary disk.

preprint2013arXiv

Rotationally-supported disks around Class I sources in Taurus: disk formation constraints

(Abridged) Disks are observed around pre-main sequence stars, but how and when they form is still heavily debated. While disks around young stellar objects have been identified through thermal dust emission, spatially and spectrally resolved molecular line observations are needed to determine their nature. We present subarcsecond observations of dust and gas toward four Class I low-mass young stellar objects in Taurus. The 13CO and C18O J=2-1 transitions at 220 GHz were observed with the Plateau de Bure Interferometer at a spatial resolution of ~0.8'' and analyzed using uv-space position velocity diagrams to determine the nature of their observed velocity radient. Rotationally supported disks (RSDs) are detected around 3 of the 4 Class I sources studied. The derived masses identify them as Stage I objects; i.e., their stellar mass is higher than their envelope and disk masses. The outer radii of the Keplerian disks toward our sample of Class I sources are <~ 100 AU. The lack of on-source C18O emission for TMR1 puts an upper limit of 50 AU on its size. Flattened structures at radii > 100 AU around these sources are dominated by infalling motion (v propto r^-1). A large-scale envelope model is required to estimate the basic parameters of the flattened structure from spatially resolved continuum data. Similarities and differences between the gas and dust disk are discussed. Combined with literature data, the sizes of the RSDs around Class I objects are best described with evolutionary models with an initial rotation of 10^-14 Hz and slow sound speeds. Based on the comparison of gas and dust disk masses, little CO is frozen out within 100 AU in these disks. RSDs with radii up to 100 AU are present around Class I embedded objects. Larger surveys of both Class 0 and I objects are needed to determine whether most disks form late or early in the embedded phase.

preprint2013arXiv

The water abundance behind interstellar shocks: results from $Herschel$/PACS and $Spitzer$/IRS observations of H$_2$O, CO, and H$_2$

We have investigated the water abundance in shock-heated molecular gas, making use of $Herschel$ measurements of far-infrared CO and H$_2$O line emissions in combination with $Spitzer$ measurements of mid-IR H$_2$ rotational emissions. We present far-infrared line spectra obtained with $Herschel$'s PACS instrument in range spectroscopy mode towards two positions in the protostellar outflow NGC 2071 and one position each in the supernova remnants W28 and 3C391. These spectra provide unequivocal detections, at one or more positions, of 12 rotational lines of water, 14 rotational lines of CO, 8 rotational lines of OH (4 lambda doublets), and 7 fine-structure transitions of atoms or atomic ions. We first used a simultaneous fit to the CO line fluxes, along with H$_2$ rotational line fluxes measured previously by $Spitzer$, to constrain the temperature and density distribution within the emitting gas; and we then investigated the water abundances implied by the observed H$_2$O line fluxes. The water line fluxes are in acceptable agreement with standard theoretical models for nondissociative shocks that predict the complete vaporization of grain mantles in shocks of velocity $v \sim 25$ km/s, behind which the characteristic gas temperature is $\sim 1300$ K and the H$_2$O/CO ratio is 1.2

preprint2013arXiv

VLT-CRIRES Survey of Rovibrational CO Emission from Protoplanetary Disks

We present a large, comprehensive survey of rovibrational CO line emission at 4.7 micron from 69 protoplanetary disks, obtained with CRIRES on the ESO Very Large Telescope at the highest available spectral resolving power (R=95,000, v=3.2 km/s). The CO fundamental band (Delta v=1) is a well-known tracer of warm gas in the inner, planet-forming regions of gas-rich disks around young stars, with the lines formed in the super-heated surfaces of the disks at radii of 0.1-10 AU. Our high spectral resolution data provide new insight into the kinematics of the inner disk gas. Pure double-peaked Keplerian profiles are surprisingly uncommon, beyond the frequency expected based on disk inclination. The majority of the profiles are consistent with emission from a disk plus a slow (few km/s) molecular disk wind. This is evidenced by analysis of different categories as well as an overall tendency for line profiles to have excess blue emission. Weak emission lines from isotopologues and vibrationally excited levels are readily detected. In general, 13CO lines trace cooler gas than the bulk 12CO emission and may arise from further out in the disk, as indicated by narrower line profiles. A high fraction of the sources show vibrationally excited emission (~50%) which is correlated with accretion luminosity, consistent with ultra-violet (UV) fluorescent excitation. Disks around early-type Herbig AeBe stars have narrower lines, on average, than their lower-mass late-type counterparts, due to their increased luminosity. Evolutionary changes in CO are also seen. Removal of the protostellar envelope between class I and II results in the disappearance of the strong absorption lines and CO ice feature characteristic of class I spectra. However, CO emission from class I and II objects is similar in detection frequency, excitation and line shape, indicating that inner disk characteristics are established early.

preprint2012arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: III. NGC 1333 IRAS 4A/4B envelope, outflow and UV heating

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

preprint2012arXiv

Detection of the simplest sugar, glycolaldehyde, in a solar-type protostar with ALMA

Glycolaldehyde (HCOCH2OH) is the simplest sugar and an important intermediate in the path toward forming more complex biologically relevant molecules. In this paper we present the first detection of 13 transitions of glycolaldehyde around a solar-type young star, through Atacama Large Millimeter Array (ALMA) observations of the Class 0 protostellar binary IRAS 16293-2422 at 220 GHz (6 transitions) and 690 GHz (7 transitions). The glycolaldehyde lines have their origin in warm (200-300 K) gas close to the individual components of the binary. Glycolaldehyde co-exists with its isomer, methyl formate (HCOOCH3), which is a factor 10-15 more abundant toward the two sources. The data also show a tentative detection of ethylene glycol, the reduced alcohol of glycolaldehyde. In the 690 GHz data, the seven transitions predicted to have the highest optical depths based on modeling of the 220 GHz lines all show red-shifted absorption profiles toward one of the components in the binary (IRAS16293B) indicative of infall and emission at the systemic velocity offset from this by about 0.2" (25 AU). We discuss the constraints on the chemical formation of glycolaldehyde and other organic species - in particular, in the context of laboratory experiments of photochemistry of methanol-containing ices. The relative abundances appear to be consistent with UV photochemistry of a CH3OH-CO mixed ice that has undergone mild heating. The order of magnitude increase in line density in these early ALMA data illustrate its huge potential to reveal the full chemical complexity associated with the formation of solar system analogs.

preprint2012arXiv

Exploring organic chemistry in planet-forming zones

Over the last few years, the chemistry of molecules other than CO in the planet-forming zones of disks is starting to be explored with Spitzer and high-resolution ground-based data. However, these studies have focused only on a few simple molecules. The aim of this study is to put observational constraints on the presence of more complex organic and sulfur-bearing molecules predicted to be abundant in chemical models of disks and to simulate high resolution spectra in view of future missions. High S/N Spitzer spectra at 10-30 micron of the near edge-on disks IRS46 and GVTau are used to search for mid-infrared absorption bands of various molecules. These disks are good laboratories because absorption studies do not suffer from low line/continuum ratios that plague emission data. Simple LTE slab models are used to infer column densities (or upper limits) and excitation temperatures. Bands of HCN, C2H2 and CO2 are clearly detected toward both sources. The HCN and C2H2 absorption arises in warm gas with Tex of 400-700 K, whereas the CO2 absorption originates in cooler gas of app. 250 K (as in Lahuis 2006). No other absorption features are detected. Limits of those molecules are determined and compared with disk models. The inferred ratios wrt. C2H2 and HCN are roughly consistent with models of the chemistry in high-T gas. Models of UV irradiated disk surfaces generally agree better than pure X-ray models. The limit on NH3/HCN implies that evaporation of NH3 containing ices is only a minor contributor. The inferred ratios also compare well with those found in comets, suggesting that part of the cometary material may derive from warm inner disk gas. High resolution simulations show that future instruments on JWST, ELTs, SOFIA and SPICA can probe up to an order of magnitude lower ratios and put important new constraints on the models, especially if pushed to high S/N ratios.

preprint2012arXiv

First detection of water vapor in a pre-stellar core

Water is a crucial molecule in molecular astrophysics as it controls much of the gas/grain chemistry, including the formation and evolution of more complex organic molecules in ices. Pre-stellar cores provide the original reservoir of material from which future planetary systems are built, but few observational constraints exist on the formation of water and its partitioning between gas and ice in the densest cores. Thanks to the high sensitivity of the Herschel Space Observatory, we report on the first detection of water vapor at high spectral resolution toward a dense cloud on the verge of star formation, the pre-stellar core L1544. The line shows an inverse P-Cygni profile, characteristic of gravitational contraction. To reproduce the observations, water vapor has to be present in the cold and dense central few thousand AU of L1544, where species heavier than Helium are expected to freeze-out onto dust grains, and the ortho:para H2 ratio has to be around 1:1 or larger. The observed amount of water vapor within the core (about 1.5x10^{-6} Msun) can be maintained by Far-UV photons locally produced by the impact of galactic cosmic rays with H2 molecules. Such FUV photons irradiate the icy mantles, liberating water wapor in the core center. Our Herschel data, combined with radiative transfer and chemical/dynamical models, shed light on the interplay between gas and solids in dense interstellar clouds and provide the first measurement of the water vapor abundance profile across the parent cloud of a future solar-type star and its potential planetary system.

preprint2012arXiv

Herschel Search for O2 Toward the Orion Bar

We report the results of a search for molecular oxygen (O2) toward the Orion Bar, a prominent photodissociation region at the southern edge of the HII region created by the luminous Trapezium stars. We observed the spectral region around the frequency of the O2 N_J = 3_3 - 1_2 transition at 487 GHz and the 5_4 - 3_4 transition at 774 GHz using the Heterodyne Instrument for the Far Infrared on the Herschel Space Observatory. Neither line was detected, but the 3sigma upper limits established here translate to a total line-of-sight O2 column density < 1.5 10^16 cm^-2 for an emitting region whose temperature is between 30K and 250 K, or < 1 10^16 cm^-2 if the O2 emitting region is primarily at a temperature of ~< 100 K. Because the Orion Bar is oriented nearly edge-on relative to our line of sight, the observed column density is enhanced by a factor estimated to be between 4 and 20 relative to the face-on value. Our upper limits imply that the face-on O2 column density is less than 4 10^15 cm^-2, a value that is below, and possibly well below, model predictions for gas with a density of 10^4 - 10^5 cm^-3 exposed to a far ultraviolet flux 10^4 times the local value, conditions inferred from previous observations of the Orion Bar. The discrepancy might be resolved if: (1) the adsorption energy of O atoms to ice is greater than 800 K; (2) the total face-on Av of the Bar is less than required for O2 to reach peak abundance; (3) the O2 emission arises within dense clumps with a small beam filling factor; or, (4) the face-on depth into the Bar where O2 reaches its peak abundance, which is density dependent, corresponds to a sky position different from that sampled by our Herschel beams.

preprint2012arXiv

Matryoshka Holes: Nested Emission Rings in the Transitional Disk Oph IRS 48

The processes that form transition disks - disks with depleted inner regions - are not well understood; possible scenarios include planet formation, grain growth and photoevaporation. Disks with spatially resolved dust holes are rare, but, in general, even less is known about the gas structure. The disk surrounding A0 star Oph IRS 48 in the nearby Rho Ophiuchus region has a 30 AU radius hole previously detected in the 18.7 micron dust continuum and in warm CO in the 5 micron fundamental ro-vibrational band. We present here Submillimeter Array 880 micron continuum imaging resolving an inner hole. However, the radius of the hole in the millimeter dust is only 13 AU, significantly smaller than measured at other wavelengths. The nesting structure of the disk is counter-intuitive, with increasingly large radii rings of emission seen in the millimeter dust (12.9 +1.7/-3.4 AU), 5 micron CO (30 AU) and 18.7 micron dust (peaking at 55 AU). We discuss possible explanations for this structure, including self-shadowing that cools the disk surface layers, photodissociation of CO, and photoevaporation. However, understanding this unusual disk within the stringent multi-wavelength spatial constraints will require further observations to search for cold atomic and molecular gas.

preprint2012arXiv

Subarcsecond resolution observations of warm water towards three deeply embedded low-mass protostars

Water is present during all stages of star formation: as ice in the cold outer parts of protostellar envelopes and dense inner regions of circumstellar disks, and as gas in the envelopes close to the protostars, in the upper layers of circumstellar disks and in regions of powerful outflows and shocks. In this paper we probe the mechanism regulating the warm gas-phase water abundance in the innermost hundred AU of deeply embedded (Class~0) low-mass protostars, and investigate its chemical relationship to other molecular species during these stages. Millimeter wavelength thermal emission from the para-H2-18O 3(1,3)-2(2,0) (Eu=203.7 K) line is imaged at high angular resolution (0.75"; 190 AU) with the IRAM Plateau de Bure Interferometer towards the deeply embedded low-mass protostars NGC 1333-IRAS2A and NGC 1333-IRAS4A. Compact H2-18O emission is detected towards IRAS2A and one of the components in the IRAS4A binary; in addition CH3OCH3, C2H5CN, and SO2 are detected. Extended water emission is seen towards IRAS2A, possibly associated with the outflow. The detections in all systems suggests that the presence of water on <100 AU scales is a common phenomenon in embedded protostars. We present a scenario in which the origin of the emission from warm water is in a flattened disk-like structure dominated by inward motions rather than rotation. The gas-phase water abundance varies between the sources, but is generally much lower than a canonical abundance of 10^-4, suggesting that most water (>96 %) is frozen out on dust grains at these scales. The derived abundances of CH3OCH3 and SO2 relative to H2-18O are comparable for all sources pointing towards similar chemical processes at work. In contrast, the C2H5CN abundance relative to H2-18O is significantly lower in IRAS2A, which could be due to different chemistry in the sources.

preprint2012arXiv

The Physical Structure of Protoplanetary Disks: the Serpens Cluster Compared with Other Regions

Spectral energy distributions are presented for 94 young stars surrounded by disks in the Serpens Molecular Cloud, based on photometry and Spitzer IRS spectra. Taking a distance to the cloud of 415 pc rather than 259 pc, the distribution of ages is shifted to lower values, in the 1-3 Myr range, with a tail up to 10 Myr. The mass distribution spans 0.2-1.2 Msun, with median mass of 0.7 Msun. The distribution of fractional disk luminosities in Serpens resembles that of the young Taurus Molecular Cloud, with most disks consistent with optically thick, passively irradiated disks in a variety of disk geometries (Ldisk/Lstar ~ 0.1). In contrast, the distributions for the older Upper Scorpius and Eta Chamaeleontis clusters are dominated by optically thin lower luminosity disks (Ldisk/Lstar ~ 0.02). This evolution in fractional disk luminosities is concurrent with that of disk fractions. The actively accreting and non-accreting stars (based on Ha data) in Serpens show very similar distributions in fractional disk luminosities, differing only in the brighter tail dominated by strongly accreting stars. In contrast with a sample of Herbig Ae/Be stars, the T Tauri stars in Serpens do not have a clear separation in fractional disk luminosities for different disk geometries: both flared and flat disks present wider, overlapping distributions. This result is consistent with previous suggestions of a faster evolution for disks around Herbig Ae/Be stars. Furthermore, the results for the mineralogy of the dust in the disk surface do not show any correlation to either stellar and disk characteristics or mean cluster age in the 1-10 Myr range probed here. A possible explanation for the lack of correlation is that the processes affecting the dust within disks have short timescales, happening repeatedly, making it difficult to distinguish long lasting evolutionary effects. [abridged]

preprint2012arXiv

The warm gas atmosphere of the HD 100546 disk seen by Herschel (Evidence of a gas-rich, carbon-poor atmosphere?)

(Abridged) Context. With the Herschel Space Observatory, lines of simple molecules (C+, O, and CO) have been observed in the atmosphere of protoplanetary disks. When combined with ground-based [CI], all principle forms of carbon can be studied. The absence of neutral carbon [CI], which is predicted by models to be strong, can then be interpreted together with [CII] and carbon monoxide. Aims. We study the gas temperature, excitation, and chemical abundance of the simple carbon-bearing species by the method of chemical-physical modeling. We explore the sensitivity of the lines to the entering parameters and constrain the region from which the line radiation emerges. Methods. Numerical models of the radiative transfer are used together with a chemical network simulation and a calculation of the gas energetics to obtain the gas temperature. We present our new model, which is based on our previous models but includes several improvements. Results. A model of the disk around the Herbig Be star HD 100546 is able to reproduce the CO ladder together with the atomic fine-structure lines of [OI] and either [CI] or [CII]. We find that the high-J lines of CO can only be reproduced by a warm atmosphere with Tgas>>Tdust. The high-J CO observable with PACS are dominated from regions within some tens of AU. Conclusions. Only a warm atmosphere with Tgas>>Tdust can reproduce the CO ladder. The CO ladder together with [O I] and the upper limit to [CI] can be reproduced by models with a high gas/dust ratio and a low abundance of volatile carbon. These models however produce too small amounts of [CII]. Models with a low gas/dust ratio and more volatile carbon also reproduce CO and [OI], are in closer agreement with observations of [CII], but overproduce [CI]. Due to the uncertain origin of the [CII] emission, we prefer the high gas/dust ratio models, indicating a low abundance of volatile carbon.

preprint2012arXiv

Warm water deuterium fractionation in IRAS 16293-2422 - The high-resolution ALMA and SMA view

Measuring the water deuterium fractionation in the inner warm regions of low-mass protostars has so far been hampered by poor angular resolution obtainable with single-dish ground- and space-based telescopes. Observations of water isotopologues using (sub)millimeter wavelength interferometers have the potential to shed light on this matter. Observations toward IRAS 16293-2422 of the 5(3,2)-4(4,1) transition of H2-18O at 692.07914 GHz from Atacama Large Millimeter/submillimeter Array (ALMA) as well as the 3(1,3)-2(2,0) of H2-18O at 203.40752 GHz and the 3(1,2)-2(2,1) transition of HDO at 225.89672 GHz from the Submillimeter Array (SMA) are presented. The 692 GHz H2-18O line is seen toward both components of the binary protostar. Toward one of the components, "source B", the line is seen in absorption toward the continuum, slightly red-shifted from the systemic velocity, whereas emission is seen off-source at the systemic velocity. Toward the other component, "source A", the two HDO and H2-18O lines are detected as well with the SMA. From the H2-18O transitions the excitation temperature is estimated at 124 +/- 12 K. The calculated HDO/H2O ratio is (9.2 +/- 2.6)*10^(-4) - significantly lower than previous estimates in the warm gas close to the source. It is also lower by a factor of ~5 than the ratio deduced in the outer envelope. Our observations reveal the physical and chemical structure of water vapor close to the protostars on solar-system scales. The red-shifted absorption detected toward source B is indicative of infall. The excitation temperature is consistent with the picture of water ice evaporation close to the protostar. The low HDO/H2O ratio deduced here suggests that the differences between the inner regions of the protostars and the Earth's oceans and comets are smaller than previously thought.

preprint2012arXiv

Water in star-forming regions with Herschel: highly excited molecular emission from the NGC 1333 IRAS 4B outflow

During the embedded phase of pre-main sequence stellar evolution, a disk forms from the dense envelope while an accretion-driven outflow carves out a cavity within the envelope. Highly excited H2O emission in spatially unresolved Spitzer/IRS spectra of a low-mass Class 0 object, NGC 1333 IRAS 4B, has previously been attributed to the envelope-disk accretion shock but could instead be produced in an outflow. As part of the survey of low-mass sources in the Water in Star Forming Regions with Herschel (WISH-LM) program, we used Herschel/PACS to obtain a far-IR spectrum and several Nyquist-sampled spectral images with to determine the origin of excited H2O emission from NGC 1333 IRAS 4B. The spectrum has high signal-to-noise in a rich forest of H2O, CO, and OH lines, providing a near-complete census of far-IR molecular emission from a Class 0 protostar. The excitation diagrams for the three molecules all require fits with two excitation temperatures, indicating the presence of two physical components. The highly excited component of H2O emission is characterized by subthermal excitation of 1500 K gas with a density of 10^5 - 10^7 cm-3, conditions that also reproduce the mid-IR H2O emission detected by Spitzer. On the other hand, a high density, low temperature gas can reproduce the H2O spectrum observed by Spitzer but underpredicts the H2O lines seen by Herschel. Nyquist-sampled spectral maps of several lines show two spatial components of H2O emission, one centered at 1200 AU south of the central source at the position of the blueshifted outflow lobe and a second centered on-source. Both spatial components of the far-IR H2O emission are consistent with emission from the outflow. The gas cooling from the IRAS 4B envelope cavity walls is dominated by far-IR H2O emission, in contrast to stronger [O I] and CO cooling from more evolved protostars. [one sentence truncated]

preprint2011arXiv

A 30 AU radius CO gas hole in the disk around the Herbig Ae star Oph IRS 48

The physical processes leading to the disappearance of disks around young stars are not well understood. A subclass of transitional disks, the so-called cold disks with large inner dust holes, provide a crucial laboratory for studying disk dissipation processes. IRS 48 has a 30 AU radius hole previously measured from dust continuum imaging at 18.7 micron. Using new optical spectra, we determine that IRS 48 is a pre-main sequence A0 star. In order to characterize this disk's gas distribution, we obtained AO-assisted VLT CRIRES high resolution (R ~100,000) spectra of the CO fundamental rovibrational band at 4.7 micron. All CO emission, including that from isotopologues and vibrationally excited molecules, is off-source and peaks at 30 AU. The gas is thermally excited to a rotational temperature of 260 K and is also strongly UV pumped, showing a vibrational excitation temperature of ~5000 K. We model the kinematics and excitation of the gas and posit that the CO emission arises from the dust hole wall. Prior imaging of UV-excited PAH molecules, usually a gas tracer, within the hole makes the large CO hole even more unexpected.

preprint2011arXiv

Chemical history of molecules in circumstellar disks

The chemical composition of a protoplanetary disk is determined not only by in situ chemical processes during the disk phase, but also by the history of the gas and dust before it accreted from the natal envelope. In order to understand the disk's chemical composition at the time of planet formation, especially in the midplane, one has to go back in time and retrace the chemistry to the molecular cloud that collapsed to form the disk and the central star. Here we present a new astrochemical model that aims to do just that. The model follows the core collapse and disk formation in two dimensions, which turns out to be a critical upgrade over older collapse models. We predict chemical stratification in the disk due to different physical conditions encountered along different streamlines. We argue that the disk-envelope accretion shock does not play a significant role for the material in the disk at the end of the collapse phase. Finally, our model suggests that complex organic species are formed on the grain surfaces at temperatures of 20 to 40 K, rather than in the gas phase in the T>100 K hot corino.

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

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.

preprint2011arXiv

Disks and Outflows in CO Rovibrational Emission from Embedded, Low-Mass Young Stellar Objects

Young circumstellar disks that are still embedded in dense molecular envelopes may differ from their older counterparts, but are historically difficult to study because emission from a disk can be confused with envelope or outflow emission. CO fundamental emission is a potentially powerful probe of the disk/wind structure within a few AU of young protostars. In this paper, we present high spectral (R=90,000) and spatial (0.3") resolution VLT/CRIRES M-band spectra of 18 low-mass young stellar objects (YSOs) with dense envelopes in nearby star-froming regions to explore the utility of CO fundamental 4.6 micron emission as a probe of very young disks. CO fundamental emission is detected from 14 of the YSOs in our sample. The emission line profiles show a range of strengths and shapes, but can generally be classified into a broad, warm component and a narrow, cool component. The broad CO emission is detected more frequently from YSOs with bolometric luminosities of <15 Lsun than those with >15 Lsun, and as with CO emission from CTTSs is attributed to the warm (~1000 K) inner AU of the disk. The CO emission from objects with high bolometric luminosity is produced in cooler (~320 K), narrow lines in 12CO and in rarer isotopologues. From some objects, the narrow lines are blueshifted by up to ~10 km/s, indicating a slow wind origin. For other sources the lines are located at the systemic velocity of the star and likely arise in the disk. For a few YSOs, spatially-extended CO and H2 S(9) emission is detected up to 2" from the central source and is attributed to interactions between the wind and surrounding molecular material. Warm CO absorption is detected in the wind of six objects with velocities up to 100 km/s, often in discrete velocity components. That the wind is partially molecular where it is launched favors ejection in a disk wind rather than a coronal or chromospheric wind.

preprint2011arXiv

Energetic processes revealed by spectrally resolved high-J CO lines in low-mass star-forming regions with Herschel-HIFI

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars in NGC1333. The observations show several energetic components including shocked and quiescent gas. Radiative transfer models are used to quantify the C18O envelope abundance which require a jump in the abundance at an evaporation temperature, T_ev ~25 K, providing new direct evidence of a CO ice evaporation zone around protostars. The abundance in the outermost part of the envelope, X_0, is within the canonical value of 2x10^-4; however the inner abundance, X_in, is found around a factor of 3-5 lower than X_0.

preprint2011arXiv

Herschel Measurements of Molecular Oxygen in Orion

We report observations of three rotational transitions of molecular oxygen (O2) in emission from the H2 Peak 1 position of vibrationally excited molecular hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines using HIFI on the Herschel Space Observatory, having velocities of 11 km s-1 to 12 km s-1 and widths of 3 km s-1. The beam-averaged column density is N(O2) = 6.5\times1016 cm-2, and assuming that the source has an equal beam filling factor for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 relative to H2 is 0.3 - 7.3\times10-6. The unusual velocity suggests an association with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. The mass of this source is ~ 10 M\odot and the dust temperature is \geq 150 K. Our preferred explanation of the enhanced O2 abundance is that dust grains in this region are sufficiently warm (T \geq 100 K) to desorb water ice and thus keep a significant fraction of elemental oxygen in the gas phase, with a significant fraction as O2. For this small source, the line ratios require a temperature \geq 180 K. The inferred O2 column density \simeq 5\times1018 cm-2 can be produced in Peak A, having N(H2) \simeq 4\times1024 cm-2. An alternative mechanism is a low-velocity (10 to 15 km s-1) C-shock, which can produce N(O2) up to 1017 cm-2.

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

On the Evolution of Dust Mineralogy, From Protoplanetary Disks to Planetary Systems

Mineralogical studies of silicate features emitted by dust grains in protoplanetary disks and Solar System bodies can shed light on the progress of planet formation. The significant fraction of crystalline material in comets, chondritic meteorites and interplanetary dust particles indicates a modification of the almost completely amorphous ISM dust from which they formed. The production of crystalline silicates thus must happen in protoplanetary disks, where dust evolves to build planets and planetesimals. Different scenarios have been proposed, but it is still unclear how and when this happens. This paper presents dust grain mineralogy of a complete sample of protoplanetary disks in the young Serpens cluster. These results are compared to those in the young Taurus region and to sources that have retained their protoplanetary disks in the older Upper Scorpius and Eta Chamaeleontis stellar clusters, using the same analysis technique for all samples. This comparison allows an investigation of the grain mineralogy evolution with time for a total sample of 139 disks. The mean cluster age and disk fraction are used as indicators of the evolutionary stage of the different populations. Our results show that the disks in the different regions have similar distributions of mean grain sizes and crystallinity fractions (~10-20%) despite the spread in mean ages. Furthermore, there is no evidence of preferential grain sizes for any given disk geometry, nor for the mean cluster crystallinity fraction to increase with mean age in the 1-8 Myr range. The main implication is that a modest level of crystallinity is established in the disk surface early on (< 1 Myr), reaching a equilibrium that is independent of what may be happening in the disk midplane. These results are discussed in the context of planet formation, in comparison with mineralogical results from small bodies in our Solar System. [Abridged]

preprint2011arXiv

Spectroscopic properties of Young Stellar Objects in the Lupus Molecular Clouds

The results of an optical spectroscopic survey of a sample of young stellar objects (YSOs) and pre-main sequence (PMS) stars in the Lupus Clouds are presented. 92 objects were observed with VLT/FLAMES. All of those objects show IR excess as discovered by the Spitzer Legacy Program "From Molecular Cores to Planet-Forming Disks" (c2d). After reduction, 54 spectra with good signal-to-noise ratio are spectrally classified. Effective temperatures and luminosities are derived for these objects, and used to construct H-R diagrams for the population. The sample consists mostly of M-type stars, with 10% K-type stars. Individual ages and masses are inferred for the objects according to theoretical evolutionary models. The mean population age is found to be between 3.6 and 4.4 Myr, depending on the model, while the mean mass is found to be ~0.3 M for either model. Together with literature data, the distribution of spectral types is found to be similar to that in Chamaeleon I and IC348. The Hα line in emission, found in 49% of the sample, is used to distinguish between classical and weak-line T Tauri stars. 56% of the objects show Hα in emission and are accreting T Tauri stars. Mass accretion rates between 10-8 and 10-11 M yr-1 are determined from the full width at 10% of the Hα peak intensity. These mass accretion rates are, within a large scatter, consistent with the Mac \propto M 2 relation found in the literature.

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

Tracing FUV Radiation in the Embedded Phase of Star Formation

Molecules containing one or a few hydrogen atoms and a heavier atom (hydrides) have been predicted to trace FUV radiation. In some chemical models, FUV emission by the central object or protostar of a star forming region greatly enhances some of the hydride abundances. Two massive regions, W3 IRS5 and AFGL 2591, have been observed in hydride lines by HIFI onboard the {\it Herschel Space Observatory}. We use published results as well as new observations of CH$^+$ towards W3 IRS5. Molecular column densities are derived from ground state absorption lines, radiative transfer modeling or rotational diagrams. Models assuming no internal FUV are compared with two-dimensional models including FUV irradiation of outflow walls. We confirm that the effect of FUV is clearly noticeable and greatly improves the fit. The most sensitive molecules to FUV irradiation are CH$^+$ and OH$^+$, enhanced in abundance by many orders of magnitude. Modeling in addition also full line radiative transfer, Bruderer et al (2010b) achieve good agreement of a two-dimensional FUV model with observations of CH$^+$ in AFGL 2591. It is concluded that CH$^+$ and OH$^+$ are good FUV tracers in star-forming regions.

preprint2011arXiv

Water in Star and Planet Forming Regions

In this paper we discuss the astronomical search for water vapor in order to understand the disposition of water in all its phases throughout the process of star and planet formation. Our ability to detect and study water vapor has recently received a tremendous boost with the successful launch and operations of the Herschel Space Observatory. Herschel spectroscopic detections of numerous transitions in a variety of astronomical objects, along with previous work by other space-based observatories, will be threaded throughout this contribution. In particular, we present observations of water tracing the earliest stage of star birth where it is predominantly frozen as ice. When a star is born the local energy release by radiation liberates ices in its surrounding envelope and powers energetic outflows that appear to be water factories. In these regions water plays an important role in the gas physics. Finally, we end with an exploration of water in planet forming disks surrounding young stars. The availability of accurate molecular data (frequencies, collisional rate coefficients, and chemical reaction rates) are crucial to analyze the observations at each of these steps.

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

A Spitzer c2d Legacy Survey to Identify and Characterize Disks with Inner Dust Holes

Understanding how disks dissipate is essential to studies of planet formation. However, identifying exactly how dust and gas dissipates is complicated due to difficulty in finding objects clearly in the transition of losing their surrounding material. We use Spitzer IRS spectra to examine 35 photometrically-selected candidate cold disks (disks with large inner dust holes). The infrared spectra are supplemented with optical spectra to determine stellar and accretion properties and 1.3mm photometry to measure disk masses. Based on detailed SED modeling, we identify 15 new cold disks. The remaining 20 objects have IRS spectra that are consistent with disks without holes, disks that are observed close to edge-on, or stars with background emission. Based on these results, we determine reliable criteria for identifying disks with inner holes from Spitzer photometry and examine criteria already in the literature. Applying these criteria to the c2d surveyed star-forming regions gives a frequency of such objects of at least 4% and most likely of order 12% of the YSO population identified by Spitzer. We also examine the properties of these new cold disks in combination with cold disks from the literature. Hole sizes in this sample are generally smaller than for previously discovered disks and reflect a distribution in better agreement with exoplanet orbit radii. We find correlations between hole size and both disk and stellar masses. Silicate features, including crystalline features, are present in the overwhelming majority of the sample although 10 micron feature strength above the continuum declines for holes with radii larger than ~7 AU. In contrast, PAHs are only detected in 2 out of 15 sources. Only a quarter of the cold disk sample shows no signs of accretion, making it unlikely that photoevaporation is the dominant hole forming process in most cases.

preprint2010arXiv

A Spitzer Survey of Protoplanetary Disk Dust in the Young Serpens Cloud: How do Dust Characteristics Evolve with Time?

We present Spitzer IRS mid-infrared (5-35 micron) spectra of a complete flux-limited sample (> 3 mJy at 8 micron) of young stellar object (YSO) candidates selected on the basis of their infrared colors in the Serpens Molecular Cloud. Spectra of 147 sources are presented and classified. Background stars (with slope consistent with a reddened stellar spectrum and silicate features in absorption), galaxies (with redshifted PAH features) and a planetary nebula (with high ionization lines) amount to 22% of contamination in this sample, leaving 115 true YSOs. Sources with rising spectra and ice absorption features, classified as embedded Stage I protostars, amount to 18% of the sample. The remaining 82% (94) of the disk sources are analyzed in terms of spectral energy distribution shapes, PAHs and silicate features. The presence, strength and shape of these silicate features are used to infer disk properties for these systems. About 8% of the disks have 30/13 micron flux ratios consistent with cold disks with inner holes or gaps, and 3% of the disks show PAH emission. Comparison with models indicates that dust grains in the surface of these disks have sizes of at least a few μ\m. The 20 micron silicate feature is sometimes seen in absence of the 10 micron feature, which may be indicative of very small holes in these disks. No significant difference is found in the distribution of silicate feature shapes and strengths between sources in clusters and in the field. Moreover, the results in Serpens are compared with other well-studied samples: the c2d IRS sample distributed over 5 clouds and a large sample of disks in the Taurus star-forming region. The remarkably similar distributions of silicate feature characteristics in samples with different environment and median ages - if significant - imply that the dust population in the disk surface results from an equilibrium between dust growth and destructive collision processes that are maintained over a few million years for any YSO population irrespective of environment.

preprint2010arXiv

c2d Spitzer IRS spectra of embedded low-mass young stars: gas-phase emission lines

A survey of mid-IR gas-phase emission lines of H2, H2O and various atoms toward a sample of 43 embedded low-mass young stars in nearby star-forming regions is presented. The sources are selected from the Spitzer "Cores to Disks" (c2d) legacy program. The environment of embedded protostars is complex both in its physical structure (envelopes, outflows, jets, protostellar disks) and the physical processes (accretion, irradiation by UV and/or X-rays, excitation through slow and fast shocks) which take place. A key point is to spatially resolve the emission in the Spitzer-IRS spectra. An optimal extraction method is used to separate both spatially unresolved (compact, up to a few 100 AU) and spatially resolved (extended, 1000 AU or more) emission from the IRS spectra. The results are compared with the c2d disk sample and literature PDR and shock models to address the physical nature of the sources. Both compact and extended emission features are observed. Warm (Tex few 100 K) H2, observed through the pure rotational H2 S(0), S(1) and S(2) lines, and [S I] 25 mu emission is observed primarily in the extended component. [S I] is observed uniquely toward truly embedded sources and not toward disks. On the other hand hot (Tex>700 K) H2, observed primarily through the S(4) line, and [Ne II] emission is seen mostly in the compact component. [Fe II] and [Si II] lines are observed in both spatial components. Hot H2O emission is found in the compact component of some sources. The observed emission on >=1000 AU scales is characteristic of PDR emission and likely originates in the outflow cavities in the remnant envelope created by the stellar wind and jets from the embedded young stars. Weak shocks along the outflow wall can also contribute. The compact emission is likely of mixed origin, comprised of optically thick circumstellar disk and/or jet/outflow emission from the protostellar object.

preprint2010arXiv

Grain growth across protoplanetary discs: 10-micron silicate feature versus millimetre slope

Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at micron wavelengths are consistent with growth to micron sizes whereas the slope of the SED at longer wavelengths traces growth up to mm sizes. We here look for a correlation between these two grain growth indicators. A large sample of T-Tauri and Herbig-Ae/Be stars was observed with the Spitzer Space Telescope at 5-13 micron; a subsample was observed at mm wavelengths. We complement this subsample with data from the literature to maximise the overlap between micron and mm observations and search for correlations. Synthetic spectra are produced to determine which processes may produce the dust evolution. Dust disc masses in the range <1 to 7 x 10^-4 MSun are obtained. Most sources have a mm spectral slope consistent with grain growth. There is a tentative correlation between the 10-micron silicate feature and the mm slope of the SED. The observed sources seem to be grouped per star-forming region in the micron-vs-mm diagram. The modelling results show that the 10-micron feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects the 10-micron feature, but hardly the mm slope. The tentative correlation between the strength of the 10-micron feature and the mm slope suggests that the inner and outer disc evolve simultaneously. Dust with a mass dominated by mm-sized grains is required to explain the shallowest mm slopes. Other processes besides grain growth may also be responsible for the removal of small grains.

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.

preprint2010arXiv

The HDO/H2O ratio in gas in the inner regions of a low-mass protostar

The HDO/H2O abundance ratio is thought to be a key diagnostic on the evolution of water during the star- and planet-formation process and thus on its origin on Earth. We here present millimeter-wavelength high angular resolution observations of the deeply embedded protostar NGC1333-IRAS4B from the Submillimeter Array targeting the 3(1,2)-2(2,1) transition of HDO at 225.6 GHz (Eu = 170 K). We do not (or only very tentatively) detect the HDO line toward the central protostar, contrasting the previous prominent detection of a line from another water isotopologue, H2-18-O, with similar excitation properties using the IRAM Plateau de Bure Interferometer. The non-detection of the HDO line provides a direct, model independent, upper limit to the HDO/H2O abundance ratio of 6e(-4) (3sigma) in the warm gas associated with the central protostar. This upper limit suggests that the HDO/H2O abundance ratio is not significantly enhanced in the inner 50 AU around the protostar relative to what is seen in comets and Earth's oceans and does not support previous suggestions of a generally enhanced HDO/H2O ratio in these systems.

preprint2010arXiv

The Spitzer c2d Survey of Weak-Line T Tauri Stars. III. The Transition from Primordial Disks to Debris Disks

We present 3.6 to 70 μm Spitzer photometry of 154 weak-line T Tauri stars (WTTS) in the Chamaeleon, Lupus, Ophiuchus and Taurus star formation regions, all of which are within 200 pc of the Sun. For a comparative study, we also include 33 classical T Tauri stars (CTTS) which are located in the same star forming regions. Spitzer sensitivities allow us to robustly detect the photosphere in the IRAC bands (3.6 to 8 μm) and the 24 μm MIPS band. In the 70 μm MIPS band, we are able to detect dust emission brighter than roughly 40 times the photosphere. These observations represent the most sensitive WTTS survey in the mid to far infrared to date, and reveal the frequency of outer disks (r = 3-50 AU) around WTTS. The 70 μm photometry for half the c2d WTTS sample (the on-cloud objects), which were not included in the earlier papers in this series, Padgett et al. (2006) and Cieza et al. (2007), are presented here for the first time. We find a disk frequency of 19% for on-cloud WTTS, but just 5% for off- cloud WTTS, similar to the value reported in the earlier works. WTTS exhibit spectral energy distributions (SEDs) that are quite diverse, spanning the range from optically thick to optically thin disks. Most disks become more tenuous than Ldisk/L* = 2 x 10^-3 in 2 Myr, and more tenuous than Ldisk/L* = 5 x 10^-4 in 4 Myr.

preprint2010arXiv

VLT/X-shooter Spectroscopy of a Dusty Planetary Nebula Discovered with Spitzer/IRS

As part of a mid-infrared spectroscopic survey of young stars with the Spitzer Space Telescope, an unclassified red emission line object was discovered. Based on its high ionization state indicated by the Spitzer spectrum, this object could either be a dusty Supernova Remnant (SNR) or a Planetary Nebula (PN). In this research note, the object is classified and the available spectroscopic data are presented to the community for further analysis. UV/optical/NIR spectra were obtained during the science verification run of the VLT/X-shooter. A large number of emission lines are identified allowing the determination of the nature of this object. The presence of strong, narrow (Delta_v ~ 8 - 74 km/s) emission lines, combined with very low line ratios of, e.g., [N II]/Halpha and [S II]/Halpha show that the object is a Planetary Nebula (PN) that lies at an undetermined distance behind the Serpens Molecular Cloud. This illustrates the potential of X-shooter as an efficient tool for constraining the nature of faint sources with unknown spectral properties or colors.

preprint2010arXiv

Water Vapor in the Inner 25 AU of a Young Disk around a Low-Mass Protostar

Water is one of the key molecules in the physical and chemical evolution of star- and planet-forming regions. We here report the first spatially resolved observation of thermal emission of (an isotopologue of) water with the Plateau de Bure Interferometer toward the deeply embedded Class 0 protostar NGC 1333-IRAS4B. The observations of the H2-18-O 3_13-2_20 transition at 203.4 GHz resolve the emission of water toward this source with an extent of about 0.2" corresponding to the inner 25 AU (radius). The H2-18-O emission reveals a tentative velocity gradient perpendicular to the extent of the protostellar outflow/jet probed by observations of CO rotational transitions and water masers. The line is narrow, about 1 km/s (FWHM), significantly less than what would be expected for emission from an infalling envelope or accretion shock, but consistent with emission from a disk seen at a low inclination angle. The water column density inferred from these data suggests that the water emitting gas is a thin warm layer containing about 25 M_Earth of material, 0.03% of the total disk mass traced by continuum observations.

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

IRC+10216's Innermost Envelope -- The eSMA's View

We used the Extended Submillimeter Array (eSMA) in its most extended configuration to investigate the innermost (within a radius of 290 R* from the star) circumstellar envelope (CSE) of IRC+10216. We imaged the CSE using HCN and other molecular lines with a beam size of 0."22 x 0."46, deeply into the very inner edge (15 R*) of the envelope where the expansion velocity is only 3 km/s. The excitation mechanism of hot HCN and KCl maser lines is discussed. HCN maser components are spatially resolved for the first time on an astronomical object. We identified two discrete regions in the envelope: a region with a radius of . 15 R*, where molecular species have just formed and the gas has begun to be accelerated (region I) and a shell region (region II) with a radius of 23 R* and a thickness of 15 R*, whose expansion velocity has reached up to 13 km/s, nearly the terminal velocity of 15 km/s. The Si$^{34}$S line detected in region I shows a large expansion velocity of 16 km/s due to strong wing components, indicating that the emission may arise from a shock region in the innermost envelope. In region II, the P.A. of the most copious mass loss direction was found to be 120 +/- 10 degrees, which may correspond to the equatorial direction of the star. Region II contains a torus-like feature. These two regions may have emerged due to significant differences in the size distributions of the dust particles in the two regions.

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.

preprint2009arXiv

PROSAC: A Submillimeter Array survey of low-mass protostars. II. The mass evolution of envelopes, disks, and stars from the Class 0 through I stages

The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through the circumstellar disk onto the central star. A sample of 20 Class 0 and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular resolution with the Submillimeter Array. Using detailed dust radiative transfer models, we have developed a framework for disentangling the continuum emission from the envelopes and disks, and from that estimated their masses. For the Class I sources in the sample, HCO+ 3-2 line emission has furthermore been observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, constraining the masses of the central stars. Both Class 0 and I protostars are surrounded by disks with typical masses of about 0.05 M_sun. No evidence is found for a correlation between the disk mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from 1 M_sun in the Class 0 stage to <0.1 M_sun in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses in the Class 0 stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the embedded protostellar stages. Theoretical models tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution and rapidly process material accreted from the larger scale envelope onto the central star.

preprint2009arXiv

Quantification of segregation dynamics in ice mixtures

(Abridged) The observed presence of pure CO2 ice in protostellar envelopes is attributed to thermally induced ice segregation, but a lack of quantitative experimental data has prevented its use as a temperature probe. Quantitative segregation studies are also needed to characterize diffusion in ices, which underpins all ice dynamics and ice chemistry. This study aims to quantify the segregation mechanism and barriers in different H2O:CO2 and H2O:CO ice mixtures covering a range of astrophysically relevant ice thicknesses and mixture ratios. The ices are deposited at 16-50 K under (ultra-)high vacuum conditions. Segregation is then monitored at 23-70 K as a function of time, through infrared spectroscopy. Thin (8-37 ML) H2O:CO2/CO ice mixtures segregate sequentially through surface processes, followed by an order of magnitude slower bulk diffusion. Thicker ices (>100 ML) segregate through a fast bulk process. The thick ices must therefore be either more porous or segregate through a different mechanism, e.g. a phase transition. The segregation dynamics of thin ices are reproduced qualitatively in Monte Carlo simulations of surface hopping and pair swapping. The experimentally determined surface-segregation rates for all mixture ratios follow the Ahrrenius law with a barrier of 1080[190] K for H2O:CO2 and 300[100] K for H2O:CO mixtures. During low-mass star formation H2O:CO2 segregation will be important already at 30[5] K. Both surface and bulk segregation is proposed to be a general feature of ice mixtures when the average bond strengths of the mixture constituents in pure ice exceeds the average bond strength in the ice mixture.

preprint2008arXiv

Detection of CI in absorption towards PKS 1830-211 with the eSMA

We report the first science observations and results obtained with the "extended" SMA (eSMA), which is composed of the SMA (Submillimeter Array), JCMT (James Clerk Maxwell Telescope) and CSO (Caltech Submillimeter Observatory). Redshifted absorptions at z=0.886 of CI (^3P_1 - ^3P_0) were observed with the eSMA with an angular resolution of 0.55"x0.22" at 1.1 mm toward the southwestern image of the remarkable lensed quasar PKS 1830-211, but not toward the northeastern component at a separation of ~1". Additionally, SMA observations of CO, 13CO and C18O (all J=4-3) were obtained toward this object: CO was also detected toward the SW component, but none of the isotopologues were. This is the first time [CI] is detected in this object, allowing the first direct determination of relative abundances of neutral atomic carbon to CO in the molecular clouds of a spiral galaxy at z>0.1. The [CI] and CO profiles can be decomposed into two and three velocity components respectively. We derive C/CO column density ratios ranging from <0.5 (representative of dense cores) to ~2.5 (close to translucent clouds values). This could indicate that we are seeing environments with different physical conditions or that we are witnessing chemical evolution of regions where C has not completely been converted into CO.

preprint2008arXiv

Optical Characterization of a New Young Stellar Population in the Serpens Molecular Cloud

We report on the results of an optical spectroscopic survey designed to confirm the youth and determine the spectral types among a sample of young stellar object (YSO) candidates in the Serpens Molecular Cloud. We observed 150 infrared excess objects, previously discovered by the Spitzer Legacy Program "From Molecular Cores to Planet-Forming Disks" (c2d), bright enough for subsequent Spitzer/IRS spectroscopy. We obtained 78 optical spectra of sufficient S/N for analysis. Extinctions, effective temperatures and luminosities are estimated for this sample, and used to construct H-R diagrams for the population. We identified 20 background giants contaminating the sample, based on their relatively high extinction, position in the H-R diagram, the lack of Halpha emission and relatively low infrared excess. Such strong background contamination (25%) is consistent with the location of Serpens being close to the Galactic plane (5degrees Galactic latitude). The remaining 58 stars (75%) were all confirmed to be young, mostly K and M-type stars that are presumed to belong to the cloud. Individual ages and masses for the YSOs are inferred based on theoretical evolutionary models. The models indicate a spread in stellar ages from 1 to 15 Myr, peaking at 2 - 6 Myr, and a mass distribution of 0.2 to 1.2 Msun with median value around 0.8 Msun. Strong H emission lines (EW[Halpha] > 3 A) have been detected in more than half of the sample (35 stars). The mass accretion rates as derived from the H line widths span a broad distribution over 4 orders of magnitude with median accretion rate of 10^-8 Msun/yr. Our analysis shows that the majority of the infrared excess objects detected in Serpens are actively accreting, young T-Tauri stars.

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.

preprint2008arXiv

The Spitzer c2d Legacy Results: Star Formation Rates and Efficiencies; Evolution and Lifetimes

(Abridged) The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instruments for five large, nearby molecular clouds. This paper combines information drawn from studies of individual clouds into a combined and updated statistical analysis of star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. Current star formation efficiencies range from 3% to 6%. Taken together, the five clouds are producing about 260 solar masses of stars per Myr. The star formation surface density is more than an order of magnitude larger than would be predicted from the Kennicutt relation used in extragalactic studies. Measured against the dense gas probed by the maps of dust continuum emission, the efficiencies are much higher, and the current stock of dense cores would be exhausted in 1.8 Myr on average. The derived lifetime for the Class I phase is 0.44 to 0.54 Myr, considerably longer than some estimates. Similarly, the lifetime for the Class 0 SED class, 0.10 to 0.16 Myr, is longer than early estimates. The great majority (90%) of young stars lie within loose clusters with at least 35 members and a stellar density of 1 solar mass per cubic pc. Accretion at the sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star of about 0.25 solar masses, given an efficiency of 0.3. Our data confirm and aggravate the "luminosity problem" for protostars. Our results strongly suggest that accretion is time variable, with prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight most luminous objects.

preprint2007arXiv

PROSAC: A Submillimeter Array Survey of Low-Mass Protostars. I. Overview of Program: Envelopes, Disks, Outflows and Hot Cores

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

preprint2007arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds VIII. Serpens Observed with MIPS

We present maps of 1.5 square degrees of the Serpens dark cloud at 24, 70, and 160\micron observed with the Spitzer Space Telescope MIPS Camera. More than 2400 compact sources have been extracted at 24um, nearly 100 at 70um, and 4 at 160um. We estimate completeness limits for our 24um survey from Monte Carlo tests with artificial sources inserted into the Spitzer maps. We compare source counts, colors, and magnitudes in the Serpens cloud to two reference data sets, a 0.50 deg^2 set on a low-extinction region near the dark cloud, and a 5.3 deg^2 subset of the SWIRE ELAIS N1 data that was processed through our pipeline. These results show that there is an easily identifiable population of young stellar object candidates in the Serpens Cloud that is not present in either of the reference data sets. We also show a comparison of visual extinction and cool dust emission illustrating a close correlation between the two, and find that the most embedded YSO candidates are located in the areas of highest visual extinction.

preprint2007arXiv

The Spitzer c2d Survey of Weak-line T Tauri Stars II: New Constraints on the Timescale for Planet Building

One of the central goals of the Spitzer Legacy Project ``From Molecular Cores to Planet-forming Disks'' (c2d) is to determine the frequency of remnant circumstellar disks around weak-line T Tauri stars (wTTs) and to study the properties and evolutionary status of these disks. Here we present a census of disks for a sample of over 230 spectroscopically identified wTTs located in the c2d IRAC (3.6, 4.5, 4.8, and 8.0 um) and MIPS (24 um) maps of the Ophiuchus, Lupus, and Perseus Molecular Clouds. We find that ~20% of the wTTs in a magnitude limited subsample have noticeable IR-excesses at IRAC wavelengths indicating the presence of a circumstellar disk. The disk frequencies we find in these 3 regions are ~3-6 times larger than that recently found for a sample of 83 relatively isolated wTTs located, for the most part, outside the highest extinction regions covered by the c2d IRAC and MIPS maps. The disk fractions we find are more consistent with those obtained in recent Spitzer studies of wTTs in young clusters such as IC 348 and Tr 37. From their location in the H-R diagram, we find that, in our sample, the wTTs with excesses are among the younger part of the age distribution. Still, up to ~50% of the apparently youngest stars in the sample show no evidence of IR excess, suggesting that the circumstellar disks of a sizable fraction of pre-main-sequence stars dissipate in a timescale of ~1 Myr. We also find that none of the stars in our sample apparently older than ~10 Myrs have detectable circumstellar disks at wavelengths < 24 um. Also, we find that the wTTs disks in our sample exhibit a wide range of properties (SED morphology, inner radius, L_DISK/L*, etc) which bridge the gaps observed between the cTTs and the debris disk regimes.

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.

preprint2006arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. III. Perseus Observed with IRAC

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

preprint2006arXiv

The Spitzer c2d survey of Large, Nearby, Interstellar Clouds. V. Chamaeleon II Observed with IRAC

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

preprint2005arXiv

Photodissociation of water in crystalline ice: a molecular dynamics study

Ultraviolet irradiation of ice is of great interest for understanding the chemistry in both atmospheric and astrophysical environments. In interstellar space, photodissociation of H2O molecules can be a driving force behind the chemistry on icy dust grains in dense, cold molecular clouds even though the flux of UV photons is extremely low. The mechanisms of such photoinduced processes are poorly understood, however. In this work the photodissociation dynamics of a water molecule in crystalline ice at 10 K is studied computationally using classical molecular dynamics. Photodissociation in the first bilayer leads mainly to H atoms desorbing (65%), while in the third bilayer trapping of H and OH dominates (51%). The kinetic energy distribution of the desorbing H atoms is much broader than that for the corresponding gas-phase photodissociation. The H atoms on average move 11 Angstroms before becoming trapped, while OH radicals typically move 2 Angstroms. In accordance with experiments a blueshift of the absorption spectrum is obtained relative to gas-phase water.

preprint2004arXiv

Organic molecules in protoplanetary disks around TTauri and HerbigAe stars

The results of single-dish observations of low- and high-J transitions of selected molecules from protoplanetary disks around two TTauri stars (LkCa15 and TWHya) and two HerbigAe stars (HD163296 and MWC480) are reported. Simple molecules such as CO, 13CO, HCO+, CN and HCN are detected. Several lines of H2CO are found toward the TTauri star LkCa15 but not in other objects. No CH3OH has been detected down to abundances of 10E-9 - 10E-8 with respect to H2. SO and CS lines have been searched for without success. Line ratios indicate that the molecular emission arises from dense 10E6 - 10E8 cm-3 and moderately warm (T ~ 20-40K) intermediate height regions of the disk atmosphere, in accordance with predictions from models of the chemistry in disks. The abundances of most species are lower than in the envelope around the solar-mass protostar IRAS 16293-2422. Freeze-out and photodissociation are likely causes of the depletion. DCO+ is detected toward TWHya, but not in other objects. The high inferred DCO+/HCO+ ratio of ~0.035 is consistent with models of the deuterium fractionation in disks which include strong depletion of CO. The inferred ionization fraction in the intermediate height regions as deduced from HCO+ is at least 10E-11 - 10E-10, comparable to that derived for the midplane from recent H2D+ observations. (abridged abstract)

preprint2003arXiv

Interferometric Observations of Formaldehyde in the Protoplanetary Disk around LkCa15

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

preprint1998arXiv

Molecular Cloud Structure in the Magellanic Clouds: Effect of Metallicity

The chemical structure of neutral clouds in low metallicity environments is examined with particular emphasis on the H to H_2 and C+ to CO transitions. We observed near-IR H_2 lines and the CO J=1-0 line from 30 Doradus and N159/N160 in the Large Magellanic Cloud and from DEM S 16, DEM S 37, and LI-SMC 36 in the Small Magellanic Cloud. We find that the H_2 emission is UV-excited and that (weak) CO emission always exists (in our surveyed regions) toward positions where H_2 and [CII] emission have been detected. Using a PDR code and a radiative transfer code, we simulate the emission of line radiation from spherical clouds and from large planar clouds. Because the [CII] emission and H_2 emission arise on the surface of the cloud and the lines are optically thin, these lines are not affected by changes in the relative sizes of the neutral cloud and the CO bearing core, while the optically thick CO emission can be strongly affected. The sizes of clouds are estimated by measuring the deviation of CO emission strength from that predicted by a planar cloud model of a given size. The average cloud column density and therefore size increases as the metallicity decreases. Our result agrees with the photoionization regulated star formation theory by Mc Kee (1989).

preprint1996arXiv

Circumstellar C2, CN, and CH+ in the optical spectra of post-AGB stars

We present optical high-resolution spectra of a sample of sixteen post-AGB stars and IRC +10216. Of the post-AGB stars, ten show C2 Phillips and Swan and CN Red System absorption, one CH+ emission, one CH+ absorption, and four without any molecules. We find typically Trot=43-399, 155-202, and 18-50 K, log N = 14.90-15.57, 14.35, and 15.03-16.47 cm-2 for C2, CH+, and CN respectively, and 0.6<N(CN)/N(C2)<11.2. We did not detect isotopic lines, which places a lower limit on the isotope ratio of 12C/13C>20. The presence of C2 and CN absorption is correlated with cold dust (Tdust<300K) and the presence of CH+ with hot dust (Tdust>300K). All objects with the unidentified 21mum emission feature exhibit C2 and CN absorption, but not all objects with C2 and CN detections exhibit a 21mum feature. The derived expansion velocity, ranging from 5 to 44 km/s, is the same as that derived from CO millimeter line emission. This unambiguously proves that these lines are of circumstellar origin and are formed in the AGB ejecta (circumstellar shell expelled during the preceding AGB phase). Furthermore there seems to be a relation between the C2 molecular column density and the expansion velocity, which is attributed to the fact that a higher carbon abundance of the dust leads to a more efficient acceleration of the AGB wind. Using simple assumptions for the location of the molecular lines and molecular abundances, mass-loss rates have been derived from the molecular absorption lines and are comparable to those obtained from CO emission lines and the infrared excess.