Chemistry of dark clouds: databases, networks, and models

A Globally Accurate Neural Network Potential Energy Surface and Quantum Dynamics Studies on Be+(2S) + H2/D2 → BeH+/BeD+ + H/D Reactions

Chemical reactions between Be+ ions and H2 molecules have significance in the fields of ultracold chemistry and astrophysics, but the corresponding dynamics studies on the ground-state reaction have not been reported because of the lack of a global potential energy surface (PES). Herein, a globally accurate ground-state BeH2+ PES is constructed using the neural network model based on 18,657 ab initio points calculated by the multi-reference configuration interaction method with the aug-cc-PVQZ basis set. On the newly constructed PES, the state-to-state quantum dynamics calculations of the Be+(2S) + H2(v0 = 0; j0 = 0) and Be+(2S) + D2(v0 = 0; j0 = 0) reactions are performed using the time-dependent wave packet method. The calculated results suggest that the two reactions are dominated by the complex-forming mechanism and the direct abstraction process at relatively low and high collision energies, respectively, and the isotope substitution has little effect on the reaction dynamics characteristics. The new PES can be used to further study the reaction dynamics of the BeH2+ system, such as the effects of rovibrational excitations and alignment of reactant molecules, and the present dynamics data could provide an important reference for further experimental studies at a finer level.

Reactivity of the Ethenium Cation (C2H5+) with Ethyne (C2H2): A Combined Experimental and Theoretical Study

The gas-phase reaction between the ethyl cation (C2H5+) and ethyne (C2H2) is re-investigated by measuring absolute reactive cross sections (CSs) and branching ratios (BRs) as a function of collision energy, in the thermal and hyperthermal energy range, via tandem-guided ion beam mass spectrometry under single collision conditions. Dissociative photoionization of C2H5Br using tuneable VUV radiation in the range 10.5-14.0 eV is employed to generate C2H5+, which has also allowed us to explore the impact of increasing (vibrational) excitation on the reactivity. Reactivity experiments are complemented by theoretical calculations, at the G4 level of theory, of the relative energies and structures of the most relevant stationary points on the reactive potential energy hypersurface (PES) and by mass-analyzed ion kinetic energy (MIKE) spectrometry experiments to probe the metastable decomposition from the [C4H7]+ PES and elucidate the underlying reaction mechanisms. Two main product channels have been identified at a centre-of-mass collision energy of ∼0.1 eV: (a) C3H3++CH4, with BR = 0.76±0.05 and (b) C4H5++H2, with BR = 0.22±0.02. A third channel giving C2H3+ in association with C2H4 is shown to emerge at both high internal excitation of C2H5+ and high collision energies. From CS measurements, energy-dependent total rate constants in the range 4.3×10-11-5.2×10-10 cm3·molecule-1·s-1 have been obtained. Theoretical calculations indicate that both channels stem from a common covalently bound intermediate, CH3CH2CHCH+, from which barrierless and exothermic pathways exist for the production of both cyclic c-C3H3+ and linear H2CCCH+ isomers of the main product channel. For the minor C4H5+ product, two isomers are energetically accessible: the three-member cyclic isomer c-C3H2(CH3)+ and the higher energy linear structure CH2CHCCH2+, but their formation requires multiple isomerization steps and passages via transition states lying only 0.11 eV below the reagents' energy, thus explaining the smaller BR. Results have implications for the modeling of hydrocarbon chemistry in the interstellar medium and the atmospheres of planets and satellites as well as in laboratory plasmas (e.g., plasma-enhanced chemical vapor deposition of carbon nanotubes and diamond-like carbon films).

Effect of temperature on the gas-phase reaction of CH3CN with OH radicals: experimental (T = 11.7-177.5 K) and computational (T = 10-400 K) kinetic study

Acetonitrile (CH3CN) is present in the interstellar medium (ISM) in a variety of environments. However, at the ultracold temperatures of the ISM, radical-molecule reactions are not widely investigated because of the experimental handicap of getting organic molecules in the gas phase by conventional techniques. The CRESU (French acronym for Reaction Kinetics in a Uniform Supersonic Flow) technique solves this problem. For this reason, we present in this work the kinetic study of the gas-phase reaction of CH3CN with one of the most ubiquitous radicals, the hydroxyl (OH) radical, as a function of temperature (11.7-177.5 K). The kinetic technique employed to investigate the CH3CN + OH reaction was the pulsed laser photolysis-laser induced fluorescence. The rate coefficient for this reaction k(T) has been observed to drastically increase from 177.5 K to 107.0 K (about 2 orders of magnitude), while the increase in k(T) from 107.0 K to 11.7 K was milder (around 4 times). The temperature dependent expressions for k(T) are provided in the two distinct T-ranges, excluding the upper limit obtained for k(177.5 K): In addition, the rate coefficients estimated by the canonical competitive unified statistical (CCUS) theory show a similar behaviour to the experimental results, when evaluated within the high-pressure limit. This is consistent with the experimentally observed independence of k(T) with total gas density at selected temperatures. Astrochemical networks, such as the KIDA database or UMIST, do not include the CH3CN + OH reaction as a potential depletion process for acetonitrile in the ISM because the current studies predict very low rate coefficients at IS temperatures. According to the model (T = 10 K), the impact of the titled reaction on the abundances of CH3CN appears to be negligible in dark molecular clouds of the ISM (∼1% of the total depletion reactions included in UMIST network). With respect to the potential formation of the CH2CN radical in those environments, even in the most favourable scenario, where this radical could be formed in a 100% yield from the CH3CN + OH reaction, this route would only contribute around 2% to the current assumed formation routes by the UMIST network.

Collisional excitation of HCNH+ by He and H2: New potential energy surfaces and inelastic rate coefficients

