Search

Context: Ultraviolet radiation (UV) influences the physics and chemistry of star-forming regions, but its properties and significance in the immediate surroundings of low-mass protostars are still poorly understood. Aims: We aim to extend the use of the CN/HCN ratio, already established for high-mass protostars, to the low-mass regime to trace and characterize the UV field around low-mass protostars on $\sim 0.6\times0.6$ pc scales. Methods: We present $5'\times5'$ maps of the Serpens Main Cloud encompassing 10 protostars observed with the EMIR receiver at the IRAM 30 m telescope in CN 1-0, HCN 1-0, CS 3-2, and some of their isotopologues. The radiative-transfer code RADEX and the chemical model Nahoon are used to determine column densities of molecules, gas temperature and density, and the UV field strength, $G_\mathrm{0}$. Results: The spatial distribution of HCN and CS are well-correlated with CO 6-5 emission that traces outflows. The CN emission is extended from the central protostars to their immediate surroundings also tracing outflows, likely as a product of HCN photodissociation. The ratio of CN to HCN total column densities ranges from $\sim$1 to 12 corresponding to G$_0$ $\approx$ $10^{1}-10^{3}$ for gas densities and temperatures typical for outflows of low-mass protostars. Conclusions: UV radiation associated with protostars and their outflows is indirectly identified in a significant part of the Serpens Main low-mass star-forming region. Its strength is consistent with the values obtained from the OH and H$_2$O ratios observed with Herschel and compared with models of UV-illuminated shocks. From a chemical viewpoint, the CN to HCN ratio is an excellent tracer of UV fields around low- and intermediate-mass star-forming regions., Comment: 32 pages, 25 figures, accepted by A\&A

In the high-mass star-forming region G35.20-0.74N, small scale (about 800 AU) chemical segregation has been observed in which complex organic molecules containing the CN group are located in a small location. We aim to determine the physical origin of the large abundance difference (about 4 orders of magnitude) in complex cyanides within G35.20-0.74 B, and we explore variations in age, gas and dust temperature, and gas density. We performed gas-grain astrochemical modeling experiments with exponentially increasing (coupled) gas and dust temperature rising from 10 to 500 K at constant H$_2$ densities of 10$^7$, 10$^8$, and 10$^9$ cm$^{-3}$. We tested the effect of varying the initial ice composition, cosmic-ray ionization rate, warm-up time (over 50, 200, and 1000 kyr), and initial (10, 15, and 25 K) and final temperatures (300 and 500 K). Varying the initial ice compositions within the observed and expected ranges does not noticeably affect the modeled abundances indicating that the chemical make-up of hot cores is determined in the warm-up stage. Complex cyanides vinyl and ethyl cyanide (CH$_2$CHCN and C$_2$H$_5$CN, respectively) cannot be produced in abundances (versus H$_2$) greater than 5x10$^{-10}$ for CH$_2$CHCN and 2x10$^{-10}$ for C$_2$H$_5$CN with a fast warm-up time (52 kyr), while the lower limit for the observed abundance of C$_2$H$_5$CN toward source B3 is 3.4x10$^{-10}$. Complex cyanide abundances are reduced at higher initial temperatures and increased at higher cosmic-ray ionization rates. Reproducing the observed abundances toward G35.20-0.74 Core B3 requires a fast warm-up at a high cosmic-ray ionization rate (1x10$^{-16}$ s$^{-1}$) at a high gas density (>10$^9$ cm$^{-3}$). G35.20-0.74 source B3 only needs to be about 2000 years older than B1/B2 for the observed chemical difference to be present. (This abstract has been shortened), 14 pages (plus 11 pages appendices), 12 figures main text

