Source author record

E. Roueff

E. Roueff 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

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

44 published item(s)

preprint2026arXiv

New Estimate for the Cosmic Ray-Induced $\rm H_2$ Photodissociation Rate in the Interstellar Medium

In the interstellar medium, cosmic rays (CRs) generate a field of ultraviolet (UV) photons via the excitation and subsequent radiative decay of $\rm H_2$ molecules. This UV field is a major agent of ionization and dissociation in the inner regions of molecular clouds that are shielded from the effects of the interstellar radiation field. In particular, the dissociation of $\rm H_2$, by far the most abundant molecule in interstellar clouds, leads to the production of atomic hydrogen which then takes part in the production of a multitude of molecules, in particular complex organics on the surfaces of interstellar dust grains. Precise knowledge of the rates of CR-induced dissociation processes is thus crucial for constructing reliable chemical models. For the present paper, we have derived a new value of $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ for the rate of $\rm H_2$ dissociation, where $ζ$ is the CR ionization rate of $\rm H_2$. This prediction contrasts a previous value from the Leiden database which overestimated the rate due to an inconsistent treatment of the $\rm H_2$ abundances and photodissociation cross sections. By running a series of chemical models, we show that the overestimated dissociation rate has a large effect on the results of chemical simulations, with the abundance of methanol being overestimated by over one order of magnitude. Hence, we strongly recommend the adoption of our new estimate $k_{\rm diss, CR}(\mbox{$\rm H_2$})=0.831ζ$ in all chemical models that include this process. Our newly derived value corresponds to $\rm H_2$ being purely in the para form ($J^{\prime\prime} = 0$). However, in the interiors of molecular clouds the $\rm H_2$ ortho-to-para ratio is low and using the rate for para-$\rm H_2$ is an adequate approximation.

preprint2022arXiv

Dissociative recombination and rotational transitions of D$_2^+$ in collisions with slow electrons

Rate coefficients for dissociative recombination and state-to-state rotational transitions of the D$_{2}^{+}$ ion induced by collisions with very low-energy electrons have been reported following our previous studies on HD$^{+}$ and H$_{2}^{+}$ [9,10]. The same molecular structure data sets, excitations ($N_{i}^{+} \rightarrow$ $N_{f}^{+}=N_{i}^{+}+2$ for $N_{i}^{+}=0$ to $10$) and de-excitations ($N_{i}^{+}$ $\rightarrow$ $N_{f}^{+}=N_{i}^{+}-2$, for $N_{i}^{+}=2$ to $10$) were used for collision energies ranging from $0.01$ meV to $0.3$ eV. Isotopic effects for dissociative recombination and rotational transitions of the vibrationally relaxed targets are presented.

preprint2022arXiv

Gas phase Elemental abundances in Molecular cloudS (GEMS) VI. A sulphur journey across star-forming regions: study of thioformaldehyde emission

In the context of the IRAM 30m Large Program GEMS, we present a study of thioformaldehyde in several starless cores located in star-forming filaments of Taurus, Perseus, and Orion. We investigate the influence of the environmental conditions on the abundances of these molecules in the cores, and the effect of time evolution. We have modelled the observed lines of H2CS, HDCS, and D2CS using the radiative transfer code RADEX. We have also used the chemical code Nautilus to model the evolution of these species depending on the characteristics of the starless cores. We derive column densities and abundances for all the cores. We also derive deuterium fractionation ratios, Dfrac, to determine and compare the evolutionary stage between different parts of each star-forming region. Our results indicate that the north region of the B213 filament in Taurus is more evolved than the south, while the north-eastern part of Perseus presents an earlier evolutionary stage than the south-western zone. Model results also show that Dfrac decreases with the cosmic-ray ionisation rate, while it increases with density and with the degree of sulphur depletion. In particular, we only reproduce the observations when the initial sulphur abundance in the starless cores is at least one order of magnitude lower than the solar elemental sulphur abundance. The progressive increase in HDCS/H2CS and D2CS/H2CS with time makes these ratios powerful tools for deriving the chemical evolutionary stage of starless cores. However, they cannot be used to derive the temperature of these regions, since both ratios present a similar evolution at two different temperature ranges (7-11 K and 15-19 K). Regarding chemistry, (deuterated) thioformaldehyde is mainly formed through gas-phase reactions (double-replacement and neutral-neutral displacement reactions), while surface chemistry plays an important role as a destruction mechanism.

preprint2021arXiv

New deuterated species in TMC-1: Detection of CH2DC4H with the QUIJOTE line survey

We report the first detection in space of the single deuterated isotopologue of methyldiacetylene, CH2DC4H. A total of 12 rotational transitions, with J = 8-12 and Ka = 0 and 1, were identified for this species in TMC-1 in the 31.0-50.4 GHz range using the Yebes 40m radio telescope. The observed frequencies allowed us to obtain, for the first time, the spectroscopic parameters of this deuterated isotopologue. We derived a column density of (5.5+/-0.2)e11 cm-2. The abundance ratio between CH3C4H and CH2DC4H is 24+/-2. This ratio is similar to that found for the CH3C3N/CH2DC3N analogue system, which is 22+/-2.We did not detect the deuterated species CH3C4D, which has already been observed in laboratory experiments. The detection of deuterated CH3C4H allows us to extend the discussion on the chemical mechanisms of deuterium fractionation at work in TMC-1 using a new gas-phase chemical model with multiply deuterated molecules. Introducing a possible deuterium exchange reaction between CH3CCH and atomic deuterium allows us to account for the CH3C4H/CH2DC4H abundance ratio.