Protonated molecules have been increasingly detected in the interstellar medium (ISM), and usually astrochemical models fail at reproducing the abundances derived from observational spectra. Rigorous interpretation of the detected interstellar emission lines requires prior calculations of collisional rate coefficients with H₂ and He, i.e., the most abundant species in the ISM. In this work, we focus on the excitation of HCNH⁺ induced by collision with H₂ and He. Therefore, we first calculate ab initio potential energy surfaces (PESs) using the explicitly correlated and standard coupled cluster method with single, double, and non-iterative triple excitation in conjunction with the augmented-correlation consistent-polarized valence triple zeta basis set. Both the HCNH⁺-H₂ and HCNH⁺-He potentials are characterized by deep global minima of 1426.60 and 271.72 cm^-1, respectively, and large anisotropies. From these PESs, we derive state-to-state inelastic cross sections for the 16 low-lying rotational energy levels of HCNH⁺ using the quantum mechanical close-coupling approach. The differences between cross sections due to ortho- and para-H₂ impacts turn out to be minor. Using a thermal average of these data, we retrieve downward rate coefficients for kinetic temperatures of up to 100 K. As it could be anticipated, differences of up to two orders of magnitude exist between the rate coefficients induced by H₂ and He collisions. We expect that our new collision data will help to improve the disagreement between abundances retrieved from observational spectra and astrochemical models.

The emergence of interstellar molecular complexity explained by interacting networks

Recent years have witnessed the detection of an increasing number of complex organic molecules in interstellar space, some of them being of prebiotic interest. Disentangling the origin of interstellar prebiotic chemistry and its connection to biochemistry and ultimately, to biology is an enormously challenging scientific goal where the application of complexity theory and network science has not been fully exploited. Encouraged by this idea, we present a theoretical and computational framework to model the evolution of simple networked structures toward complexity. In our environment, complex networks represent simplified chemical compounds and interact optimizing the dynamical importance of their nodes. We describe the emergence of a transition from simple networks toward complexity when the parameter representing the environment reaches a critical value. Notably, although our system does not attempt to model the rules of real chemistry nor is dependent on external input data, the results describe the emergence of complexity in the evolution of chemical diversity in the interstellar medium. Furthermore, they reveal an as yet unknown relationship between the abundances of molecules in dark clouds and the potential number of chemical reactions that yield them as products, supporting the ability of the conceptual framework presented here to shed light on real scenarios. Our work reinforces the notion that some of the properties that condition the extremely complex journey from the chemistry in space to prebiotic chemistry and finally, to life could show relatively simple and universal patterns.

Infrared Spectroscopic Study on Swift-Ion Irradiation of Solid N2O-H2O Samples: Synthesis of N-O Bearing Species in Astrophysical Ices

As of early 2022, only six species bearing an N-O bond have been detected toward cold molecular clouds and regions of star formation. It is not clear yet if the small number of N-O bond species found in the interstellar medium so far stems from physical and technological limitations of astronomical detection techniques, or whether in fact molecules that bear an N-O bond are for some reason rare in these objects of the interstellar medium. Astronomical N-O bearing molecules are important because they are part of astrochemical models which propose that they are precursors of hydroxylamine (NH₂OH), a species linked to the formation of prebiotic amino acids in space. The aim of this study is the better understanding of the open question of the interstellar synthesis of N-O bearing species. We have analyzed by infrared spectroscopy an astrophysically relevant polar ice mixture of N₂O:H₂O processed by 90 MeV ¹³⁶Xe²³⁺ ions, which can mimic the physicochemical processes triggered by cosmic rays in water-covered interstellar ice grains. The results show the formation of N₂O₃ and of H₂O₂, but no HN-O species of any kind were detected. Such findings are discussed in light of recent studies from our group and from the literature.

Neural network potential energy surface and quantum dynamics studies for the Ca+(2S) + H2 → CaH+ + H reaction

Reactive collisions of Ca+ ion with H2 molecule play a crucial role in ultracold chemistry, quantum information and other cutting-edge fields, and have been widely concerned experimentally, but the corresponding...

The sulphur depletion problem in molecular clouds: The H2S case

Sulphur is one of the most abundant elements in the Universe and plays a crucial role in biological systems. However, sulphuretted molecules in the ISM are not as abundant as expected and there is no clear answer of where the missing Sulphur is yet. To shed light onto this open question, we focus our attention on the chemistry of H2S, thought to be an important reservoir of Sulphur and formed mainly by grain-phase reactions. To understand the formation of H2S, the growth of ices, and the chemical desorption process, we study the CO, CH3OH, N2H+, and H2S abundances towards Barnard 1b, a Sulphur-rich cloud hosting a first Larson core. We look for correlations between gas-phase abundances of H2S and CH3OH that better constrain the location of the CO snowline in dark cores. Finally, this provides additional data to benchmark models for a deeper insight on the chemical desorption process and its efficiency.

X-ray-Induced CO2 Formation via CO Reaction with TiO2 at Cryogenic Temperature

Cosmic rays, γ-rays, and X-ray-driven reactions on interstellar grains have been the main focus in understanding the extraterrestrial synthesis of complex organic molecules. Herein, using temperature-programmed desorption and X-ray photoelectron spectroscopy, we report for the first time X-ray-induced CO₂ production from a CO-covered rutile TiO₂(110) surface at 110 K, in addition to X-ray-induced CO desorption. The X-ray-induced CO₂ production was found to follow a Mars and van Krevelen (MvK) mechanism via the reaction between adsorbed CO and surface lattice oxygen of TiO₂. Defects of TiO₂ suppress the X-ray-induced CO₂ production but promote the X-ray-induced CO desorption. These findings suggest a novel X-ray-induced CO-to-CO₂ reaction on oxide surfaces at cryogenic temperature likely occurring in the interstellar medium and also warn us to keep in mind the possibility of X-rays inducing chemical reactions during structural characterizations of oxides with X-ray-involved techniques.

The sulphur saga in TMC-1: Discovery of HCSCN and HCSCCH

