Search

Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration throughout the prototypical pre-stellar core L1544. Methods: We use single-dish observations of c-C$_3$H$_2$, c-H$^{13}$CC$_2$H, c-C$_3$HD, and c-C$_3$D$_2$ towards the pre-stellar core L1544, performed at the IRAM 30m telescope. We derive the column density and deuterium fraction maps, and compare these observations with non-LTE radiative transfer simulations. Results: The highest deuterium fractions are found close to the dust peak at the centre of L1544, where the increased abundance of H$_2$D$^+$ ions drives the deuteration process. The peak values are N(c-C$_3$HD)/N(c-C$_3$H$_2)=0.17\pm0.01$, N(c-C$_3$D$_2$)/N(c-C$_3$H$_2)=0.025\pm0.003$ and N(c-C$_3$D$_2$)/N(c-C$_3$HD$)=0.16\pm0.03$, which is consistent with previous single point observations. The distributions of c-C$_3$HD and c-C$_3$D$_2$ indicate that the deuterated forms of c-C$_3$H$_2$ in fact trace the dust peak and not the c-C$_3$H$_2$ peak. Conclusions: The N(c-C$_3$D$_2$)/N(c-C$_3$HD) map confirms that the process of deuteration is more efficient towards the centre of the core and demonstrates that carbon-chain molecules are still present at high densities. This is likely caused by an increased abundance of He$^+$ ions destroying CO, which increases the amount of carbon atoms in the gas phase., Comment: 13 pages, 9 figures

We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center., Comment: 20 pages, 8 figures, 3 tables. ApJ accepted

Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules., Comment: Accepted for publication on The Astrophysical Journal

Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars., Comment: 7 pages, 1 table, 5 figures, A&A Letters, in press

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations., Comment: 12 pages, 8 figures, accepted for publication in A&A

We report extensive theoretical calculations on the rotation-inversion excitation of interstellar ammonia (NH3) due to collisions with atomic and molecular hydrogen (both para- and ortho-H2). Close-coupling calculations are performed for total energies in the range 1-2000 cm-1 and rotational cross sections are obtained for all transitions among the lowest 17 and 34 rotation-inversion levels of ortho- and para-NH3, respectively. Rate coefficients are deduced for kinetic temperatures up to 200 K. Propensity rules for the three colliding partners are discussed and we also compare the new results to previous calculations for the spherically symmetrical He and para-H2 projectiles. Significant differences are found between the different sets of calculations. Finally, we test the impact of the new rate coefficients on the calibration of the ammonia thermometer. We find that the calibration curve is only weakly sensitive to the colliding partner and we confirm that the ammonia thermometer is robust.

Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K., Comment: Accepted for publication in A&A

New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks., Comment: 8 pages, 2 tables, 4 figures, accepted for publication in MNRAS

The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios.

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

The availability of collisional rate coefficients with H$_2$ is a pre-requisite for interpretation of observations of molecules whose energy levels are populated under non local thermodynamical equilibrium conditions. In the current study, we present collisional rate coefficients for the NH$_2$D / para--H$_2$($J_2 = 0,2$) collisional system, for energy levels up to $J_\tau = 7_7$ ($E_u$$\sim$735 K) and for gas temperatures in the range $T = 5-300$K. The cross sections are obtained using the essentially exact close--coupling (CC) formalism at low energy and at the highest energies, we used the coupled--states (CS) approximation. For the energy levels up to $J_\tau = 4_2$ ($E_u$$\sim$215 K), the cross sections obtained through the CS formalism are scaled according to a few CC reference points. These reference points are subsequently used to estimate the accuracy of the rate coefficients for higher levels, which is mainly limited by the use of the CS formalism. Considering the current potential energy surface, the rate coefficients are thus expected to be accurate to within 5\% for the levels below $J_\tau = 4_2$, while we estimate an accuracy of 30\% for higher levels.

Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest shock region B1 of the L1157 molecular outflow. We observed several ortho- and para-H2O transitions using HIFI and PACS instruments on board Herschel, providing a detailed picture of the kinematics and spatial distribution of the gas. We performed a LVG analysis to derive the physical conditions of H2O shocked material, and ultimately obtain its abundance. We detected 13 H2O lines probing a wide range of excitation conditions. PACS maps reveal that H2O traces weak and extended emission associated with the outflow identified also with HIFI in the o-H2O line at 556.9 GHz, and a compact (~10") bright, higher-excitation region. The LVG analysis of H2O lines in the bow-shock show the presence of two gas components with different excitation conditions: a warm (Tkin~200-300 K) and dense (n(H2)~(1-3)x10^6 cm-3) component with an assumed extent of 10" and a compact (~2"-5") and hot, tenuous (Tkin~900-1400 K, n(H2)~10^3-10^4 cm-3) gas component, which is needed to account for the line fluxes of high Eu transitions. The fractional abundance of the warm and hot H2O gas components is estimated to be (0.7-2)x10^{-6} and (1-3)x10^{-4}, respectively. Finally, we identified an additional component in absorption in the HIFI spectra of H2O lines connecting with the ground state level, probably arising from the photodesorption of icy mantles of a water-enriched layer at the edges of the cloud., Comment: Accepted for publication in A&A. 12 pages, 9 figures, 4 tables

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources., Comment: Accepted for publication in ApJS

We present the first detection of hydrogen chlorine in a protostellar shock, by observing the fundamental transition at 626 GHz with the Herschel HIFI spectrometer. We detected two of the three hyperfine lines, from which we derived a line opacity < 1. Using a non-LTE LVG code, we constrained the HCl column density, temperature and density of the emitting gas. The hypothesis that the emission originates in the molecular cloud is ruled out, as it would imply a too dense gas. Conversely, assuming that the emission originates in the 10"-15" size shocked gas previously observed at the IRAM PdB interferometer, we obtain: N(HCl)=0.7-2 x 10(13) cm-2, temperature > 15 K and density > 3 x 10(5) cm-3}. Combining with the Herschel HIFI CO(5-4) observations allows to further constrain the gas density and temperature, 10(5)-10(6) cm-3 and 120-250 K, as well as the HCl column density, 2 x 10(13) cm-2, and, finally, abundance: 3-6 x 10(-9). The estimated HCl abundance is consistent with that previously observed in low- and high- mass protostars. This puzzling result in the L1157-B1 shock, where species from volatile and refractory grains components are enhanced, suggests either that HCl is not the main reservoir of chlorine in the gas phase, against previous chemical models predictions, or that the elemental chlorine abundance is low in L1157-B1. Astrochemical modelling suggests that HCl is in fact formed in the gas phase, at low temperatures, prior to the occurance of the shock, and that the latter does not enhance its abundance., Comment: ApJ, in press

