Source author record

Ruud Visser

Ruud Visser 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

27works
6topics
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

27 published item(s)

preprint2020arXiv

Missing water in Class I protostellar disks

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

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

Hydrocarbon emission rings in protoplanetary disks induced by dust evolution

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

preprint2015arXiv

Chemical tracers of episodic accretion in low-mass protostars

Aims: Accretion rates in low-mass protostars can be highly variable in time. Each accretion burst is accompanied by a temporary increase in luminosity, heating up the circumstellar envelope and altering the chemical composition of the gas and dust. This paper aims to study such chemical effects and discusses the feasibility of using molecular spectroscopy as a tracer of episodic accretion rates and timescales. Methods: We simulate a strong accretion burst in a diverse sample of 25 spherical envelope models by increasing the luminosity to 100 times the observed value. Using a comprehensive gas-grain network, we follow the chemical evolution during the burst and for up to 10^5 yr after the system returns to quiescence. The resulting abundance profiles are fed into a line radiative transfer code to simulate rotational spectra of C18O, HCO+, H13CO+, and N2H+ at a series of time steps. We compare these spectra to observations taken from the literature and to previously unpublished data of HCO+ and N2H+ 6-5 from the Herschel Space Observatory. Results: The bursts are strong enough to evaporate CO throughout the envelope, which in turn enhances the abundance of HCO+ and reduces that of N2H+. After the burst, it takes 10^3-10^4 yr for CO to refreeze and for HCO+ and N2H+ to return to normal. The chemical effects of the burst remain visible in the rotational spectra for as long as 10^5 yr after the burst has ended, highlighting the importance of considering luminosity variations when analyzing molecular line observations in protostars. The spherical models are currently not accurate enough to derive robust timescales from single-dish observations. As follow-up work, we suggest that the models be calibrated against spatially resolved observations in order to identify the best tracers to be used for statistically significant source samples.

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

Molecule sublimation as a tracer of protostellar accretion: Evidence for accretion bursts from high angular resolution C18O images

The accretion histories of embedded protostars are an integral part of descriptions of their physical and chemical evolution. In particular, are the accretion rates smoothly declining from the earlier toward later stages or in fact characterized by variations such as intermittent bursts? We aim to characterize the impact of possible accretion variations in a sample of embedded protostars by measuring the size of the inner regions of their envelopes where CO is sublimated and relate those to their temperature profiles dictated by their current luminosities. Using observations from the Submillimeter Array we measure the extents of the emission from the C18O isotopologue toward 16 deeply embedded protostars. We compare these measurements to the predicted extent of the emission given the current luminosities of the sources through dust and line radiative transfer calculations. Eight out of sixteen sources show more extended C18O emission than predicted by the models. The modeling shows that the likely culprit for these signatures is sublimation due to increases in luminosities of the sources by about a factor five or more during the recent 10,000 years - the time it takes for CO to freeze-out again on dust grains. For four of those sources the increase would have had to have been a factor 10 or more. The compact emission seen toward the other half of the sample suggests that C18O only sublimates when the temperature exceeds 30 K - as one would expect if CO is mixed with H2O in the grain ice-mantles. The small-number statistics from this survey suggest that protostars undergo significant bursts about once every 20,000 years. This also illustrates the importance of taking the physical evolutionary histories into account for descriptions of the chemical structures of embedded protostars.

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.

preprint2014arXiv

ALMA observations of the kinematics and chemistry of disc formation

Context: The R CrA cloud hosts a handful of Class 0/I low-mass young stellar objects. The chemistry and physics at scales $>500$ AU in this cloud are dominated by the irradiation from the nearby Herbig Be star R CrA. The luminous large-scale emission makes it necessary to use high-resolution spectral imaging to study the chemistry and dynamics of the inner envelopes and discs of the protostars. Aims: We aim to better understand the structure of the inner regions of these protostars and, in particular, the interplay between the chemistry and the presence of discs. Methods: Using Atacama Large Millimeter/submillimeter Array (ALMA) high-resolution spectral imaging interferometry observations, we study the molecular line and dust continuum emission at submillimetre wavelengths. Results: We detect dust continuum emission from four circumstellar discs around Class 0/I objects within the R CrA cloud. Towards IRS7B we detect C$^{17}$O emission showing a rotation curve consistent with a Keplerian disc with a well-defined edge that gives a good estimate for the disc radius at 50 AU. We derive the central object mass to $2.3M_{\odot}$ and the disc mass to $0.024M_{\odot}$. The observations are also consistent with a model of material infalling under conservation of angular momentum; however, this model provides a worse fit to the data. We also report a likely detection of faint CH$_3$OH emission towards this point source, as well as more luminous CH$_3$OH emission in an outflow orthogonal to the major axis of the C$^{17}$O emission. Conclusions: The faint CH$_3$OH emission seen towards IRS7B can be explained by a flat density profile of the inner envelope caused by the disc with a radius $\lesssim50$ AU. We propose that the regions of the envelopes where complex organic molecules are present in Class 0/I young stellar objects can become quenched as the disc grows.

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

Star Formation: Chemistry as a Probe of Embedded Protostars

The embedded phase of star formation is the crucial phase where most of the stellar mass is assembled. Velocity-resolved spectra reveal an infalling envelope, bipolar outflows, and perhaps an infant circumstellar disk -- all locked together in a cosmic dance of gravitational collapse and magnetic winds. Densities and temperatures change by orders of magnitude as the protostar evolves, driving a chemistry as exotic as it is fascinating. I will review two examples of how to exploit chemistry and molecular spectroscopy to study the physics of low-mass star formation: energetic feedback and episodic accretion.

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

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

Fundamental aspects of episodic accretion chemistry explored with single-point models

We explore a set of single-point chemical models to study the fundamental chemical aspects of episodic accretion in low-mass embedded protostars. Our goal is twofold: (1) to understand how the repeated heating and cooling of the envelope affects the abundances of CO and related species; and (2) to identify chemical tracers that can be used as a novel probe of the timescales and other physical aspects of episodic accretion. We develop a set of single-point models that serve as a general prescription for how the chemical composition of a protostellar envelope is altered by episodic accretion. The main effect of each accretion burst is to drive CO ice off the grains in part of the envelope. The duration of the subsequent quiescent stage (before the next burst hits) is similar to or shorter than the freeze-out timescale of CO, allowing the chemical effects of a burst to linger long after the burst has ended. We predict that the resulting excess of gas-phase CO can be observed with single-dish or interferometer facilities as evidence of an accretion burst in the past 10^3 - 10^4 yr.

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

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

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

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

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

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

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

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

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

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.

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.