We report the detection, for the first time in space, of cyano thioformaldehyde (HCSCN) and propynethial (HCSCCH) towards the starless core TMC-1. Cyano thioformaldehyde presents a series of prominent a- and b-type lines, which are the strongest previously unassigned features in our Q-band line survey of TMC-1. Remarkably, HCSCN is four times more abundant than cyano formaldehyde (HCOCN). On the other hand, HCSCCH is five times less abundant than propynal (HCOCCH). Surprisingly, we find an abundance ratio HCSCCH/HCSCN of ∼ 0.25, in contrast with most other ethynyl-cyanide pairs of molecules for which the CCH-bearing species is more abundant than the CN-bearing one. We discuss the formation of these molecules in terms of neutral-neutral reactions of S atoms with CH2CCH and CH2CN radicals as well as of CCH and CN radicals with H2CS. The calculated abundances for the sulphur-bearing species are, however, significantly below the observed values, which points to an underestimation of the abundance of atomic sulphur in the model or to missing formation reactions, such as ion-neutral reactions.

Cumulene carbenes in TMC-1: Astronomical discovery of l-H2C5 ★

We report the first detection in space of the cumulene carbon chain l-H2C5. A total of eleven rotational transitions, with Jup = 7-10 and Ka = 0 and 1, were detected in TMC-1 in the 31.0-50.4 GHz range using the Yebes 40m radio telescope. We derive a column density of (1.8±0.5)×1010 cm-2. In addition, we report observations of other cumulene carbenes detected previously in TMC-1, to compare their abundances with the newly detected cumulene carbene chain. We find that l-H2C5 is ~4.0 times less abundant than the larger cumulene carbene l-H2C6, while it is ~300 and ~500 times less abundant than the shorter chains l-H2C3 and l-H2C4. We discuss the most likely gas-phase chemical routes to these cumulenes in TMC-1 and stress that chemical kinetics studies able to distinguish between different isomers are needed to shed light on the chemistry of C n H2 isomers with n > 3.

Molecular dynamics reveals formation path of benzonitrile and other molecules in conditions relevant to the interstellar medium

Polycyclic aromatic hydrocarbons and polycyclic aromatic nitrogen heterocycles are believed to be widespread in different areas of the interstellar medium. However, the astronomical detection of specific aromatic molecules is extremely challenging. As a result, only a few aromatic molecules have been identified, and very little is known about how they are formed in different areas of the interstellar medium. Recently, McGuire et al. [Science 359, 202-205 (2018)] detected the simple aromatic molecule benzonitrile in Taurus Molecular Cloud-1. Although benzonitrile has been observed, the molecular mechanism for its formation is still unknown. In this study, we use quantum chemistry and ab initio molecular dynamics to model ionization processes of van der Waals clusters containing cyanoacetylene and acetylene molecules. We demonstrate computationally that the clusters' ionization leads to molecular formation. For pure cyanoacetylene clusters, we observe bond formation among two and three monomer units, whereas in mixed clusters, bond formation can also occur in up to four units. We show that the large amount of energy available to the system after ionization can lead to barrier crossing and the formation of complex molecules. Our study reveals the rich chemistry that is observed upon ionization of the clusters, with a wide variety of molecules being formed. Benzonitrile is among the observed molecules, and we study the potential energy path for its formation. These results also offer insights that can guide astronomers in their search for aromatic molecules in the interstellar medium.

Pure hydrocarbon cycles in TMC-1: Discovery of ethynyl cyclopropenylidene, cyclopentadiene and indene

We report the detection for the first time in space of three new pure hydrocarbon cycles in TMC-1: c-C3HCCH (ethynyl cyclopropenylidene), c-C5H6 (cyclopentadiene) and c-C9H8 (indene). We derive a column density of 3.1 × 1011 cm-2 for the former cycle and similar values, in the range (1-2) × 1013 cm-2, for the two latter molecules. This means that cyclopentadiene and indene, in spite of their large size, are exceptionally abundant, only a factor of five less abundant than the ubiquitous cyclic hydrocarbon c-C3H2. The high abundance found for these two hydrocarbon cycles, together with the high abundance previously found for the propargyl radical (CH2CCH) and other hydrocarbons like vinyl and allenyl acetylene (Agúndez et al. 2021; Cernicharo et al. 2021a,b), start to allow us to quantify the abundant content of hydrocarbon rings in cold dark clouds and to identify the intermediate species that are probably behind the in situ bottom-up synthesis of aromatic cycles in these environments. While c-C3HCCH is most likely formed through the reaction between the radical CCH and c-C3H2, the high observed abundances of cyclopentadiene and indene are difficult to explain through currently proposed chemical mechanisms. Further studies are needed to identify how are five- and six-membered rings formed under the cold conditions of clouds like TMC-1.

O-bearing complex organic molecules at the cyanopolyyne peak of TMC-1: detection of C2H3CHO, C2H3OH, HCOOCH3, and CH3OCH3

We report the detection of the oxygen-bearing complex organic molecules propenal (C2H3CHO), vinyl alcohol (C2H3OH), methyl formate (HCOOCH3), and dimethyl ether (CH3OCH3) toward the cyanopolyyne peak of the starless core TMC-1. These molecules are detected through several emission lines in a deep Q-band line survey of TMC-1 carried out with the Yebes 40m telescope. These observations reveal that the cyanopolyyne peak of TMC-1, which is the prototype of cold dark cloud rich in carbon chains, contains also O-bearing complex organic molecules like HCOOCH3 and CH3OCH3, which have been previously seen in a handful of cold interstellar clouds. In addition, this is the first secure detection of C2H3OH in space and the first time that C2H3CHO and C2H3OH are detected in a cold environment, adding new pieces in the puzzle of complex organic molecules in cold sources. We derive column densities of (2.2 ± 0.3) × 1011 cm™2, (2.5 ± 0.5) × 1012 cm-2, (1.1 ± 0.2) × 1012 cm-2, and (2.5 ± 0.7) × 1012 cm-2 for C2H3CHO, C2H3OH, HCOOCH3, and CH3OCH3, respectively. Interestingly, C2H3OH has an abundance similar to that of its well known isomer acetaldehyde (CH3CHO), with C2H3OH/CH3CHO ~ 1 at the cyanopolyyne peak. We discuss potential formation routes to these molecules and recognize that further experimental, theoretical, and astronomical studies are needed to elucidate the true mechanism of formation of these O-bearing complex organic molecules in cold interstellar sources.

