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. T
his 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.
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.
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.
Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is
to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.
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 l
ayers 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.
Aims. Photodissociation by UV light is an important destruction mechanism for CO in many astrophysical environments, ranging from interstellar clouds to protoplanetary disks. The aim of this work is to gain a better understanding of the depth depende
nce and isotope-selective nature of this process. Methods. We present a photodissociation model based on recent spectroscopic data from the literature, which allows us to compute depth-dependent and isotope-selective photodissociation rates at higher accuracy than in previous work. The model includes self-shielding, mutual shielding and shielding by atomic and molecular hydrogen, and it is the first such model to include the rare isotopologues C17O and 13C17O. We couple it to a simple chemical network to analyse CO abundances in diffuse and translucent clouds, photon-dominated regions, and circumstellar disks. Results. The photodissociation rate in the unattenuated interstellar radiation field is 2.6e-10 s^-1, 30% higher than currently adopted values. Increasing the excitation temperature or the Doppler width can reduce the photodissociation rates and the isotopic selectivity by as much as a factor of three for temperatures above 100 K. The model reproduces column densities observed towards diffuse clouds and PDRs, and it offers an explanation for both the enhanced and the reduced N(12CO)/N(13CO) ratios seen in diffuse clouds. The photodissociation of C17O and 13C17O shows almost exactly the same depth dependence as that of C18O and 13C18O, respectively, so 17O and 18O are equally fractionated with respect to 16O. This supports the recent hypothesis that CO photodissociation in the solar nebula is responsible for the anomalous 17O and 18O abundances in meteorites.