preprint2020arXiv

Detection of deuterated molecules, but not of lithium hydride, in the z=0.89 absorber toward PKS1830-211

Deuterium and lithium are light elements of high cosmological and astrophysical importance. In this work we report the first detection of deuterated molecules and a search for lithium hydride, 7LiH, at redshift z=0.89 in the spiral galaxy intercepting the line of sight to the quasar PKS1830-211. We used ALMA to observe several submillimeter lines of ND, NH2D, and HDO, and their related isotopomers NH2, NH3, and H2^18O, in absorption against the southwest image of the quasar, allowing us to derive XD/XH abundance ratios. The absorption spectra mainly consist of two distinct narrow velocity components for which we find remarkable differences. One velocity component shows XD/XH abundances that is about 10 times larger than the primordial elemental D/H ratio, and no variability of the absorption profile during the time span of our observations. [...] The second component has XD/XH abundances that are 100 times larger than the primordial D/H ratio, a deepening of the absorption by a factor of two within a few months, and a rich chemical composition, with relative enhancements of N2H+, CH3OH, SO2, and complex organic molecules. We therefore speculate that this component is associated with the analog of a Galactic dark cloud, while the first component is likely more diffuse. Our search for the 7LiH (1--0) line was unsuccessful and we derive an upper limit 7LiH/H2 = 4 x 10^-13 (3sigma) in the z=0.89 absorber toward PKS1830-211. Besides, with ALMA archival data, we could not confirm the previous tentative detections of this line in the z=0.68 absorber toward B0218+357; we derive an upper limit 7LiH/H2 = 5 x 10^-11 (3sigma), although this is less constraining than our limit toward PKS1830-211. We conclude that, as in the Milky Way, only a tiny fraction of lithium nuclei are possibly bound in LiH in these absorbers at intermediate redshift.

preprint2020arXiv

Gas phase Elemental abundances in Molecular cloudS (GEMS). II. On the quest for the sulphur reservoir in molecular clouds: the $H_{2}S$ case

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

preprint2020arXiv

Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction

CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance.

preprint2016arXiv

A milestone toward understanding PDR properties in the extreme environment of LMC-30Dor

More complete knowledge of galaxy evolution requires understanding the process of star formation and interaction between the interstellar radiation field and the interstellar medium in galactic environments traversing a wide range of physical parameter space. Here we focus on the impact of massive star formation on the surrounding low metallicity ISM in 30 Doradus in the Large Magellanic Cloud. A low metal abundance, as is the case of some galaxies of the early universe, results in less ultra-violet shielding for the formation of the molecular gas necessary for star formation to proceed. The half-solar metallicity gas in this region is strongly irradiated by the super star cluster R136, making it an ideal laboratory to study the structure of the ISM in an extreme environment. Our spatially resolved study investigates the gas heating and cooling mechanisms, particularly in the photo-dissociation regions where the chemistry and thermal balance are regulated by far-ultraviolet photons (6 eV< hν<13.6 eV). We present Herschel observations of far-infrared fine-structure lines obtained with PACS and SPIRE/FTS. We have combined atomic fine-structure lines from Herschel and Spitzer observations with ground-based CO data to provide diagnostics on the properties and the structure of the gas by modeling it with the Meudon PDR code. We derive the spatial distribution of the radiation field, the pressure, the size, and the filling factor of the photodissociated gas and molecular clouds. We find a range of pressure of ~ 10^5 - 1.7x10^6 cm^{-3} K and a range of incident radiation field G_UV ~ 10^2 - 2.5x10^4 through PDR modeling. Assuming a plane-parallel geometry and a uniform medium, we find a total extinction of 1-3 mag , which correspond to a PDR cloud size of 0.2 to 3pc, with small CO depth scale of 0.06 to 0.5pc. We also determine the three dimensional structure of the gas. (Abridged)

preprint2016arXiv

A New Reference Chemical Composition for TMC-1

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi & Kaifu and the values reported in the review from Agundez & Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

preprint2016arXiv

Collisional excitation of doubly and triply deuterated ammonia ND$_2$H and ND$_3$ by H$_2$

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.

preprint2016arXiv

N$_2$H$^+$ and N$^{15}$NH$^+$ towards the prestellar core 16293E in L1689N

Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statistically significant studies. In particular, the N$_2$H$^+$ molecule currently shows the largest spread of $^{14}$N/$^{15}$N ratios in high-mass star forming regions. However, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ was obtained in only two low-mass star forming regions (L1544 and B1b). The current work extends this sample to a third dark cloud. We targeted the 16293E prestellar core, where the N$^{15}$NH$^+$ $J$=1-0 line was detected. Using a model previously developed for the physical structure of the source, we solved the molecular excitation with a non-local radiative transfer code. For that purpose, we computed specific collisional rate coefficients for the N$^{15}$NH$^+$-H$_2$ collisional system. As a first step of the analysis, the N$_2$H$^+$ abundance profile was constrained by reproducing the N$_2$H$^+$ $J$=1-0 and 3-2 maps. A scaling factor was then applied to this profile to match the N$^{15}$NH$^+$ $J$=1-0 spectrum. We derive a column density ratio N$_2$H$^+$ / N$^{15}$NH$^+$ = $330^{+170}_{-100}$. The current estimate $\sim$330 agrees with the value typical of the elemental isotopic ratio in the local ISM. It is however lower than in some other cores, where values as high as 1300 have been reported.