TMC-1, the starless core sulfur factory: Discovery of NCS, HCCS, H2CCS, H2CCCS, and C4S and detection of C5S

We report the detection of the sulfur-bearing species NCS, HCCS, H2CCS, H2CCCS, and C4S for the first time in space. These molecules were found towards TMC-1 through the observation of several lines for each species. We also report the detection of C5S for the first time in a cold cloud through the observation of five lines in the 31-50 GHz range. The derived column densities are N(NCS) = (7.8±0.6)×1011 cm-2, N(HCCS) = (6.8±0.6)×1011 cm-2, N(H2CCS) = (7.8±0.8)×1011 cm-2, N(H2CCCS) = (3.7±0.4)×1011 cm-2, N(C4S) = (3.8±0.4)×1010 cm-2, and N(C5S) = (5.0±1.0)×1010 cm-2. The observed abundance ratio between C3S and C4S is 340, that is to say a factor of approximately one hundred larger than the corresponding value for CCS and C3S. The observational results are compared with a state-of-the-art chemical model, which is only partially successful in reproducing the observed abundances. These detections underline the need to improve chemical networks dealing with S-bearing species.

Discovery of CH2CHCCH and detection of HCCN, HC4N, CH3CH2CN, and, tentatively, CH3CH2CCH in TMC-1

We present the discovery in TMC-1 of vinyl acetylene, CH2CHCCH, and the detection, for the first time in a cold dark cloud, of HCCN, HC4N, and CH3CH2CN. A tentative detection of CH3CH2CCH is also reported. The column density of vinyl acetylene is (1.2±0.2)×1013 cm-2, which makes it one of the most abundant closed-shell hydrocarbons detected in TMC-1. Its abundance is only three times lower than that of propylene, CH3CHCH2. The column densities derived for HCCN and HC4N are (4.4±0.4)×1011 cm-2 and (3.7±0.4)×1011 cm-2, respectively. Hence, the HCCN/HC4N abundance ratio is 1.2±0.3. For ethyl cyanide we derive a column density of (1.1 ±0.3)×1011 cm-2. These results are compared with a state-of-the-art chemical model of TMC-1, which is able to account for the observed abundances of these molecules through gas-phase chemical routes.

Discovery of allenyl acetylene, H2CCCHCCH, in TMC-1. A study of the isomers of C5H4

We present the discovery in TMC-1 of allenyl acetylene, H2CCCHCCH, through the observation of nineteen lines with a signal-to-noise ratio ~4-15. For this species, we derived a rotational temperature of 7±1K and a column density of 1.2±0.2×1013 cm-2. The other well known isomer of this molecule, methyl diacetylene (CH3C4H), has also been observed and we derived a similar rotational temperature, Tr=7.0±0.3 K, and a column density for its two states (A and E) of 6.5±0.3×1012 cm-2. Hence, allenyl acetylene and methyl diacetylene have a similar abundance. Remarkably, their abundances are close to that of vinyl acetylene (CH2CHCCH). We also searched for the other isomer of C5H4, HCCCH2CCH (1.4-Pentadiyne), but only a3σ upper limit of 2.5×1012 cm-2 to the column density can be established. These results have been compared to state-of-the-art chemical models for TMC-1, indicating the important role of these hydrocarbons in its chemistry. The rotational parameters of allenyl acetylene have been improved by fitting the existing laboratory data together with the frequencies of the transitions observed in TMC-1.

Discovery of the propargyl radical (CH2CCH) in TMC-1: one of the most abundant radicals ever found and a key species for cyclization to benzene in cold dark clouds

We present the first identification in interstellar space of the propargyl radical (CH2CCH). This species was observed in the cold dark cloud TMC-1 using the Yebes 40m telescope. The six strongest hyperfine components of the 20,2-10,1 rotational transition, lying at 37.46 GHz, were detected with signal-to-noise ratios in the range 4.6-12.3 σ. We derive a column density of 8.7 × 1013 cm-2 for CH2CCH, which translates to a fractional abundance relative to H2 of 8.7 × 10-9. This radical has a similar abundance to methyl acetylene, with an abundance ratio CH2CCH/CH3CCH close to one. The propargyl radical is thus one of the most abundant radicals detected in TMC-1, and it is probably the most abundant organic radical with a certain chemical complexity ever found in a cold dark cloud. We constructed a gas-phase chemical model and find calculated abundances that agree with, or fall two orders of magnitude below, the observed value depending on the poorly constrained low-temperature reactivity of CH2CCH with neutral atoms. According to the chemical model, the propargyl radical is essentially formed by the C + C2H4 reaction and by the dissociative recombination of C3Hn + ions with n = 4-6. The propargyl radical is believed to control the synthesis of the first aromatic ring in combustion processes, and it probably plays a key role in the synthesis of large organic molecules and cyclization processes to benzene in cold dark clouds.

A study of C4H3N isomers in TMC-1: line by line detection of HCCCH2CN

