Source author record

A. Faure

A. Faure 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

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

30 published item(s)

preprint2022arXiv

Absolute measurements of state-to-state rotational energy transfer between CO and H2 at interstellar temperatures

Experimental measurements and theoretical calculations of state-to-state rate coefficients for rotational energy transfer of CO in collision with H$_2$ are reported at the very low temperatures prevailing in dense interstellar clouds (5 - 20 K). Detailed agreement between quantum state-selected experiments performed in cold supersonic flows using time-resolved infrared - vacuum-ultraviolet double resonance spectroscopy and close-coupling quantum scattering calculations confirms the validity of the calculations for collisions between the two most abundant molecules in the interstellar medium.

preprint2022arXiv

Collisional excitation and non-LTE modelling of interstellar chiral propylene oxide

The first set of theoretical cross sections for propylene oxide (CH3CHCH2O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H2 are deduced in the temperature range 5-30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made towards the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (~ 10 K) and translucent (nH ~ 2000 cm-3) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is Ntot ~ 3e12 cm-2, corresponding to a fractional abundance relative to total hydrogen of ~ 2.5e-11. The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.

preprint2022arXiv

The ortho-to-para ratio of water in interstellar clouds

The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involving up to seven protons. The abundances and ortho-to-para ratios (OPRs) of gas-phase water and water ions (H$_2$O$^+$ and H$_3$O$^+$) are computed under the steady-state conditions representative of a dark molecular cloud and during the early phase of gravitational collapse of a prestellar core. The model incorporates the freezing of the molecules on to grains, simple grain surface chemistry and cosmic-ray induced and direct desorption of ices. The predicted OPRs are found to deviate significantly from both thermal and statistical values and to be independent of temperature below $\sim $30~K. The OPR of H$_2$O is shown to lie between 1.5 and 2.6, depending on the spin-state of H$_2$, in good agreement with values derived in translucent clouds with relatively high extinction. In the prestellar core collapse calculations, the OPR of H$_2$O is shown to reach the statistical value of 3 in regions with severe depletion ($n_{\rm H}> 10^7$~cm$^{-3}$). We conclude that a low water OPR ($\lesssim 2.5$) is consistent with gas-phase ion-neutral chemistry and reflects a gas with OPR(H$_2)\lesssim 1$. Available OPR measurements in protoplanetary disks and comets are finally discussed.

preprint2021arXiv

Sulfur gas-phase abundance in dense cores

The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. However, atomic sulfur, which could be the main gas-phase carrier, cannot be observed directly in cold cores. We present measurements of the nitrogen sulfide (NS) radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that $n(\rm S)/n_{\rm H}$ is close, or even equal, to the cosmic abundance of sulfur, 14$\times 10^{-6}$, demonstrating that sulfur is not depleted and is atomic, which is in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime also need to be explored from their chemical and dynamical perspectives.

preprint2020arXiv

Ammonium salts are a reservoir of nitrogen on a cometary nucleus and possibly on some asteroids

The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.

preprint2020arXiv

The effect of CO-H2O collisions in the rotational excitation of cometary CO

We present the first accurate rate coefficients for the rotational excitation of CO by H2O in the kinetic temperature range 5-100 K. The statistical adiabatic channel method (SACM) is combined with a high-level rigid-rotor CO-H2O intermolecular potential surface. Transitions among the first 11 rotational levels of CO and the first 8 rotational levels of both para-H2O and ortho-H2O are considered. Our rate coefficients are compared to previous data from the literature and they are also incorporated in a simple non-LTE model of cometary coma including collision-induced transitions, solar radiative pumping and radiative decay. We find that the uncertainties in the collision data have significant influence on the CO population distribution for H2O densities in the range 10^3-10^8 cm^-3. We also show that the rotational distribution of H2O plays an important role in CO excitation (owing to correlated energy transfer in both CO and H2O), while the impact of the ortho-to-para ratio of H2O is found to be negligible.

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

Collisional excitation of HC3N by para- and ortho-H2

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.

preprint2016arXiv

Detection of protonated formaldehyde in the prestellar core L1689B

Complex organic molecules (COMs) are detected in many regions of the interstellar medium, including prestellar cores. However, their formation mechanisms in cold (~10 K) cores remain to this date poorly understood. The formyl radical HCO is an important candidate precursor for several O-bearing terrestrial COMs in cores, as an abundant building block of many of these molecules. Several chemical routes have been proposed to account for its formation, both on grain surfaces, as an incompletely hydrogenated product of H addition to frozen-out CO molecules, or in the gas phase, either the product of the reaction between H2CO and a radical, or as a product of dissociative recombination of protonated formaldehyde H2COH+. The detection and abundance determination of H2COH+, if present, could provide clues as to whether this latter scenario might apply. We searched for protonated formaldehyde H2COH+ in the prestellar core L1689B using the IRAM 30m telescope. The H2COH+ ion is unambiguously detected, for the first time in a cold (~10 K) source. The derived abundance agrees with a scenario in which the formation of H2COH+ results from the protonation of formaldehyde. We use this abundance value to constrain the branching ratio of the dissociative recombination of H2COH+ towards the HCO channel to ~10-30%. This value could however be smaller if HCO can be efficiently formed from gas-phase neutral-neutral reactions, and we stress the need for laboratory measurements of the rate constants of these reactions at 10 K. Given the experimental difficulties in measuring branching ratios experimentally, observations can bring valuable constraints on these values, and provide a useful input for chemical networks.

preprint2016arXiv

History of the solar-type protostar IRAS16293-2422 as told by the cyanopolyynes

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.

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

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

C$_2$H observations toward the Orion Bar

C$_2$H is one of the first radicals to be detected in the interstellar medium. Its higher rotational transitions have recently become available with the Herschel Space Observatory. We aim to constrain the physical parameters of the C$_2$H emitting gas toward the Orion Bar. We analyse the C$_2$H line intensities measured toward the Orion Bar CO$^+$ Peak and Herschel/HIFI maps of C$_2$H, CH, and HCO$^+$, and a NANTEN map of [CI]. We interpret the observed C$_2$H emission using radiative transfer and PDR models. Five rotational transitions of C$_2$H have been detected in the HIFI frequency range toward the CO$^+$ peak. A single component rotational diagram gives a rotation temperature of ~64 K and a beam-averaged C$_2$H column density of 4$\times$10$^{13}$ cm$^{-2}$. The measured transitions cannot be explained by any single parameter model. According to a non-LTE model, most of the C$_2$H column density produces the lower-$N$ C$_2$H transitions and traces a warm ($T_{\rm{kin}}$ ~ 100-150 K) and dense ($n$(H$_2$)~10$^5$-10$^6$ cm$^{-3}$) gas. A small fraction of the C$_2$H column density is required to reproduce the intensity of the highest-$N$ transitions ($N$=9-8 and N=10-9) originating from a high density ($n$(H$_2$)~5$\times$10$^6$ cm$^{-3}$) hot ($T_{\rm{kin}}$ ~ 400 K) gas. The total beam-averaged C$_2$H column density in the model is 10$^{14}$ cm$^{-2}$. Both the non-LTE radiative transfer model and a simple PDR model representing the Orion Bar with a plane-parallel slab of gas and dust suggest, that C$_2$H cannot be described by a single pressure component, unlike the reactive ion CH$^+$, which was previously analysed toward the Orion Bar CO$^+$ peak. The physical parameters traced by the higher rotational transitions ($N$=6-5,...,10-9) of C$_2$H may be consistent with the edges of dense clumps exposed to UV radiation near the ionization front of the Orion Bar.

preprint2015arXiv

Hydrogen isotope exchanges between water and methanol in interstellar ices

The deuterium fractionation of gas-phase molecules in hot cores is believed to reflect the composition of interstellar ices. The deuteration of methanol is a major puzzle, however, because the isotopologue ratio [CH2DOH]/[CH3OD], which is predicted to be equal to 3 by standard grain chemistry models, is much larger (~20) in low-mass hot corinos and significantly lower (~1) in high-mass hot cores. This dichotomy in methanol deuteration between low-mass and massive protostars is currently not understood. In this study, we report a simplified rate equation model of the deuterium chemistry occurring in the icy mantles of interstellar grains. We apply this model to the chemistry of hot corinos and hot cores, with IRAS 16293-2422 and the Orion~KL Compact Ridge as prototypes, respectively. The chemistry is based on a statistical initial deuteration at low temperature followed by a warm-up phase during which thermal hydrogen/deuterium (H/D) exchanges occur between water and methanol. The exchange kinetics is incorporated using laboratory data. The [CH2DOH]/[CH3OD] ratio is found to scale inversely with the D/H ratio of water, owing to the H/D exchange equilibrium between the hydroxyl (-OH) functional groups of methanol and water. Our model is able to reproduce the observed [CH2DOH]/[CH3OD] ratios provided that the primitive fractionation of water ice [HDO]/[H2O] is ~ 2% in IRAS 16293-2422 and ~0.6% in Orion~KL. We conclude that the molecular D/H ratios measured in hot cores may not be representative of the original mantles because molecules with exchangeable deuterium atoms can equilibrate with water ice during the warm-up phase.

preprint2015arXiv

Ortho-to-para ratio of NH2. Herschel-HIFI observations of ortho- and para-NH2 rotational transitions towards W31C, W49N, W51 and G34.3+0.1

We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores and surrounding envelopes in order to investigate the excitation and possible emission of the ground state rotational transitions of ortho-NH2 N_(K_a,K_c} J=1_(1,1) 3/2 - 0_(0,0) 1/2 and para-NH2 2_(1,2) 5/2 - 1_(0,1) 3/2$ used in the OPR calculations. Our best estimate of the average OPR in the envelopes lie above the high temperature limit of three for W49N, specifically 3.5 with formal errors of \pm0.1, but for W31C, W51, and G34.3+0.1 we find lower values of 2.5\pm0.1, 2.7\pm0.1, and 2.3\pm0.1, respectively. Such low values are strictly forbidden in thermodynamical equilibrium since the OPR is expected to increase above three at low temperatures. In the translucent interstellar gas towards W31C, where the excitation effects are low, we find similar values between 2.2\pm0.2 and 2.9\pm0.2. In contrast, we find an OPR of 3.4\pm0.1 in the dense and cold filament connected to W51, and also two lower limits of >4.2 and >5.0 in two other translucent gas components towards W31C and W49N. At low temperatures (T \lesssim 50 K) the OPR of H2 is <10^-1, far lower than the terrestrial laboratory normal value of three. In such a "para-enriched H2" gas, our astrochemical models can reproduce the variations of the observed OPR, both below and above the thermodynamical equilibrium value, by considering nuclear-spin gas-phase chemistry. The models suggest that values below three arise in regions with temperatures >20-25 K, depending on time, and values above three at lower temperatures.

preprint2015arXiv

The 2014 KIDA network for interstellar chemistry

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

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

Collisional excitation of water by hydrogen atoms

We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of both ortho- and para-H$_2$O and for temperatures in the range T = 5-1500 K. These rate coefficients are subsequently compared to the values previously published for the H$_2$O / He and H$_2$O / H$_2$ collisional systems. It is shown that no simple relation exists between the three systems and that specific calculations are thus mandatory.

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

A KInetic Database for Astrochemistry (KIDA)

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.

preprint2012arXiv

On the influence of collisional rate coefficients on the water vapour excitation

Water is a key molecule in many astrophysical studies. Its high dipole moment makes this molecule to be subthermally populated under the typical conditions of most astrophysical objects. This motivated the calculation of various sets of collisional rate coefficients (CRC) for H$_2$O (with He or H$_2$) which are necessary to model its rotational excitation and line emission. We performed accurate non--local non--LTE radiative transfer calculations using different sets of CRC in order to predict the line intensities from transitions that involve the lowest energy levels of H$_2$O (E $<$ 900 K). The results obtained from the different CRC sets are then compared using line intensity ratio statistics. For the whole range of physical conditions considered in this work, we obtain that the intensities based on the quantum and QCT CRC are in good agreement. However, at relatively low H$_2$ volume density ($n$(H$_2$) $<$ 10$^7$ cm$^{-3}$) and low water abundance ($χ$(H$_2$O) $<$ 10$^{-6}$), these physical conditions being relevant to describe most molecular clouds, we find differences in the predicted line intensities of up to a factor of $\sim$ 3 for the bulk of the lines. Most of the recent studies interpreting early Herschel Space Observatory spectra used the QCT CRC. Our results show that although the global conclusions from those studies will not be drastically changed, each case has to be considered individually, since depending on the physical conditions, the use of the QCT CRC may lead to a mis--estimate of the water vapour abundance of up to a factor of $\sim$ 3.

preprint2011arXiv

Detection of HF emission from the Orion Bar

The clumpy density structure of photon-dominated regions is well established, but the physical properties of the clumps and of the surrounding interclump medium are only approximately known. The aim of this paper is to constrain the physical and chemical conditions in the Orion Bar, a prototypical nearby photon-dominated region. We present observations of the HF J=1-0 line, which appears in emission toward the Orion Bar, and compare the brightness of the line to non-LTE radiative transfer calculations. The large width of the HF line suggests an origin of the emission in the interclump gas, but collisional excitation by H2 in the interclump gas underpredicts the observed line intensity by factors of 3-5. In contrast, an origin of the line in the dense clumps requires a density of ~10^9 cm^-3, 10-100 times higher than previous estimates, which is unlikely. However, electron impact excitation reproduces our observations for T = 100 K and n(e) = 10 cm^-3, as expected for the interclump gas. We conclude that HF emission is a signpost of molecular gas with a high electron density. Similar conditions may apply to active galactic nuclei where HF also appears in emission.

preprint2011arXiv

First detection of Hydrogen Chloride towards protostellar shocks

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.

preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

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.

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.

preprint2010arXiv

The rotational excitation of HCN and HNC by He: New insights on the HCN/HNC abundance ratio in molecular clouds

Modeling of molecular emission from interstellar clouds requires the calculation of rates for excitation by collisions with the most abundant species. The present paper focuses on the calculation of rate coefficients for rotational excitation of the HCN and HNC molecules in their ground vibrational state in collision with He. The calculations are based on new two-dimensional potential energy surfaces obtained from highly correlated \textit{ab initio} calculations. Calculations of pure rotational (de)excitation cross sections of HCN and HNC by He were performed using the essentially exact close-coupling method. Cross sections for transitions among the 8 first rotational levels of HCN and HNC were calculated for kinetic energies up to 1000 cm$^{-1}$. These cross sections were used to determine collisional rate constants for temperatures ranging from 5 K to 100 K. A propensity for even $Δj$ transitions is observed in the case of HCN--He collisions whereas a propensity for odd $Δj$ transitions is observed in the case of HNC--He collisions. The consequences for astrophysical models are evaluated and it is shown that the use of HCN rate coefficients to interpret HNC observations can lead to significant inaccuracies in the determination of the HNC abundance, in particular in cold dark clouds for which the new HNC rates show that the $j=1-0$ line of this species will be more easily excited by collisions than HCN. An important result of the new HNC-He rates is that the HNC/HCN abundance ratio derived from observations in cold clouds has to be revised from $>$1 to $\simeq$1, in good agreement with detailed chemical models available in the literature.

preprint2009arXiv

Physical structure and water line spectrum predictions of the Intermediate Mass protostar OMC2-FIR4

Aims. Intermediate Mass (IM) stars are an important component of our Galaxy, as they significantly contribute to the interstellar FUV field and, consequently, play an important role in the energy balance of the ISM. Despite their importance, very little is known about their formation process and only a few studies have been devoted to characterize the first phases in the evolution of intermediate mass protostars. Here we consider in great detail the case of the brightest and closest known young IM protostar: FIR4 in the OMC2 component of the Orion molecular cloud complex. Methods. We analyzed the available continuum emission (maps and SED) through one-dimensional dust radiative transfer calculations. We ran large grids of models to find the envelope model that best fits the data. The derived dust density and temperature profiles have been then used to compute the gas temperature profile, equating gas cooling and heating terms across the envelope. Last, we computed the water line spectrum for various possible values of water abundance. Results. The luminosity of FIR4 has been reevaluated to 1000 Lo, making FIR4 definitively an Intermediate Mass protostar. The envelope surrounding FIR4 has a relatively shallow density power law index, ~ 0.6. The most surprising result is that the gas and dust are thermally decoupled in the interior of the envelope, where the dust ices sublimate at 100 K. This has important consequences in the interpretation of the line data. We provide the predictions for the water spectrum, and discuss in detail the lines which will be observed by the Herschel Space Observatory .

preprint2006arXiv

First Measurements of the Electron Density Enhancements Expected in C-shocks

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.