preprint2016arXiv

Star Formation and Feedback: A Molecular Outflow-Prestellar Core Interaction in L1689N

We present Herschel, ALMA Compact Array (ACA), and Caltech Submillimeter Observatory (CSO) observations of the prestellar core in L1689N, which has been suggested to be interacting with a molecular outflow driven by the nearby solar type protostar IRAS 16293-2422. This source is characterized by some of the highest deuteration levels seen in the interstellar medium. The change in the NH2D line velocity and width across the core provides clear evidence of an interaction with the outflow, traced by the high-velocity water emission. Quiescent, cold gas, characterized by narrow line widths is seen in the NE part of the core, while broader, more disturbed line profiles are seen in the W/SW part. Strong N2D+ and ND3 emission is detected with the ACA, extending S/SW from the peak of the single-dish NH2D emission. The ACA data also reveal the presence a compact dust continuum source, with a mean size of ~1100 au, a central density of (1-2) 10^7 cm-3, and a mass of 0.2-0.4 Msun. The dust emission peak is displaced ~5" to the south with respect to the N2D+ and ND3 emission, as well as the single-dish dust continuum peak, suggesting that the northern, quiescent part of the core is characterized by spatially extended continuum emission, which is resolved out by the interferometer. We see no clear evidence of fragmentation in this quiescent part of the core, which could lead to a second generation of star formation, although a weak dust continuum source is detected in this region in the ACA data.

preprint2016arXiv

The first CO+ image: Probing the HI/H2 layer around the ultracompact HII region Mon R2

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

preprint2016arXiv

The NH$_2$D hyperfine structure revealed by astrophysical observations

The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.

preprint2015arXiv

Herschel HIFI Observations of the Sgr A +50 km/s Cloud. Deep Searches for O2 in Emission and Foreground Absorption

To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).

preprint2015arXiv

Isotopic fractionation of carbon, deuterium and nitrogen : a full chemical study

Context. The increased sensitivity and high spectral resolution of millimeter telescopes allow the detection of an increasing number of isotopically substituted molecules in the interstellar medium. The 14N/ 15N ratio is difficult to measure directly for carbon containing molecules. Aims. We want to check the underlying hypothesis that the 13C/ 12C ratio of nitriles and isonitriles is equal to the elemental value via a chemical time dependent gas phase chemical model. Methods. We have built a chemical network containing D, 13C and 15N molecular species after a careful check of the possible fractionation reactions at work in the gas phase. Results. Model results obtained for 2 different physical conditions corresponding respectively to a moderately dense cloud in an early evolutionary stage and a dense depleted pre-stellar core tend to show that ammonia and its singly deuterated form are somewhat enriched in 15N, in agreement with observations. The 14N/ 15N ratio in N2H+ is found to be close to the elemental value, in contrast to previous models which obtain a significant enrichment, as we found that the fractionation reaction between 15N and N2H+ has a barrier in the entrance channel. The large values of the N2H+/15NNH+ and N2H+/ N15NH+ ratios derived in L1544 cannot be reproduced in our model. Finally we find that nitriles and isonitriles are in fact significantly depleted in 13C, questioning previous interpretations of observed C15N, HC15N and H15NC abundances from 13C containing isotopologues.

preprint2015arXiv

Probing non polar interstellar molecules through their protonated form: Detection of protonated cyanogen (NCCNH+)

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

preprint2015arXiv

Spatially resolved l-c3h+ emission in the horsehead photodissociation region: Further evidence for a top-down hydrocarbon chemistry

Small hydrocarbons, such as C2H, C3H and C3H2 are more abundant in photo-dissociation regions (PDRs) than expected based on gas-phase chemical models. To explore the hydrocarbon chemistry further, we observed a key intermediate species, the hydrocarbon ion l-C3H+, in the Horsehead PDR with the Plateau de Bure Interferometer at high-angular resolution (6''). We compare with previous observations of C2H and c-C3H2 at similar angular resolution and new gas-phase chemical model predictions to constrain the dominant formation mechanisms of small hydrocarbons in low-UV flux PDRs. We find that, at the peak of the HCO emission (PDR position), the measured l-C3H+, C2H and c-C3H2 abundances are consistent with current gas-phase model predictions. However, in the first PDR layers, at the 7.7 mum PAH band emission peak, which are more exposed to the radiation field and where the density is lower, the C2H and c-C3H2 abundances are underestimated by an order of magnitude. At this position, the l-C3H+ abundance is also underpredicted by the model but only by a factor of a few. In addition, contrary to the model predictions, l-C3H+ peaks further out in the PDR than the other hydrocarbons, C2H and c-C3H2. This cannot be explained by an excitation effect. Current gas-phase photochemical models thus cannot explain the observed abundances of hydrocarbons, in particular in the first PDR layers. Our observations are consistent with a top-down hydrocarbon chemistry, in which large polyatomic molecules or small carbonaceous grains are photo-destroyed into smaller hydrocarbon molecules/precursors.

preprint2014arXiv

Collisional excitation of singly deuterated ammonia NH$_2$D by H$_2$

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_τ= 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_τ= 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_τ= 4_2$, while we estimate an accuracy of 30\% for higher levels.

preprint2014arXiv

Deuteration around the ultracompact HII region Mon R2

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.

preprint2014arXiv