We present Yebes 40m telescope observations of the three most stable C4H3N isomers towards the cyanopolyyne peak of TMC-1. We have detected 13 transitions from CH3C3N (A and E species), 16 lines from CH2CCHCN, and 27 lines (a-type and b-type) from HCCCH2CN. We thus provide a robust confirmation of the detection of HCCCH2CN and CH2CCHCN in space. We have constructed rotational diagrams for the three species, and obtained rotational temperatures between 4-8 K and similar column densities for the three isomers, in the range (1.5-3)×1012 cm-2. Our chemical model provides abundances of the order of the observed ones, although it overestimates the abundance of CH3CCCN and underestimates that of HCCCH2CN. The similarity of the observed abundances of the three isomers suggests a common origin, most probably involving reactions of the radical CN with the unsaturated hydrocarbons methyl acetylene and allene. Studies of reaction kinetics at low temperature and further observations of these molecules in different astronomical sources are needed to draw a clear picture of the chemistry of C4H3N isomers in space.

Space and laboratory discovery of HC3 . ASTRONOMY AND ASTROPHYSICS 2021;646:L3. [PMID: 33824540 PMCID: PMC7610522 DOI: 10.1051/0004-6361/202040013]

We report the detection in TMC-1 of the protonated form of C3S. The discovery of the cation HC3S+ was carried through the observation of four harmonically related lines in the Q band using the Yebes 40m radiotelescope, and is supported by accurate ab initio calculations and laboratory measurements of its rotational spectrum. We derive a column density N(HC3S+) = (2.0 ± 0.5) × 1011 cm-2, which translates to an abundance ratio C3S/HC3S+ of 65 ± 20. This ratio is comparable to the CS/HCS+ ratio (35 ± 8) and is a factor of about ten larger than the C3O/HC3O+ ratio previously found in the same source. However, the abundance ratio HC3O+/HC3S+ is 1.0 ± 0.5, while C3O/C3S is just ~ 0.11. We also searched for protonated C2S in TMC-1, based on ab initio calculations of its spectroscopic parameters, and derive a 3σ upper limit of N(HC2S+)≤ 9×1011 cm-2 and a C2S/HC2S+ ≥ 60. The observational results are compared with a state-of-the-art gas-phase chemical model and conclude that HC3S+ is mostly formed through several pathways: proton transfer to C3S, reaction of S+ with c-C3H2, and reaction between neutral atomic sulfur and the ion C3H+ 3.

Tentative detection of HC5NH+ in TMC-1

Using the Yebes 40m radio telescope, we report the detection of a series of seven lines harmonically related with a rotational constant B 0=1295.81581 ± 0.00026 MHz and a distortion constant D 0 = 27.3 ± 0.5 Hz towards the cold dense cloud TMC-1. Ab initio calculations indicate that the best possible candidates are the cations HC5NH+ and NC4NH+. From a comparison between calculated and observed rotational constants and other arguments based on proton affinities and dipole moments, we conclude that the best candidate for a carrier of the observed lines is the protonated cyanodiacetylene cation, HC5NH+. The HC5N/HC5NH+ ratio derived in TMC-1 is 240, which is very similar to the HC3N/HC3NH+ ratio. Results are discussed in the framework of a chemical model for protonated molecules in cold dense clouds.

Gas-phase kinetics of CH3CHO with OH radicals between 11.7 and 177.5 K

Gas-phase reactions in the interstellar medium (ISM) are a source of molecules in this environment. The knowledge of the rate coefficient for neutral-neutral reactions as a function of temperature, k(T), is essential to improve astrochemical models. In this work, we have experimentally measured k(T) for the reaction between the OH radical and acetaldehyde, both present in many sources of the ISM. Laser techniques coupled to a CRESU system were used to perform the kinetic measurements. The obtained modified Arrhenius equation is k(T = 11.7-177.5 K) = (1.2 ± 0.2) × 10-11 (T/300 K)-(1.8±0.1) exp-{(28.7 ± 2.5)/T} cm3 molecule-1 s-1. The k(T) value of the title reaction has been measured for the first time below 60 K. No pressure dependence of k(T) was observed at ca. 21, 50, 64 and 106 K. Finally, a pure gas-phase model indicates that the title reaction could become the main CH3CO formation pathway in dark molecular clouds, assuming that CH3CO is the main reaction product at 10 K.

Interstellar nitrile anions: Detection of C3N- and C5N- in TMC-1. ASTRONOMY AND ASTROPHYSICS 2020;641:L9. [PMID: 33173234 PMCID: PMC7116340 DOI: 10.1051/0004-6361/202039231]

We report on the first detection of C3N- and C5N- towards the cold dark core TMC-1 in the Taurus region, using the Yebes 40 m telescope. The observed C3N/C3N- and C5N/C5N- abundance ratios are ~140 and ~2, respectively; that is similar to those found in the circumstellar envelope of the carbon-rich star IRC +10216. Although the formation mechanisms for the neutrals are different in interstellar (ion-neutral reactions) and circumstellar clouds (photodissociation and radical-neutral reactions), the similarity of the C3N/C3N- and C5N/C5N- abundance ratios strongly suggests a common chemical path for the formation of these anions in interstellar and circumstellar clouds. We discuss the role of radiative electronic attachment, reactions between N atoms and carbon chain anions C n -, and that of H- reactions with HC3N and HC5N as possible routes to form C n N-. The detection of C5N- in TMC-1 gives strong support for assigning to this anion the lines found in IRC +10216, as it excludes the possibility of a metal-bearing species, or a vibrationally excited state. New sets of rotational parameters have been derived from the observed frequencies in TMC-1 and IRC +10216 for C5N- and the neutral radical C5N.

Synthesis and Characterization of Cyanobutadiene Isomers-Molecules of Astrochemical Significance