It is not well known what drives the chemistry of a protostellar envelope, in particular the role of the stellar mass and the outflows on its chemical enrichment. We study the chemical structure of NGC 1333 IRAS 4A in order to (i) investigate the influence of the outflows on the chemistry, (ii) constrain the age of our object, (iii) compare it with a typical high-mass protostellar envelope. In our analysis we use JCMT line mapping and HIFI pointed spectra. To study the influence of the outflow on the degree of deuteration, we compare JCMT maps of HCO+ and DCO+ with non-LTE (RADEX) models in a region that spatially covers the outflow activity of IRAS 4A. To study the envelope chemistry, we derive empirical molecular abundance profiles for the observed species using the radiative transfer code (RATRAN) and adopting a 1D dust density/temperature profile from the literature. We compare our best-fit observed abundance profiles with the predictions from the time dependent gas grain chemical code (ALCHEMIC). The CO, HCN, HNC and CN abundance require an enhanced UV field which points towards an outflow cavity. The abundances (wrt H2) are 1 to 2 orders of magnitude lower than those observed in the high mass protostellar envelope (AFGL 2591), while they are found to be similar within factors of a few with respect to CO. Differences in UV radiation may be responsible for such chemical differentiation, but temperature differences seem a more plausible explanation. The CH3OH modeled abundance profile points towards an age of > 4x10^4 yrs for IRAS 4A. The spatial distribution of H2D+ differs from that of other deuterated species, indicating an origin from a foreground colder layer (, 26 pages, 23 Figures, 6 tables, accepted for publication in A&A

We have used the Odin satellite to obtain strip maps of the ground-state rotational transitions of ortho-water and ortho-ammonia, as well as CO(5-4) and 13CO(5-4) across the PDR, and H218O in the central position. A physi-chemical inhomogeneous PDR model was used to compute the temperature and abundance distributions for water, ammonia and CO. A multi-zone escape probability method then calculated the level populations and intensity distributions. These results are compared to a homogeneous model computed with an enhanced version of the RADEX code. H2O, NH3 and 13CO show emission from an extended PDR with a narrow line width of ~3 kms. Like CO, the water line profile is dominated by outflow emission, however, mainly in the red wing. The PDR model suggests that the water emission mainly arises from the surfaces of optically thick, high density clumps with n(H2)>10^6 cm^-3 and a clump water abundance, with respect to H2, of 5x10^-8. The mean water abundance in the PDR is 5x10^-9, and between ~2x10^-8 -- 2x10^-7 in the outflow derived from a simple two-level approximation. Ammonia is also observed in the extended clumpy PDR, likely from the same high density and warm clumps as water. The average ammonia abundance is about the same as for water: 4x10^-9 and 8x10^-9 given by the PDR model and RADEX, respectively. The similarity of water and ammonia PDR emission is also seen in the almost identical line profiles observed close to the bright rim. Around the central position, ammonia also shows some outflow emission although weaker than water in the red wing. Predictions of the H2O(110-101) and (111-000) antenna temperatures across the PDR are estimated with our PDR model for the forthcoming observations with the Herschel Space Observatory., 13 pages, 14 figures, 10 tables. Accepted for publication in Astronomy & Astrophysics 14 November 2008

The Delta-variance analysis is an efficient tool for measuring the structural scaling behaviour of interstellar turbulence in astronomical maps. In paper I we proposed essential improvements to the Delta-variance analysis. In this paper we apply the improved Delta-variance analysis to i) a hydrodynamic turbulence simulation with prominent density and velocity structures, ii) an observed intensity map of rho Oph with irregular boundaries and variable uncertainties of the different data points, and iii) a map of the turbulent velocity structure in the Polaris Flare affected by the intensity dependence on the centroid velocity determination. The tests confirm the extended capabilities of the improved Delta-variance analysis. Prominent spatial scales were accurately identified and artifacts from a variable reliability of the data were removed. The analysis of the hydrodynamic simulations showed that the injection of a turbulent velocity structure creates the most prominent density structures are produced on a scale somewhat below the injection scale. The new analysis of a rho Oph continuum map reveals an intermediate stage in the molecular cloud evolution showing both signatures of the typical molecular cloud scaling behaviour and the formation of condensed cores. When analysing the velocity structure of the Polaris Flare we show that a universal power law connects scales from 0.03 pc to 3 pc. However, a plateau in the Delta-variance spectrum around 5 pc indicates that the visible large-scale velocity gradient is not converted directly into a turbulent cascade., Comment: Accepted for publication in A&A, Section 6