H$_2$ Lyman and Werner band lines and their sensitivity for a variation of the proton-electron mass ratio in the gravitational potential of white dwarfs

Recently an accurate analysis of absorption spectra of molecular hydrogen, observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope, in the photosphere of white dwarf stars GD133 and GD29-38 was published in a Letter [Phys. Rev. Lett. 113, 123002 (2014)], yielding a constraint on a possible dependence of the proton-electron mass ratio on a gravitational field of strength 10,000 times that at the Earth's surface. In the present paper further details of that study are presented, in particular a re-evaluation of the spectrum of the $B^1Σ_u^+ - X^1Σ_g^+ (v',v'')$ Lyman bands relevant for the prevailing temperatures (12,000 - 14,000 K) of the photospheres. An emphasis is on the calculation of so-called $K_i$-coefficients, that represent the sensitivity of each individual line to a possible change in the proton-electron mass ratio. Such calculations were performed by semi-empirical methods and by ab initio methods providing accurate and consistent values. A full listing is provided for the molecular physics data on the Lyman bands (wavelengths $λ_i$, line oscillator strengths $f_i$, radiative damping rates $Γ_i$, and sensitivity coefficients $K_i$) as required for the analyses of H$_2$-spectra in hot dwarf stars. A similar listing of the molecular physics parameters for the $C^1Π_u - X^1Σ_g^+ (v',v'')$ Werner bands is provided for future use in the analysis of white dwarf spectra.

preprint2014arXiv

Photodissociation of interstellar ArH+

Aims. Following the recent detection of 36ArH+ in the Crab nebula spectrum, we have computed the photodissociation rate of ArH+ in order to constrain the physical processes at work in this environment. Methods. Photodissociation cross sections of ArH+ are computed in an ab initio approach including explicit account of spin-orbit coupling. Results. We report the photodissociation cross section of ArH+ as a function of wavelength. Photodissociation probabilities are derived for different impinging radiation fields.The photodissociation probability of for a very small unshielded cloud surrounded on all sides by the unshielded InterStellar Radiation Field (ISRF) model described by Draine (1978) is equal to 9.9e-12 s-1 and 1.9e-9 s-1 in the Crab nebula conditions. The dependence on the visual extinction is obtained by using the Meudon Photon Dominated Region (PDR) code and corresponding analytical fits are provided. Conclusions. These data will help to produce a realistic chemical network to interpret the observations. Photodissociation of ArH+ is found to be moderate and the presence of this molecular ion is mainly dependent on the molecular fraction

preprint2014arXiv

Tentative Detection of the Nitrosylium Ion in Space

We report the tentative detection in space of the nitrosylium ion, NO$^+$. The observations were performed towards the cold dense core Barnard 1-b. The identification of the NO$^+$ $J$=2--1 line is supported by new laboratory measurements of NO$^+$ rotational lines up to the $J$=8--7 transition (953207.189\,MHz), which leads to an improved set of molecular constants: $B_0 = 59597.1379(62)$\,MHz, $D_0 = 169.428(65)$\,kHz, and $eQq_0(\textrm{N}) = -6.72(15)$\,MHz. The profile of the feature assigned to NO$^+$ exhibits two velocity components at 6.5 and 7.5 km s$^{-1}$, with column densities of $1.5 \times 10^{12}$ and $6.5\times10^{11}$ cm$^{-2}$, respectively. New observations of NO and HNO, also reported here, allow to estimate the following abundance ratios: $X$(NO)/$X$(NO$^+$)$\simeq511$, and $X$(HNO)/$X$(NO$^+$)$\simeq1$. This latter value provides important constraints on the formation and destruction processes of HNO. The chemistry of NO$^+$ and other related nitrogen-bearing species is investigated by the means of a time-dependent gas phase model which includes an updated chemical network according to recent experimental studies. The predicted abundance for NO$^+$ and NO is found to be consistent with the observations. However, that of HNO relative to NO is too high. No satisfactory chemical paths have been found to explain the observed low abundance of HNO. HSCN and HNCS are also reported here with an abundance ratio of $\simeq1$. Finally, we have searched for NNO, NO$_2$, HNNO$^+$, and NNOH$^+$, but only upper limits have been obtained for their column density, except for the latter for which we report a tentative 3-$σ$ detection.

preprint2013arXiv

CH2D+, the Search for the Holy Grail

CH2D+, the singly deuterated counterpart of CH3+, offers an alternative way to mediate formation of deuterated species at temperatures of several tens of K, as compared to the release of deuterated species from grains. We report a longstanding observational search for this molecular ion, whose rotational spectroscopy is not yet completely secure. We summarize the main spectroscopic properties of this molecule and discuss the chemical network leading to the formation of CH2D+, with explicit account of the ortho/para forms of H2, H3+ and CH3+. Astrochemical models support the presence of this molecular ion in moderately warm environments at a marginal level.

preprint2013arXiv

Detection of the Ammonium Ion in Space

We report on the detection of a narrow feature at 262816.73 MHz towards Orion and the cold prestellar core B1-bS, that we attribute to the 1(0)-0(0) line of the deuterated Ammonium ion, NH3D+. The observations were performed with the IRAM 30m radio telescope. The carrier has to be a light molecular species as it is the only feature detected over 3.6 GHz of bandwidth. The hyperfine structure is not resolved indicating a very low value for the electric quadrupolar coupling constant of Nitrogen which is expected for NH3D+ as the electric field over the N nucleus is practically zero. Moreover, the feature is right at the predicted frequency for the 1(0)-0(0) transition of the Ammonium ion, 262817(6) MHz (3sigma), using rotational constants derived from new infrared data obtained in our laboratory in Madrid. The estimated column density is 1.1(0.2)e12 cm-2. Assuming a deuterium enhancement similar to that of NH2D, we derive N(NH4+) sim 2.6e13 cm-2, i.e., an abundance for Ammonium of a few 1e(-11).

preprint2013arXiv

Hydrogen chloride in diffuse interstellar clouds along the line of sight to W31C (G10.6-0.4)

We report the detection of hydrogen chloride, HCl, in diffuse molecular clouds on the line of sight towards the star-forming region W31C (G10.6-0.4). The J = 1-0 lines of the two stable HCl isotopologues, H35Cl and H37Cl, are observed using the 1b receiver of the Heterodyne Instrument for the Far-Infrared (HIFI) aboard the Herschel Space Observatory. The HCl line is detected in absorption, over a wide range of velocities associated with diffuse clouds along the line of sight to W31C. The analysis of the absorption strength yields a total HCl column density of few 10^13 cm^-2, implying that HCl accounts for ~0.6 % of the total gasphase chlorine, which exceeds by a factor of ~6 the theoretical model predictions. This result is comparable to those obtained from the chemically-related species H2Cl+ and HCl+, for which large column densities have also been reported on the same line of sight. The source of discrepancy between models and observations is still unknown; however, the detection of these Cl-bearing molecules, provides key constraints for the chlorine chemistry in the diffuse gas.

preprint2013arXiv

Nitrogen isotopic ratios in Barnard 1: a consistent study of the N2H+, NH3, CN, HCN and HNC isotopologues

The 15N isotopologue abundance ratio measured today in different bodies of the solar system is thought to be connected to 15N-fractionation effects that would have occured in the protosolar nebula. The present study aims at putting constraints on the degree of 15N-fractionation that occurs during the prestellar phase, through observations of D, 13C and 15N-substituted isotopologues towards B1b. Both molecules from the nitrogen hydride family, i.e. N2H+ and NH3, and from the nitrile family, i.e. HCN, HNC and CN, are considered in the analysis. As a first step, we model the continuum emission in order to derive the physical structure of the cloud, i.e. gas temperature and H2 density. These parameters are subsequently used as an input in a non-local radiative transfer model to infer the radial abundances profiles of the various molecules. Our modeling shows that all the molecules are affected by depletion onto dust grains, in the region that encompasses the B1-bS and B1-bN cores. While high levels of deuterium fractionation are derived, we conclude that no fractionation occurs in the case of the nitrogen chemistry. Independently of the chemical family, the molecular abundances are consistent with 14N/15N~300, a value representative of the elemental atomic abundances of the parental gas. The inefficiency of the 15N-fractionation effects in the B1b region can be linked to the relatively high gas temperature ~17K which is representative of the innermost part of the cloud. Since this region shows signs of depletion onto dust grains, we can not exclude the possibility that the molecules were previously enriched in 15N, earlier in the B1b history, and that such an enrichment could have been incorporated into the ice mantles. It is thus necessary to repeat this kind of study in colder sources to test such a possibility.

preprint2013arXiv

The IRAM-30m line survey of the Horsehead PDR: III. High abundance of complex (iso-)nitrile molecules in UV-illuminated gas

Complex (iso-)nitrile molecules, such as CH3CN and HC3N, are relatively easily detected in our Galaxy and in other galaxies. We constrain their chemistry through observations of two positions in the Horsehead edge: the photo-dissociation region (PDR) and the dense, cold, and UV-shielded core just behind it. We systematically searched for lines of CH3CN, HC3N, C3N, and some of their isomers in our sensitive unbiased line survey at 3, 2, and 1mm. We derived column densities and abundances through Bayesian analysis using a large velocity gradient radiative transfer model. We report the first clear detection of CH3NC at millimeter wavelength. We detected 17 lines of CH3CN at the PDR and 6 at the dense core position, and we resolved its hyperfine structure for 3 lines. We detected 4 lines of HC3N, and C3N is clearly detected at the PDR position. We computed new electron collisional rate coefficients for CH3CN, and we found that including electron excitation reduces the derived column density by 40% at the PDR position. While CH3CN is 30 times more abundant in the PDR than in the dense core, HC3N has similar abundance at both positions. The isomeric ratio CH3NC/CH3CN is 0.15+-0.02. In the case of CH3CN, pure gas phase chemistry cannot reproduce the amount of CH3CN observed in the UV-illuminated gas. We propose that CH3CN gas phase abundance is enhanced when ice mantles of grains are destroyed through photo-desorption or thermal-evaporation in PDRs, and through sputtering in shocks. (abridged)

preprint2013arXiv

The IRAM-30m line survey of the Horsehead PDR: IV. Comparative chemistry of H2CO and CH3OH

Aims. We investigate the dominant formation mechanism of H2CO and CH3OH in the Horsehead PDR and its associated dense core. Methods. We performed deep integrations of several H2CO and CH3OH lines at two positions in the Horsehead, namely the PDR and dense core, with the IRAM-30m telescope. In addition, we observed one H2CO higher frequency line with the CSO telescope at both positions. We determine the H2CO and CH3OH column densities and abundances from the single-dish observations complemented with IRAM-PdBI high-angular resolution maps (6") of both species. We compare the observed abundances with PDR models including either pure gas-phase chemistry or both gas-phase and grain surface chemistry. Results. We derive CH3OH abundances relative to total number of hydrogen atoms of ~1.2e-10 and ~2.3e-10 in the PDR and dense core positions, respectively. These abundances are similar to the inferred H2CO abundance in both positions (~2e-10). We find an abundance ratio H2CO/CH3OH of ~2 in the PDR and ~1 in the dense core. Pure gas-phase models cannot reproduce the observed abundances of either H2CO or CH3OH at the PDR position. Both species are therefore formed on the surface of dust grains and are subsequently photodesorbed into the gas-phase at this position. At the dense core, on the other hand, photodesorption of ices is needed to explain the observed abundance of CH3OH, while a pure gas-phase model can reproduce the observed H2CO abundance. The high-resolution observations show that CH3OH is depleted onto grains at the dense core. CH3OH is thus present in an envelope around this position, while H2CO is present in both the envelope and the dense core itself. Conclusions. Photodesorption is an efficient mechanism to release complex molecules in low FUV-illuminated PDRs, where thermal desorption of ice mantles is ineffective.

preprint2012arXiv

Multi-line detection of O2 toward rho Oph A

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

preprint2012arXiv

The hyperfine structure in the rotational spectrum of CF+

Context. CF+ has recently been detected in the Horsehead and Orion Bar photo-dissociation regions. The J=1-0 line in the Horsehead is double-peaked in contrast to other millimeter lines. The origin of this double-peak profile may be kinematic or spectroscopic. Aims. We investigate the effect of hyperfine interactions due to the fluorine nucleus in CF+ on the rotational transitions. Methods. We compute the fluorine spin rotation constant of CF+ using high-level quantum chemical methods and determine the relative positions and intensities of each hyperfine component. This information is used to fit the theoretical hyperfine components to the observed CF+ line profiles, thereby employing the hyperfine fitting method in GILDAS. Results. The fluorine spin rotation constant of CF+ is 229.2 kHz. This way, the double-peaked CF+ line profiles are well fitted by the hyperfine components predicted by the calculations. The unusually large hyperfine splitting of the CF+ line therefore explains the shape of the lines detected in the Horsehead nebula, without invoking intricate kinematics in the UV-illuminated gas.

preprint2012arXiv

The IRAM-30m line survey of the Horsehead PDR: I. CF+ as a tracer of C+ and a measure of the Fluorine abundance

C+ is a key species in the interstellar medium but its 158 μm fine structure line cannot be observed from ground-based telescopes. Current models of fluorine chemistry predict that CF+ is the second most important fluorine reservoir, in regions where C+ is abundant. We detected the J = 1-0 and J = 2-1 rotational lines of CF+ with high signal-to-noise ratio towards the PDR and dense core positions in the Horsehead. Using a rotational diagram analysis, we derive a column density of N(CF+) = (1.5 - 2.0) \times 10^12 cm^-2. Because of the simple fluorine chemistry, the CF+ column density is proportional to the fluorine abundance. We thus infer the fluorine gas-phase abundance to be F/H = (0.6 - 1.5) \times 10^-8. Photochemical models indicate that CF+ is found in the layers where C+ is abundant. The emission arises in the UV illuminated skin of the nebula, tracing the outermost cloud layers. Indeed, CF+ and C+ are the only species observed to date in the Horsehead with a double peaked line profile caused by kinematics. We therefore propose that CF+, which is detectable from the ground, can be used as a proxy of the C+ layers.

preprint2012arXiv

The IRAM-30m line survey of the Horsehead PDR: II. First detection of the l-C3H+ hydrocarbon cation

We present the first detection of the l-C3H+ hydrocarbon in the interstellar medium. The Horsehead WHISPER project, a millimeter unbiased line survey at two positions, namely the photo-dissociation region (PDR) and the nearby shielded core, revealed a consistent set of eight unidentified lines toward the PDR position. Six of them are detected with a signal-to-noise ratio from 6 to 19, while the two last ones are tentatively detected. Mostly noise appears at the same frequency toward the dense core, located less than 40" away. We simultaneously fit 1) the rotational and centrifugal distortion constants of a linear rotor, and 2) the Gaussian line shapes located at the eight predicted frequencies. The observed lines can be accurately fitted with a linear rotor model, implying a 1Sigma ground electronic state. The deduced rotational constant value is Be= 11244.9512 +/- 0.0015 MHz, close to that of l-C3H. We thus associate the lines to the l-C3H+ hydrocarbon cation, which enables us to constrain the chemistry of small hydrocarbons. A rotational diagram is then used to infer the excitation temperature and the column density. We finally compare the abundance to the results of the Meudon PDR photochemical model.

preprint2011arXiv

H2CO in the Horsehead PDR: Photo-desorption of dust grain ice mantles

Aims: For the first time we investigate the role of the grain surface chemistry in the Horsehead Photo-dissociation region (PDR). Methods: We performed deep observations of several H2CO rotational lines toward the PDR and its associated dense-core in the Horsehead nebula, where the dust is cold ($\Td \simeq 20-30$ K). We complemented these observations with a map of the p-H2CO $~3_{03}-2_{02}$ line at 218.2 GHz (with 12" angular resolution). We determine the H2CO abundances using a detailed radiative transfer analysis and compare these results with PDR models that include either pure gas-phase chemistry or both gas-phase and grain surface chemistry. Results: The H2CO abundances ($\simeq 2-3 \times 10^{-10}$) with respect to H-nuclei are similar in the PDR and dense-core. In the dense-core the pure gas-phase chemistry model reproduces the observed H2CO abundance. Thus, surface processes do not contribute significantly to the gas-phase H2CO abundance in the core. In contrast, the formation of H2CO on the surface of dust grains and subsequent photo-desorption into the gas-phase are needed in the PDR to explain the observed gas-phase H2CO abundance, because the gas-phase chemistry alone does not produce enough H2CO. The assignments of different formation routes are strengthen by the different measured ortho-to-para ratio of H2CO: the dense-core displays the equilibrium value ($\sim3$) while the PDR displays an out-of-equilibrium value ($\sim2$). Conclusions: Photo-desorption of H2CO ices is an efficient mechanism to release a significant amount of gas-phase H2CO into the Horsehead PDR.

preprint2011arXiv

Observing a column-dependent zeta in dense interstellar sources: the case of the Horsehead Nebula

Context: Observations of small carbon-bearing molecules such as CCH, C4H, c-C3H2, and HCO in the Horsehead Nebula have shown these species to have higher abundances towards the edge of the source than towards the center. Aims: Given the determination of a wide range of values for zeta (s-1), the total ionization rate of hydrogen atoms, and the proposal of a column-dependent zeta(N_H), where N_H is the total column of hydrogen nuclei, we desire to determine if the effects of zeta(N_H) in a single object with spatial variation can be observable. We chose the Horsehead Nebula because of its geometry and high density. Method: We model the Horsehead Nebula as a near edge-on photon dominated region (PDR), using several choices for zeta, both constant and as a function of column. The column-dependent zeta functions are determined by a Monte Carlo model of cosmic ray penetration, using a steep power-law spectrum and accounting for ionization and magnetic field effects. We consider a case with low-metal elemental abundances as well as a sulfur-rich case. Results: We show that use of a column-dependent zeta(N_H) of 5(-15) s-1 at the surface and 7.5(-16) s-1 at Av = 10 on balance improves agreement between measured and theoretical molecular abundances, compared with constant values of zeta.

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

Nitrogen Isotopic Fractionation in Interstellar Ammonia

Using the Green Bank Telescope (GBT), we have obtained accurate measurements of the $^{14}$N/$^{15}$N isotopic ratio in ammonia in two nearby cold, dense molecular clouds, Barnard~1 and NGC 1333. The $^{14}$N/$^{15}$N ratio in Barnard~1, $334 \pm 50$ (3$σ$), is particularly well constrained and falls in between the local interstellar medium/proto-solar value of $\sim 450$ and the terrestrial atmospheric value of 272. The NGC 1333 measurement is consistent with the Barnard~1 result, but has a larger uncertainty. We do not see evidence for the very high $^{15}$N enhancements seen in cometary CN. Sensitive observations of a larger, carefully selected sample of prestellar cores with varying temperatures and gas densities can significantly improve our understanding of the nitrogen fractionation in the local interstellar medium and its relation to the isotopic ratios measured in various solar system reservoirs.

preprint2010arXiv

Reaction Networks For Interstellar Chemical Modelling: Improvements and Challenges

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

preprint2010arXiv

Statistical universal branching ratios for cosmic ray dissociation, photodissociation, and dissociative recombination of the C(n=2-10), C(n=2-4)H and C3H2 neutral and cationic species

Fragmentation branching ratios of electronically excited molecular species are of first importance for the modeling of gas phase interstellar chemistry. Despite experimental and theoretical efforts that have been done during the last two decades there is still a strong lack of detailed information on those quantities for many molecules such as Cn, CnH or C3H2. Our aim is to provide astrochemical databases with more realistic branching ratios for Cn (n=2 to 10), CnH (n=2 to 4), and C3H2 molecules that are electronically excited either by dissociative recombination, photodissociation, or cosmic ray processes, when no detailed calculations or measurements exist in literature. High velocity collision in an inverse kinematics scheme was used to measure the complete fragmentation pattern of electronically excited Cn (n=2 to 10), CnH (n=2 to 4), and C3H2 molecules. Branching ratios of dissociation where deduced from those experiments. The full set of branching ratios was used as a new input in chemical models and branching ratio modification effects observed in astrochemical networks that describe the dense cold Taurus Molecular Cloud-1 and the photon dominated Horse Head region. The comparison between the branching ratios obtained in this work and other types of experiments showed a good agreement. It was interpreted as the signature of a statistical behavior of the fragmentation. The branching ratios we obtained lead to an increase of the C3 production together with a larger dispersion of the daughter fragments. The introduction of these new values in the photon dominated region model of the Horse Head nebula increases the abundance of C3 and C3H, but reduces the abundances of the larger Cn and hydrocarbons at a visual extinction Av smaller than 4.

preprint2010arXiv

The chemistry of vibrationally excited H2 in the interstellar medium

The internal energy available in vibrationally excited H2 molecules can be used to overcome or diminish the activation barrier of various chemical reactions of interest for molecular astrophysics. In this article we investigate in detail the impact on the chemical composition of interstellar clouds of the reactions of vibrationally excited H2 with C+, He+, O, OH, and CN, based on the available chemical kinetics data. It is found that the reaction of H2 (v>0) and C+ has a profound impact on the abundances of some molecules, especially CH+, which is a direct product and is readily formed in astronomical regions with fractional abundances of vibrationally excited H2, relative to ground state H2, in excess of 10^(-6), independently of whether the gas is hot or not. The effects of these reactions on the chemical composition of the diffuse clouds zeta Oph and HD 34078, the dense PDR Orion Bar, the planetary nebula NGC 7027, and the circumstellar disk around the B9 star HD 176386 are investigated through PDR models. We find that formation of CH+ is especially favored in dense and highly FUV illuminated regions such as the Orion Bar and the planetary nebula NGC 7027, where column densities in excess of 10^(13) cm^(-2) are predicted. In diffuse clouds, however, this mechanism is found to be not efficient enough to form CH+ with a column density close to the values derived from astronomical observations.

preprint2009arXiv

CO emission and variable CH and CH+ absorption towards HD34078: evidence for a nascent bow shock ?

The runaway star HD34078, initially selected to investigate small scale structure in a foreground diffuse cloud has been shown to be surrounded by highly excited H2. We first search for an association between the foreground cloud and HD34078. Second, we extend previous investigations of temporal absorption line variations (CH, CH+, H2) in order to better characterize them. We have mapped the CO(2-1) emission at 12 arcsec resolution around HD34078's position, using the 30 m IRAM antenna. The follow-up of CH and CH+ absorption lines has been extended over 5 more years. In parallel, CH absorption towards the reddened star Zeta Per have been monitored to check the homogeneity of our measurements. Three more FUSE spectra have been obtained to search for N(H2) variations. CO observations show a pronounced maximum near HD34078's position, clearly indicating that the star and diffuse cloud are associated. The optical spectra confirm the reality of strong, rapid and correlated CH and CH+ fluctuations. On the other hand, N(H2, J=0) has varied by less than 5 % over 4 years. We also discard N(CH) variations towards Zeta Per at scales less than 20 AU. Observational constraints from this work and from 24 micron dust emission appear to be consistent with H2 excitation but inconsistent with steady-state bow shock models and rather suggest that the shell of compressed gas surrounding HD34078, is seen at an early stage of the interaction. The CH and CH+ time variations as well as their large abundances are likely due to chemical structure in the shocked gas layer located at the stellar wind/ambient cloud interface. Finally, the lack of variations for both N(H2, J=0) towards HD34078 and N(CH) towards Zeta Per suggests that quiescent molecular gas is not subject to pronounced small-scale structure.

preprint2009arXiv

Deuterium chemistry in the Orion Bar PDR - "warm" chemistry starring CH2D+

High levels of deuterium fractionation in gas-phase molecules are usually associated with cold regions, such as prestellar cores. Significant fractionation ratios are also observed in hot environments such as hot cores or hot corinos, where they are believed to be produced by the evaporation of the icy mantles surrounding dust grains, and thus are remnants of a previous cold (either gas-phase or grain surface) chemistry. The recent detection of DCN towards the Orion Bar, in a clump at a characteristic temperature of 70K, has shown that high deuterium fractionation can also be detected in PDRs. The Orion Bar clumps thus appear as a good environment for the observational study of deuterium fractionation in luke-warm gas, allowing to validate chemistry models in a different temperature range, where dominating fractionation processes are predicted to be different than in cold gas (< 20K). We aimed at studying observationally in detail the chemistry at work in the Orion Bar PDR, to understand if DCN is produced by ice mantle evaporation, or is the result of warm gas-phase chemistry, involving the CH2D+ precursor ion (which survives higher temperatures than the usual H2D+ precursor). Using the APEX and the IRAM 30m telescopes, we targetted selected deuterated species towards two clumps in the Orion Bar. We confirmed the detection of DCN and detected two new deuterated molecules (DCO+ and HDCO) towards one clump in the Orion Bar PDR. Significant deuterium fractionations are found for HCN and H2CO, but a low fractionation in HCO+. We also give upper limits for other molecules relevant for the deuterium chemistry. (...) We show evidence that warm deuterium chemistry driven by CH2D+ is at work in the clumps.

preprint2008arXiv

Molecular cloud chemistry and the importance of dielectronic recombination

Dielectronic recombination (DR) of singly charged ions is a reaction pathway that is commonly neglected in chemical models of molecular clouds. In this study we include state-of-the-art DR data for He$^+$, C$^+$, N$^+$, O$^+$, Na$^+$, and Mg$^+$ in chemical models used to simulate dense molecular clouds, protostars, and diffuse molecular clouds. We also update the radiative recombination (RR) rate coefficients for H$^+$, He$^+$, C$^+$, N$^+$, O$^+$, Na$^+$, and Mg$^+$ to the current state-of-the-art values. The new RR data has little effect on the models. However, the inclusion of DR results in significant differences in gas-grain models of dense, cold molecular clouds for the evolution of a number of surface and gas-phase species. We find differences of a factor of 2 in the abundance for 74 of the 655 species at times of $10^4$--$10^6$ years in this model when we include DR. Of these 74 species, 16 have at least a factor of 10 difference in abundance. We find the largest differences for species formed on the surface of dust grains. These differences are due primarily to the addition of C$^+$ DR, which increases the neutral C abundance, thereby enhancing the accretion of C onto dust. These results may be important for the warm-up phase of molecular clouds when surface species are desorbed into the gas phase. We also note that no reliable state-of-the-art RR or DR data exist for Si$^+$, P$^+$, S$^+$, Cl$^+$, and Fe$^+$. Modern calculations for these ions are needed to better constrain molecular cloud models.

preprint2007arXiv

A PDR-Code Comparison Study

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. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. 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 selfconsistently. 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, 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.