Four cyanobutadiene isomers of considerable interest to the organic chemistry, molecular spectroscopy, and astrochemistry communities were synthesized in good yields and isolated as pure compounds: (E)-1-cyano-1,3-butadiene (E-1), (Z)-1-cyano-1,3-butadiene (Z-1), 4-cyano-1,2-butadiene (2), and 2-cyano-1,3-butadiene (3). A diastereoselective synthesis was developed to generate (E)-1-cyano-1,3-butadiene (1) (10:1 E/Z) via tandem S_N2 and E2' reactions. The potential energy surfaces of the E2' reactions leading to (E)- and (Z)-1-cyano-1,3-butadiene (1) were analyzed by density functional theory calculations, and the observed diastereoselectivity was rationalized in the context of the Curtin-Hammett principle. The preparation of pure samples of these reactive compounds enables measurement of their laboratory rotational spectra, which are the critical data needed to search for these species in space by radioastronomy.

Gas phase Elemental abundances in Molecular cloudS (GEMS) II. On the quest for the sulphur reservoir in molecular clouds: the H2S case

CONTEXT

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.

AIMS

Our goal is to investigate the H2S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir.

METHODS

Using millimeter observations of CS, SO, H2S, and their isotopologues, we determine the physical conditions and H2S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H2S abundance.

RESULTS

Our modeling shows that chemical desorption is the main source of gas-phase H2S in dark cores. The measured H2S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when n H > 2 × 104. This change in the desorption rate is consistent with the formation of thick H2O and CO ice mantles on grain surfaces. The observed SO and H2S 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.

CONCLUSIONS

The gaseous H2S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H2S. The behavior of the observed H2S abundance suggests a changing desorption efficiency, which would probe the snowline in these cold cores. Our model, however, highly overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, we can only conclude that our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

A global CHIPR potential energy surface for ground-state C3H and exploratory dynamics studies of reaction C2 + CH → C3 + H

A full-dimensional global potential-energy surface (PES) is first reported for ground-state doublet C₃H using the combined-hyperbolic-inverse-power-representation (CHIPR) method and accurate ab initio energies extrapolated to the complete basis set limit. The PES is based on a many-body expansion-type development where the two-body and three-body energy terms are from our previously reported analytic potentials for C₂H(²A') and C₃(¹A',³A'), while the effective four-body one is calibrated using an extension of the CHIPR formalism for tetratomics. The final form is shown to accurately reproduce all known stationary structures of the PES, some of which are unreported thus far, and their interconversion pathways. Moreover, it warrants by built-in construction the appropriate permutational symmetry and describes in a physically reasonable manner all long-range features and the correct asymptotic behavior at dissociation. Exploratory quasi-classical trajectory calculations for the reaction C₂ + CH → C₃ + H are also performed, yielding thermalized rate coefficients for temperatures up to 4000 K.

Suprathermal particles in astrochemistry

The kinetic Monte Carlo method, used in astrochemistry to investigate suprathermal (hot) particles at the molecular level, is presented. Different modifications of this method, aimed at studying the influence of suprathermal particles in the processes occurring in gas and dust envelopes surrounding astrophysical objects — prestellar and protostellar cores of molecular clouds, planets, their moons, and comets in the Solar and extrasolar planetary systems — are considered. The important role of the fraction of suprathermal particles in astrochemical applications of this approach is demonstrated. The presence of these particles leads to local changes in the chemical composition; causes non-thermal emissions in gas and dust envelopes; enhances the chemical exchange between the gas and dust fractions of envelope; leads to the formation of extended hot coronae of planets; increases non-thermal atmospheric losses, thus determining the evolution of planetary atmosphere on astronomical time scales; and facilitates the formation of complex molecules in gas and dust envelopes of astrophysical objects.

The bibliography includes 146 references.

Gas-Phase Vibrational Spectroscopy of the Hydrocarbon Cations l-C3H+, HC3H+, and c-C3H2+: Structures, Isomers, and the Influence of Ne-Tagging

We report the first gas-phase vibrational spectra of the hydrocarbon ions C₃H⁺ and C₃H₂⁺. The ions were produced by electron impact ionization of allene. Vibrational spectra of the mass-selected ions tagged with Ne were recorded using infrared predissociation spectroscopy in a cryogenic ion trap instrument using the intense and widely tunable radiation of a free electron laser. Comparison of high-level quantum chemical calculations and resonant depletion measurements revealed that the C₃H⁺ ion is exclusively formed in its most stable linear isomeric form, whereas two isomers were observed for C₃H₂⁺. Bands of the energetically favored cyclic c-C₃H₂⁺ are in excellent agreement with calculated anharmonic frequencies, whereas for the linear open-shell HCCCH⁺ (²Π_g) a detailed theoretical description of the spectrum remains challenging because of Renner–Teller and spin–orbit interactions. Good agreement between theory and experiment, however, is observed for the frequencies of the stretching modes for which an anharmonic treatment was possible. In the case of linear l-C₃H⁺, small but non-negligible effects of the attached Ne on the ion fundamental band positions and the overall spectrum were found.

Abundances of sulphur molecules in the Horsehead nebula First NS+ detection in a photodissociation region

CONTEXT

Sulphur is one of the most abundant elements in the Universe (S/H∼1.3×10 -5 ) and plays a crucial role in biological systems on Earth. The understanding of its chemistry is therefore of major importance.

AIMS

Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER (Wide-band High-resolution Iram-30m Surveys at two positions with Emir Receivers) millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead.

METHODS

The Monte Carlo Markov Chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances.

RESULTS

A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor of > two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers.

CONCLUSIONS

This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of ∼ a few 10-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 10-8, only ∼0.1% of the cosmic sulphur abundance.

A sensitive λ 3 mm line survey of L483: A broad view of the chemical composition of a core around a Class 0 object

An exhaustive chemical characterization of dense cores is mandatory to our understanding of chemical composition changes from a starless to a protostellar stage. However, only a few sources have had their molecular composition characterized in detail. Here we present a λ 3 mm line survey of L483, a dense core around a Class 0 protostar, which was observed with the IRAM 30m telescope in the 80-116 GHz frequency range. We detected 71 molecules (140 including different isotopologs), most of which are present in the cold and quiescent ambient cloud according to their narrow lines (FWHM ~0.5 km s-1) and low rotational temperatures (≲10 K). Of particular interest among the detected molecules are the cis isomer of HCOOH, the complex organic molecules HCOOCH3, CH3OCH3, and C2H5OH, a wide variety of carbon chains, nitrogen oxides like N2O, and saturated molecules like CH3SH, in addition to eight new interstellar molecules (HCCO, HCS, HSC, NCCNH+, CNCN, NCO, H2NCO+, and NS+) whose detection has already been reported. In general, fractional molecular abundances in L483 are systematically lower than in TMC-1 (especially for carbon chains), tend to be higher than in L1544 and B1-b, and are similar to those in L1527. Apart from the overabundance of carbon chains in TMC-1, we find that L483 does not have a marked chemical differentiation with respect to starless/prestellar cores like TMC-1 and L1544, although it does chemically differentiate from Class 0 hot corino sources like IRAS 16293-2422. This fact suggests that the chemical composition of the ambient cloud of some Class 0 sources could be largely inherited from the dark cloud starless/prestellar phase. We explore the use of potential chemical evolutionary indicators, such as the HNCO/C3S, SO2/C2S, and CH3SH/C2S ratios, to trace the prestellar/protostellar transition. We also derived isotopic ratios for a variety of molecules, many of which show isotopic ratios close to the values for the local interstellar medium (remarkably all those involving 34S and 33S), while there are also several isotopic anomalies like an extreme depletion in 13C for one of the two isotopologs of c-C3H2, a drastic enrichment in 18O for SO and HNCO (SO being also largely enriched in 17O), and different abundances for the two 13C substituted species of C2H and the two 15N substituted species of N2H+. We report the first detection in space of some minor isotopologs like c-C3D. The exhaustive chemical characterization of L483 presented here, together with similar studies of other prestellar and protostellar sources, should allow us to identify the main factors that regulate the chemical composition of cores along the process of formation of low-mass protostars.

Gas phase Elemental abundances in Molecular cloudS (GEMS): I. The prototypical dark cloud TMC 1

GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3 mm and 2 mm) and the 40m Yebes telescope (1.3 cm and 7 mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from A V ~ 3 to ~20 mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of 1-5×103 cm-3; and ii) the dense phase, located at A V > 10 mag, with molecular hydrogen densities >104 cm-3. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone (A V ~ 3 mag) where C/H ~ 8×10-5 and C/O~0.8-1, until the beginning of the dense phase at A V ~ 10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H ~ (0.4 - 2.2) ×10-6, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H ~8×10-8). Based on our chemical modeling, we constrain the value of ζ H2 to ~ (0.5 - 1.8) ×10-16 s-1 in the translucent cloud.