State-to-state differential cross sections (DCSs) for rotationally inelastic scattering of H2O by H2 have been measured at 71.2 meV (574 cm-1) and 44.8 meV (361 cm-1) collision energy using crossed molecular beams combined with velocity map imaging. A molecular beam containing variable compositions of the (J = 0, 1, 2) rotational states of hydrogen collides with a molecular beam of argon seeded with water vapor that is cooled by supersonic expansion to its lowest para or ortho rotational levels (JKaKc= 000 and 101, respectively). Angular speed distributions of fully specified rotationally excited final states are obtained using velocity map imaging. Relative integral cross sections are obtained by integrating the DCSs taken with the same experimental conditions. Experimental state-specific DCSs are compared with predictions from fully quantum scattering calculations on the most complete H2O-H2 potential energy surface. Comparison of relative total cross sections and state-specific DCSs show excellent agreement with theory in almost all details, Comment: 46 pages

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1., Comment: 5 pages, 5 figures, accepted the A&A HIFI Special Issue as a letter

The interpretation of water line emission from existing observations and future HIFI/Herschel data requires a detailed knowledge of collisional rate coefficients. Among all relevant collisional mechanisms, the rotational (de)excitation of H2O by H2 molecules is the process of most interest in interstellar space. To determine rate coefficients for rotational de-excitation among the lowest 45 para and 45 ortho rotational levels of H2O colliding with both para and ortho-H2 in the temperature range 20-2000 K. Rate coefficients are calculated on a recent high-accuracy H2O-H2 potential energy surface using quasi-classical trajectory calculations. Trajectories are sampled by a canonical Monte-Carlo procedure. H2 molecules are assumed to be rotationally thermalized at the kinetic temperature. By comparison with quantum calculations available for low lying levels, classical rates are found to be accurate within a factor of 1-3 for the dominant transitions, that is those with rates larger than a few 10^{-12}cm^{3}s^{-1}. Large velocity gradient modelling shows that the new rates have a significant impact on emission line fluxes and that they should be adopted in any detailed population model of water in warm and hot environments., Comment: 8 pages, 2 figures, 1 table (the online material (4 tables) can be obtained upon request to afaure@obs.ujf-grenoble.fr)

Cross sections for the rotational (de)excitation of CO by ground state para- and ortho-H_2 are obtained using quantum scattering calculations for collision energies between 1 and 520 cm^{-1}. A new CO-H_2 potential energy surface is employed and its quality is assessed by comparison with explicitly correlated CCSD(T)-R12 calculations. Rate constants for rotational levels of CO up to 5 and temperatures in the range 5-70 K are deduced. The new potential is found to have a strong influence on the resonance structure of the cross sections at very low collision energies. As a result, the present rates at 10 K differ by up to 50% with those obtained by \citet{flower01} on a previous, less accurate, potential energy surface., Comment: 8 pages, including 2 "online material" pages, 2 figures