Aims. We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution.Methods. A number of benchmark models have been created, covering low and high gas densities n = 10(3), 10(5.5) cm(-3) and far ultraviolet intensities chi = 10, 10(5) in units of the Draine field (FUV: 6 Results. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes such as the effect of PAHs on the electron density or the temperature dependence of the dissociation of CO by cosmic ray induced secondary photons, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly. Nevertheless, the remaining spread in the computed observables such as the atomic fine-structure line intensities serves as a warning that there is still a considerable uncertainty when interpreting astronomical data with our models.

Aims. We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. Methods. A number of benchmark models have been created, covering low and high gas densities n = 10^3, 10^(5.5) cm^(-3) and far ultraviolet intensities χ = 10, 10^5 in units of the Draine field (FUV: 6 < h υ < 13.6 eV). The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature self consistently by solving the thermal balance, thus testing the modeling of the heating and cooling mechanisms accounting for the detailed energy balance throughout the clouds. Results. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes such as the effect of PAHs on the electron density or the temperature dependence of the dissociation of CO by cosmic ray induced secondary photons, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly. Nevertheless, the remaining spread in the computed observables such as the atomic fine-structure line intensities serves as a warning that there is still a considerable uncertainty when interpreting astronomical data with our models.

The ISO spectra of the bilobal planetary nebula Hb 5 are presented. These spectra are combined with the spectra in the visual wavelength region to obtain a complete, extinction corrected, spectrum. The chemical composition of the nebula is then calculated in several ways. First by directly calculating and adding individual ion abundances, assuming that all the ionic lines are formed in an ionized region surrounding the ionizing star. Secondly by building an "end-to-end model" nebula in which we have included a neutral region and a photodissociation region (PDR) beyond the ionized nebula. In this way we attempt to interpret the molecular hydrogen lines observed by ISO in a more self-consistent way. In the final analysis, the model is found to be basically heuristic, but gives new insights about the PDR and the PN. The implications of these are discussed.

We present near-IR speckle images of 21 massive Young Stellar Objects (YSOs) associated with outflows. The aim of this study is to search for sub-arcsecond reflection nebulae associated with the outflow cavity. We find that 6 of the massive YSOs show a conical nebula which can be interpreted in terms of reflected light from the dusty walls of the outflow cavity. In all cases, the small scale structures seen in our images are compared with outflow indicators found in the literature. No clear correlation is found between the presence of the reflection nebulosity and any property such as degree of embeddedness. We also note that 3 of the sources show close companions, one of them belonging also to the sample with conical nebula., Comment: 17 pages, 16 figures, accepted for publication in A&A

We present ground-based 2.1 to 4.2 microns observations of four low-mass protostars. We searched for the 4.1 microns OD stretch band, characteristic of solid HDO in grain mantles. We did not detect solid HDO in any of the four sources, but we derive 3-sigma upper limits from 0.5% to 2% for the HDO/H2O ratio depending on the source. These ratios provide strong constraints to solid-state deuteration models when compared to deuterium fractionation values observed in the gas phase. We discuss various scenarios that could lead to such a low water deuteration compared to the high formaldehyde and methanol deuteration observed in the gas phase., 8 pages, 6 figures Accepted for publication in A&A

We present near-IR speckle images of 21 massive Young Stellar Objects (YSOs) associated with outflows. The aim of this study is to search for sub-arcsecond reflection nebulae associated with the outflow cavity. We find that 6 of the massive YSOs show a conical nebula which can be interpreted in terms of reflected light from the dusty walls of the outflow cavity. In all cases, the small scale structures seen in our images are compared with outflow indicators found in the literature. No clear correlation is found between the presence of the reflection nebulosity and any property such as degree of embeddedness. We also note that 3 of the sources show close companions, one of them belonging also to the sample with conical nebula.

(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., Accepted for publication in A&A, 17 pages with appendix; 14 figures

The understanding of the formation process of massive stars (>8 Msun) is limited, due to theoretical complications and observational challenges. We investigate the physical structure of the large-scale (~10^4-10^5 AU) molecular envelope of the high-mass protostar AFGL2591 using spectral imaging in the 330-373 GHz regime from the JCMT Spectral Legacy Survey. Out of ~160 spectral features, this paper uses the 35 that are spatially resolved. The observed spatial distributions of a selection of six species are compared with radiative transfer models based on a static spherically symmetric structure, a dynamic spherical structure, and a static flattened structure. The maps of CO and its isotopic variations exhibit elongated geometries on scales of ~100", and smaller scale substructure is found in maps of N2H+, o-H2CO, CS, SO2, CCH, and methanol lines. A velocity gradient is apparent in maps of all molecular lines presented here, except SO, SO2, and H2CO. We find two emission peaks in warm (Eup~200K) methanol separated by 12", indicative of a secondary heating source in the envelope. The spherical models are able to explain the distribution of emission for the optically thin H13CO+ and C34S, but not for the optically thick HCN, HCO+, and CS, nor for the optically thin C17O. The introduction of velocity structure mitigates the optical depth effects, but does not fully explain the observations, especially in the spectral dimension. A static flattened envelope viewed at a small inclination angle does slightly better. We conclude that a geometry of the envelope other than an isotropic static sphere is needed to circumvent line optical depth effects. We propose that this could be achieved in envelope models with an outflow cavity and/or inhomogeneous structure at scales smaller than ~10^4 AU. The picture of inhomogeneity is supported by observed substructure in at least six species., 17 pages; accepted for publication in A&A