The chemistry of disks around T Tauri and Herbig Ae/Be stars

CONTEXT

Infrared and (sub-)mm observations of disks around T Tauri and Herbig Ae/Be stars point to a chemical differentiation between both types of disks, with a lower detection rate of molecules in disks around hotter stars.

AIMS

To investigate the underlying causes of the chemical differentiation indicated by observations we perform a comparative study of the chemistry of T Tauri and Herbig Ae/Be disks. This is one of the first studies to compare chemistry in the outer regions of these two types of disks.

METHODS

We developed a model to compute the chemical composition of a generic protoplanetary disk, with particular attention to the photochemistry, and applied it to a T Tauri and a Herbig Ae/Be disk. We compiled cross sections and computed photodissociation and photoionization rates at each location in the disk by solving the FUV radiative transfer in a 1+1D approach using the Meudon PDR code and adopting observed stellar spectra.

RESULTS

The warmer disk temperatures and higher ultraviolet flux of Herbig stars compared to T Tauri stars induce some differences in the disk chemistry. In the hot inner regions, H2O, and simple organic molecules like C2H2, HCN, and CH4 are predicted to be very abundant in T Tauri disks and even more in Herbig Ae/Be disks, in contrast with infrared observations that find a much lower detection rate of water and simple organics toward disks around hotter stars. In the outer regions, the model indicates that the molecules typically observed in disks, like HCN, CN, C2H, H2CO, CS, SO, and HCO+, do not have drastic abundance differences between T Tauri and Herbig Ae disks. Some species produced under the action of photochemistry, like C2H and CN, are predicted to have slightly lower abundances around Herbig Ae stars due to a narrowing of the photochemically active layer. Observations indeed suggest that these radicals are somewhat less abundant in Herbig Ae disks, although in any case the inferred abundance differences are small, of a factor of a few at most. A clear chemical differentiation between both types of disks concerns ices. Owing to the warmer temperatures of Herbig Ae disks, one expects snowlines lying farther away from the star and a lower mass of ices compared to T Tauri disks.

CONCLUSIONS

The global chemical behavior of T Tauri and Herbig Ae/Be disks is quite similar. The main differences are driven by the warmer temperatures of the latter, which result in a larger reservoir or water and simple organics in the inner regions and a lower mass of ices in the outer disk.

Discovery of Interstellar Isocyanogen (CNCN): Further Evidence that Dicyanopolyynes are Abundant in Space

It is thought that dicyanopolyynes could be potentially abundant interstellar molecules, although their lack of dipole moment makes it impossible to detect them through radioastronomical techniques. Recently, the simplest member of this chemical family, cyanogen (NCCN), was indirectly probed for the first time in interstellar space through the detection of its protonated form toward the dense clouds L483 and TMC-1. Here we present a second firm evidence of the presence of NCCN in interstellar space, namely the detection of the metastable and polar isomer isocyanogen (CNCN). This species has been identified in L483 and tentatively in TMC-1 by observing various rotational transitions in the λ 3 mm band with the IRAM 30m telescope. We derive beam-averaged column densities for CNCN of 1.6 × 1012 cm-2 in L483 and 9 × 1011 cm-2 in TMC-1, which imply fractional abundances relative to H2 in the range (5 - 9) × 10-11. While the presence of NCCN in interstellar clouds seems out of doubt owing to the detection of NCCNH+ and CNCN, putting tight constraints on its abundance is still hampered by the poor knowledge of the chemistry that links NCCN with NCCNH+ and especially with CNCN. We estimate that NCCN could be fairly abundant, in the range 10-9-10-7 relative to H2, as other abundant nitriles like HCN and HC3N.

Optical Control of Reactions between Water and Laser-Cooled Be+ Ions

We investigate reactions between laser-cooled Be⁺ ions and room-temperature water molecules using an integrated ion trap and high-resolution time-of-flight mass spectrometer. This system allows simultaneous measurement of individual reaction rates that are resolved by reaction product. The rate coefficient of the Be⁺(²S_1/2) + H₂O → BeOH⁺ + H reaction is measured for the first time and is found to be approximately two times smaller than predicted by an ion-dipole capture model. Zero-point-corrected quasi-classical trajectory calculations on a highly accurate potential energy surface for the ground electronic state reveal that the reaction is capture-dominated, but a submerged barrier in the product channel lowers the reactivity. Furthermore, laser excitation of the ions from the ²S_1/2 ground state to the ²P_3/2 state opens new reaction channels, and we report the rate and branching ratio of the Be⁺(²P_3/2) + H₂O → BeOH⁺ + H and H₂O⁺ + Be reactions. The excited-state reactions are nonadiabatic in nature.

Detection of interstellar HCS and its metastable isomer HSC: new pieces in the puzzle of sulfur chemistry

We present the first identification in interstellar space of the thioformyl radical (HCS) and its metastable isomer HSC. These species were detected toward the molecular cloud L483 thanks to observations carried out with the IRAM 30m telescope in the λ 3 mm band. We derive beam-averaged column densities of 7 × 1012 cm-2 for HCS and 1.8 × 1011 cm-2 for HSC, which translate to fractional abundances relative to H2 of 2 × 10-10 and 6 × 10-12, respectively. Although the amount of sulfur locked by these radicals is low, their detection allows to put interesting constraints on the chemistry of sulfur in dark clouds. Interestingly, the H2CS/HCS abundance ratio is found to be quite low, ~ 1, in contrast with the oxygen analogue case, in which the H2CO/HCO abundance ratio is around 10 in dark clouds. Moreover, the radical HCS is found to be more abundant than its oxygen analogue, HCO. The metastable species HOC, the oxygen analogue of HSC, has not been yet observed in space. These observational constraints are confronted with the outcome of a recent model of the chemistry of sulfur in dark clouds. The model underestimates the fractional abundance of HCS by at least one order of magnitude, overestimates the H2CS/HCS abundance ratio, and does not provide an abundance prediction for the metastable isomer HSC. These observations should prompt a revision of the chemistry of sulfur in interstellar clouds.

A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

The O(¹D) + H₂ reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(¹D) + D₂ reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.

Discovery of the Ubiquitous Cation NS+ in Space Confirmed by Laboratory Spectroscopy

We report the detection in space of a new molecular species which has been characterized spectroscopically and fully identified from astrophysical data. The observations were carried out with the 30m IRAM telescopea. The molecule is ubiquitous as its J=2→1 transition has been found in cold molecular clouds, prestellar cores, and shocks. However, it is not found in the hot cores of Orion-KL and in the carbon-rich evolved star IRC+10216. Three rotational transitions in perfect harmonic relation J' = 2/3/5 have been identified in the prestellar core B1b. The molecule has a 1Σ electronic ground state and its J=2→1 transition presents the hyperfine structure characteristic of a molecule containing a nucleus with spin 1. A careful analysis of possible carriers shows that the best candidate is NS+. The derived rotational constant agrees within 0.3-0.7% with ab initio calculations. NS+ was also produced in the laboratory to unambiguously validate the astrophysical assignment. The observed rotational frequencies and determined molecular constants confirm the discovery of the nitrogen sulfide cation in space. The chemistry of NS+ and related nitrogen-bearing species has been analyzed by means of a time-dependent gas phase model. The model reproduces well the observed NS/NS+ abundance ratio, in the range 30-50, and indicates that NS+ is formed by reactions of the neutral atoms N and S with the cations SH+ and NH+, respectively.

Is the gas-phase OH+H2CO reaction a source of HCO in interstellar cold dark clouds? A kinetic, dynamic and modelling study

Chemical kinetics of neutral-neutral gas-phase reactions at ultralow temperatures is a fascinating research subject with important implications on the chemistry of complex organic molecules in the interstellar medium (T∼10-100K). Scarce kinetic information is currently available for this kind of reactions at T<200 K. In this work we use the CRESU (Cinétique de Réaction en Ecoulement Supersonique Uniforme, which means Reaction Kinetics in a Uniform Supersonic Flow) technique to measure for the first time the rate coefficients (k) of the gas-phase OH+H2CO reaction between 22 and 107 K. k values greatly increase from 2.1×10-11 cm3 s-1 at 107 K to 1.2×10-10 cm3 s-1 at 22 K. This is also confirmed by quasi-classical trajectories (QCT) at collision energies down to 0.1 meV performed using a new full dimension and ab initio potential energy surface, recently developed which generates highly accurate potential and includes long range dipole-dipole interactions. QCT calculations indicate that at low temperatures HCO is the exclusive product for the OH+H2CO reaction. In order to revisit the chemistry of HCO in cold dense clouds, k is reasonably extrapolated from the experimental results at 10K (2.6×10-10 cm3 s-1). The modeled abundances of HCO are in agreement with the observations in cold dark clouds for an evolving time of 105-106 yrs. The different sources of production of HCO are presented and the uncertainties in the chemical networks discussed. This reaction can be expected to be a competitive process in the chemistry of prestellar cores. The present reaction is shown to account for a few percent of the total HCO production rate. Extensions to photodissociation regions and diffuse clouds environments are also commented.