-
Quantum study of the CH$_3^+$ photodissociation in full dimension Neural Networks potential energy surfaces
Authors:
Pablo del Mazo-Sevillano,
Alfredo Aguado,
Javier R. Goicoechea,
Octavio Roncero
Abstract:
CH$_3^+$, a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction {\it ab initio} calculations are done, and full dimensional potential energy surfaces are developed for the three lower electronic…
▽ More
CH$_3^+$, a cornerstone intermediate in interstellar chemistry, has recently been detected for the first time by the James Webb Space Telescope. The photodissociation of this ion is studied here. Accurate explicitly correlated multi-reference configuration interaction {\it ab initio} calculations are done, and full dimensional potential energy surfaces are developed for the three lower electronic states, with a fundamental invariant neural network method. The photodissociation cross section is calculated using a full dimensional quantum wave packet method, in heliocentric Radau coordinates. The wave packet is represented in angular and radial grids allowing to reduce the number of points physically accessible, requiring to push up the spurious states appearing when evaluating the angular kinetic terms, through a projection technique. The photodissociation spectra, when employed in astrochemical models to simulate the conditions of the Orion Bar, results in a lesser destruction of CH$_3^+$ compared to that obtained when utilizing the recommended values in the kinetic database for astrochemistry (KIDA).
△ Less
Submitted 23 April, 2024;
originally announced April 2024.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) VIII. Unlocking the CS chemistry: the CH + S$\rightarrow$ CS + H and C$_2$ + S$\rightarrow$ CS + C reactions
Authors:
Carlos M. R. Rocha,
Octavio Roncero,
Niyazi Bulut,
Piotr Zuchowski,
David Navarro-Almaida,
Asuncion Fuente,
Valentine Wakelam,
Jean-Christophe Loison,
Evelyne Roueff,
Javier R. Goicoechea,
Gisela Esplugues,
Leire Beitia-Antero,
Paola Caselli,
Valerio Lattanzi,
Jaime Pineda,
Romane Le Gal,
Marina Rodriguez-Baras,
Pablo Riviere-Marichalar
Abstract:
We revise the rates of reactions CH + S -> CS + H and C_2 + S -> CS + C, important CS formation routes in dark and diffuse warm gas. We performed ab initio calculations to characterize the main features of all the electronic states correlating to the open shell reactants. For CH+S we have calculated the full potential energy surfaces for the lowest doublet states and the reaction rate constant wit…
▽ More
We revise the rates of reactions CH + S -> CS + H and C_2 + S -> CS + C, important CS formation routes in dark and diffuse warm gas. We performed ab initio calculations to characterize the main features of all the electronic states correlating to the open shell reactants. For CH+S we have calculated the full potential energy surfaces for the lowest doublet states and the reaction rate constant with a quasi-classical method. For C_2+S, the reaction can only take place through the three lower triplet states, which all present deep insertion wells. A detailed study of the long-range interactions for these triplet states allowed to apply a statistic adiabatic method to determine the rate constants. This study of the CH + S reaction shows that its rate is nearly independent on the temperature in a range of 10-500 K with an almost constant value of 5.5 10^{-11} cm^3/s at temperatures above 100~K. This is a factor \sim 2-3 lower than the value obtained with the capture model. The rate of the reaction C_2 + S depends on the temperature taking values close to 2.0 10^{-10} cm^3/s at low temperatures and increasing to 5. 10^{-10} cm^3/s for temperatures higher than 200~K. Our modeling provides a rate higher than the one currently used by factor of \sim 2. These reactions were selected for involving open-shell species with many degenerate electronic states, and the results obtained in the present detailed calculations provide values which differ a factor of \sim 2-3 from the simpler classical capture method. We have updated the sulphur network with these new rates and compare our results in the prototypical case of TMC1 (CP). We find a reasonable agreement between model predictions and observations with a sulphur depletion factor of 20 relative to the sulphur cosmic abundance, but it is not possible to fit all sulphur-bearing molecules better than a factor of 10 at the same chemical time.
△ Less
Submitted 1 July, 2023;
originally announced July 2023.
-
The magnesium paradigm in IRC+10216: Discovery of MgC$_4$H$^+$, MgC$_3$N$^+$, MgC$_6$H$^+$, and MgC$_5$N$^+$
Authors:
J. Cernicharo,
C. Cabezas,
J. R. Pardo,
M. Agúndez,
O. Roncero,
B. Tercero,
N. marcelino,
M. Guélin,
Y. Endo,
P. deVicente
Abstract:
We found four series of harmonically related lines in IRC\,+10216 with the Yebes\,40m and IRAM\,30m telescopes. The first series corresponds to a molecule with a rotational constant, $B$, of 1448.5994$\pm$0.0013 MHz and a distortion constant, $D$, of 63.45$\pm$1.15 Hz and covers upper quantum numbers from $J_u$=11 up to 33 (B1449). The second series is fitted with $B$=1446.9380$\pm$0.0098 MHz and…
▽ More
We found four series of harmonically related lines in IRC\,+10216 with the Yebes\,40m and IRAM\,30m telescopes. The first series corresponds to a molecule with a rotational constant, $B$, of 1448.5994$\pm$0.0013 MHz and a distortion constant, $D$, of 63.45$\pm$1.15 Hz and covers upper quantum numbers from $J_u$=11 up to 33 (B1449). The second series is fitted with $B$=1446.9380$\pm$0.0098 MHz and $D$=91$\pm$23 Hz and covers upper quantum numbers from $J_u$=11 up to 17 (B1447). The third series is fitted with $B$=598.7495$\pm$0.0011 MHz and D=6.13$\pm$0.43 Hz and covers quantum numbers from $J_u$=26 up to 41 (B599). Finally, the frequencies of the last series of lines can be reproduced with $B$=594.3176$\pm$0.0026 MHz and $D$=4.92$\pm$1.16 Hz (B594). The large values of $D$ point toward four metal-bearing carriers. After exploring all plausible candidates containing Na, Al, Mg, and other metals, our ab initio calculations indicate that the cations MgC$_4$H$^+$, MgC$_3$N$^+$, MgC$_6$H$^+$, and MgC$_5$N$^+$ must be the carriers of B1449, B1447, B599, and B594, respectively. These cations could be formed by the radiative association of Mg$^+$ with C$_4$H, C$_3$N, C$_6$H, and C$_5$N, respectively. We calculated the radiative association rate coefficient of Mg$^+$ with C$_4$H, C$_3$N, C$_6$H, and C$_5$N and incorporated them in our chemical model. The results confirm that the Mg-bearing cations can be formed through these radiative association reactions in the outer layers of IRC\,+10216. This is the first time that cationic metal-bearing species have been found in space. These results provide a new paradigm on the reactivity of ionized metals with abundant radicals and open the door for further characterization of similar species in metal-rich astrophysical environments.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
The chemistry of H2NC in the interstellar medium and the role of the C + NH3 reaction
Authors:
M. Agundez,
O. Roncero,
N. Marcelino,
C. Cabezas,
B. Tercero,
J. Cernicharo
Abstract:
We carried out an observational search for the recently discovered molecule H2NC, and its more stable isomer H2CN, toward eight cold dense clouds (L1544, L134N, TMC-2, Lupus-1A, L1489, TMC-1 NH3, L1498, and L1641N) and two diffuse clouds (B0415+379 and B0355+508) in an attempt to constrain its abundance in different types of interstellar regions and shed light on its formation mechanism. We detect…
▽ More
We carried out an observational search for the recently discovered molecule H2NC, and its more stable isomer H2CN, toward eight cold dense clouds (L1544, L134N, TMC-2, Lupus-1A, L1489, TMC-1 NH3, L1498, and L1641N) and two diffuse clouds (B0415+379 and B0355+508) in an attempt to constrain its abundance in different types of interstellar regions and shed light on its formation mechanism. We detected H2NC in most of the cold dense clouds targeted, 7 out of 8, while H2CN was only detected in 5 out of 8 clouds. The column densities derived for both H2NC and H2CN are in the range 1e11-1e12 cm-2 and the abundance ratio H2NC/H2CN varies between 0.51 and >2.7. The metastable isomer H2NC is therefore widespread in cold dense clouds where it is present with an abundance similar to that of H2CN. We did not detect either H2NC or H2CN in any of the two diffuse clouds targeted, which does not allow to shed light on how the chemistry of H2NC and H2CN varies between dense and diffuse clouds. We found that the column density of H2NC is correlated with that of NH3, which strongly suggests that these two molecules are chemically linked, most likely ammonia being a precursor of H2NC through the C + NH3 reaction. We performed electronic structure and statistical calculations which show that both H2CN and H2NC can be formed in the C + NH3 reaction through two different channels involving two different transition states which lie very close in energy. The predicted product branching ratio H2NC/H2CN is very method dependent but values between 0.5 and 0.8 are the most likely ones. Therefore, both the astronomical observations and the theoretical calculations support that the reaction C + NH3 is the main source of H2NC in interstellar clouds.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
Vibrational, non-adiabatic and isotopic effects in the dynamics of the H$_2$ + H$_2^+$ $\rightarrow$ H$_3^+$ + H reaction: application to plasma modeling
Authors:
P. del Mazo-Sevillano,
D. Félix-González,
A. Aguado,
C. Sanz-Sanz,
D. -H. Kwon,
O. Roncero
Abstract:
The title reaction is studied using a quasi-classical trajectory method for collision energies between 0.1 meV and 10 eV, considering the vibrational excitation of H$_2^+$ reactant. A new potential energy surface is developed based on a Neural Network many body correction of a triatomics-in-molecules potential, which significantly improves the accuracy of the potential up to energies of 17 eV, hig…
▽ More
The title reaction is studied using a quasi-classical trajectory method for collision energies between 0.1 meV and 10 eV, considering the vibrational excitation of H$_2^+$ reactant. A new potential energy surface is developed based on a Neural Network many body correction of a triatomics-in-molecules potential, which significantly improves the accuracy of the potential up to energies of 17 eV, higher than in other previous fits.The effect of the fit accuracy and the non-adiabatic transitions on the dynamics are analyzed in detail.The reaction cross section for collision energies above 1 eV increases significantly with the increasing of the vibrational excitation of H$_2^+$($v'$), for values up to $v'$=6. The total reaction cross section (including the double fragmentation channel) obtained for $v'$=6 matches the new experimental results obtained by Savic, Schlemmer and Gerlich [Chem. Phys. Chem. 21 (13), 1429.1435(2020)]. The differences among several experimental setups, for collision energies above 1 eV, showing cross sections scattered/dispersed over a rather wide interval, can be explained by the differences in the vibrational excitations obtained in the formation of H$_2^+$ reactants. On the contrary, for collision energies below 1 eV, the cross section is determined by the long range behavior of the potential and do not depend strongly on the vibrational state of H$_2^+$. In addition in this study, the calculated reaction cross sections are used in a plasma model and compared with previous results. We conclude that the efficiency of the formation of H$_3^+$ in the plasma is affected by the potential energy surface used.
△ Less
Submitted 3 March, 2023;
originally announced March 2023.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) VII. Sulfur elemental abundance
Authors:
A. Fuente,
P. Rivière-Marichalar,
L. Beitia-Antero,
P. Caselli,
V. Wakelam,
G. Esplugues,
M. Rodríguez-Baras,
D. Navarro-Almaida,
M. Gerin,
C. Kramer,
R. Bachiller,
J. R. Goicoechea,
I. Jiménez-Serra,
J. C. Loison,
A. Ivlev,
R. Martín-Doménech,
S. Spezzano,
O. Roncero,
G. Muñoz-Caro,
S. Cazaux,
N. Marcelino
Abstract:
Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals o…
▽ More
Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals of this program. We have carried out an extensive chemical modeling of the fractional abundances of CO, HCO$^+$, HCN, HNC, CS, SO, H$_2$S, OCS, and HCS$^+$ to determine the sulfur depletion toward the 244 positions in the GEMS database. These positions sample visual extinctions from A$_V$ $\sim$ 3 mag to $>$50 mag, molecular hydrogen densities ranging from a few 10$^3$~cm$^{-3}$ to 3$\times$10$^6$~cm$^{-3}$, and T$_k$ $\sim$ 10$-$35 K. Most of the positions in Taurus and Perseus are best fitted assuming early-time chemistry, t=0.1 Myr, $ζ_{H_2}$$\sim$ (0.5$-$1)$\times$10$^{-16}$ s$^{-1}$, and [S/H]$\sim$1.5$\times$10$^{-6}$. On the contrary, most of the positions in Orion are fitted with t=1~Myr and $ζ_{H_2}$$\sim$ 10$^{-17}$ s$^{-1}$. Moreover, $\sim$40% of the positions in Orion are best fitted assuming the undepleted sulfur abundance, [S/H]$\sim$1.5$\times$10$^{-5}$. Our results suggest that sulfur depletion depends on the environment. While the abundances of sulfur-bearing species are consistent with undepleted sulfur in Orion, a depletion factor of $\sim$20 is required to explain those observed in Taurus and Perseus. We propose that differences in the grain charge distribution in the envelopes of the studied clouds might explain these variations. The shocks associated with past and ongoing star formation could also contribute to enhance [S/H] in Orion.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Spin-orbit transitions in the N$^+$($^3P_{J_A}$) + H$_2$ $\rightarrow$ NH$^+$($X^2Π$, $^4Σ^-$)+ H($^2S$) reaction, using adiabatic and mixed quantum-adiabatic statistical approaches
Authors:
Susana Gomez-Carrasco,
Daniel Felix-Gonzalez,
Alfredo Aguado,
Octavio Roncero
Abstract:
The cross section and rate constants for the title reaction are calculated for all the spin-orbit states of N$^+$($^3P_{J_A}$) using two statistical approaches, one purely adiabatic and the other one mixing quantum capture for the entrance channel and adiabatic treatment for the products channel. This is made by using a symmetry adapted basis set combining electronic (spin and orbital) and nuclear…
▽ More
The cross section and rate constants for the title reaction are calculated for all the spin-orbit states of N$^+$($^3P_{J_A}$) using two statistical approaches, one purely adiabatic and the other one mixing quantum capture for the entrance channel and adiabatic treatment for the products channel. This is made by using a symmetry adapted basis set combining electronic (spin and orbital) and nuclear angular momenta in the reactants channel. To this aim, accurate {\it ab initio} calculations are performed separately for reactants and products. In the reactants channel, the three lowest electronic states (without spin-orbit couplings) have been diabatized, and the spin-orbit couplings have been introduced through a model localizing the spin-orbit interactions in the N$^+$ atom, which yields accurate results as compared to {\it ab initio} calculations including spin-orbit couplings. For the products, eleven purely adiabatic spin-orbit states have been determined with {\it ab initio} calculations. The reactive rate constants thus obtained are in very good agreement with the available experimental data for several ortho-H$_2$ fractions, assuming a thermal initial distribution of spin-orbit states. The rate constants for selected spin-orbit $J_A$ states are obtained, to provide a proper validation of the spin-orbit effects to obtain the experimental rate constants.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
The role of highly vibrationally excited H2 initiating the N chemistry: Quantum study and 3-sigma detection of NH emission in the Orion Bar PDR
Authors:
Javier R. Goicoechea,
Octavio Roncero
Abstract:
The formation of hydrides by gas-phase reactions between H2 and a heavy element atom is a very selective process. Reactions with ground-state neutral carbon, oxygen, nitrogen, and sulfur atoms are very endoergic and have high energy barriers because the H2 molecule has to be fragmented before a hydride bond is formed. In cold interstellar clouds, these barriers exclude the formation of CH, OH, NH,…
▽ More
The formation of hydrides by gas-phase reactions between H2 and a heavy element atom is a very selective process. Reactions with ground-state neutral carbon, oxygen, nitrogen, and sulfur atoms are very endoergic and have high energy barriers because the H2 molecule has to be fragmented before a hydride bond is formed. In cold interstellar clouds, these barriers exclude the formation of CH, OH, NH, and SH radicals through hydrogen abstraction reactions. Here we study a very energetically unfavorable process, the reaction of N(4S) atoms with H2 molecules. We calculated the reaction rate coefficient for H2 in different vibrational levels, using quantum methods for v=0-7 and quasi-classical methods up to v=12. Owing to the high energy barrier, these rate coefficients increase with v and also with the gas temperature. We implemented the new rates in the Meudon PDR code and studied their effect on models with different ultraviolet (UV) illumination conditions. In strongly UV-irradiated dense gas (Orion Bar conditions), the presence of H2 in highly vibrationally excited levels (v>7) enhances the NH abundance by two orders of magnitude (at the PDR surface) compared to models that use the thermal rate coefficient for reaction N(4S) + H2 -> NH + H. The increase in NH column density across the PDR is a factor of ~25. We explore existing Herschel/HIFI observations of the Orion Bar and Horsehead PDRs. We report a 3-sigma emission feature at the ~974 GHz frequency of the NH N_J=1_2-0_1 line toward the Bar. The emission level implies N(NH)~10^13 cm^-2, which is consistent with PDR models using the new rate coefficients for reactions between N and UV-pumped H2. This formation route dominates over hydrogenation reactions involving the less abundant N+ ion. JWST observations will soon quantify the amount and reactivity of UV-pumped H2 in many interstellar and circumstellar environments.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
Formation of phosphorus monoxide (PO) in the interstellar medium: insights from quantum-chemical and kinetic calculations
Authors:
Juan García de la Concepción,
Cristina Puzzarini,
Vincenzo Barone,
Izaskun Jiménez-Serra,
Octavio Roncero
Abstract:
In recent years, phosphorus monoxide (PO) -- an important molecule for prebiotic chemistry -- has been detected in star-forming regions and in the comet 67P/Churyumov-Gerasimenko. These studies have revealed that, in the interstellar medium, PO is systematically the most abundant P-bearing species, with abundances that are $\sim$1-3 times greater than those derived for phosphorus nitride (PN), the…
▽ More
In recent years, phosphorus monoxide (PO) -- an important molecule for prebiotic chemistry -- has been detected in star-forming regions and in the comet 67P/Churyumov-Gerasimenko. These studies have revealed that, in the interstellar medium, PO is systematically the most abundant P-bearing species, with abundances that are $\sim$1-3 times greater than those derived for phosphorus nitride (PN), the second most abundant P-containing molecule. The reason why PO is more abundant than PN remains still unclear. Experimental studies with phosphorus in the gas phase are not available, probably because of the difficulties in dealing with its compounds. Therefore, the reactivity of atomic phosphorus needs to be investigated using reliable computational tools. To this end, state-of-the-art quantum-chemical computations have been employed to evaluate accurate reaction rates and branching ratios for the P + OH $\rightarrow$ PO + H and P + H$_2$O $\rightarrow$ PO + H$_2$ reactions in the framework of a master equation approach based on ab-initio transition state theory. The hypothesis that OH and H${_2}$O can be potential oxidizing agents of atomic phosphorus is based on the ubiquitous presence of H${_2}$O in the ISM. Its destruction then produces OH, which is another very abundant species. While the reaction of atomic phosphorus in its gound state with water is not a relevant source of PO because of emerged energy barriers, the P + OH reaction represents an important formation route of PO in the interstellar medium. Our kinetic results show that this reaction follow an Arrhenius behavior, and thus its rate coefficients alpha=2.28$\times$10$^{-10}$ cm${^3}$ molecule$^{-1}$ s$^{-1}$, beta=0.16 and gamma=0.37 K increase by increasing the temperature.
△ Less
Submitted 19 August, 2021;
originally announced August 2021.
-
Vibrational effects in the quantum dynamics of the H + D_2^+ charge transfer reaction
Authors:
O. Roncero,
V. Andrianarijaona,
A. Aguado,
C. Sanz-Sanz
Abstract:
The H + D_2^+(v=0,1 and 2) charge transfer reaction is studied using an accurate wave packet method, using recently proposed coupled diabatic potential energy surfaces. The state-to-state cross section is obtained for three different channels: non-reactive charge transfer, reactive charge transfer, and exchange reaction. The three processes proceed via the electronic transition from the first exci…
▽ More
The H + D_2^+(v=0,1 and 2) charge transfer reaction is studied using an accurate wave packet method, using recently proposed coupled diabatic potential energy surfaces. The state-to-state cross section is obtained for three different channels: non-reactive charge transfer, reactive charge transfer, and exchange reaction. The three processes proceed via the electronic transition from the first excited to the ground electronic state. The cross section for the three processes increases with the initial vibrational excitation. The non-reactive charge transfer process is the dominant channel, whose branching ratio increases with collision energy, and it compares well with experimental measurements at collision energies around 0.5 eV. For lower energies the experimental cross section is considerably higher, suggesting that it corresponds to higher vibrational excitation of D_2^+(v) reactants. Further experimental studies of this reaction and isotopic variants are needed, where conditions are controlled to obtain a better analysis of the vibrational effects of the D_2^+ reagents.
△ Less
Submitted 24 June, 2021;
originally announced June 2021.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS). IV. Observational results and statistical trends
Authors:
M. Rodríguez-Baras,
A. Fuente,
P. Riviére-Marichalar,
D. Navarro-Almaida,
P. Caselli,
M. Gerin,
C. Kramer,
E. Roueff,
V. Wakelam,
G. Esplugues,
S. García-Burillo,
R. Le Gal,
S. Spezzano,
T. Alonso-Albi,
R. Bachiller,
S. Cazaux,
B. Commercon,
J. R. Goicoechea,
J. C. Loison,
S. P. Treviño-Morales,
O. Roncero,
I. Jiménez-Serra,
J. Laas,
A. Hacar,
J. Kirk
, et al. (11 additional authors not shown)
Abstract:
Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S…
▽ More
Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies.
△ Less
Submitted 25 February, 2021;
originally announced February 2021.
-
Near-resonant effects in the quantum dynamics of the H+H$_2^+$ $\rightarrow$ H$_2$+ H$^+$ charge transfer reaction and isotopic variants
Authors:
Cristina Sanz-Sanz,
Alfredo Aguado,
Octavio Roncero
Abstract:
The non-adiabatic quantum dynamics of the H+H$_2^+$ $\rightarrow$ H$_2$+ H$^+$ charge transfer reactions, and some isotopic variants, is studied with an accurate wave packet method. A recently developed $3\times$3 diabatic potential model is used, which is based on very accurate {\it ab initio} calculations and includes the long-range interactions for ground and excited states. It is found that fo…
▽ More
The non-adiabatic quantum dynamics of the H+H$_2^+$ $\rightarrow$ H$_2$+ H$^+$ charge transfer reactions, and some isotopic variants, is studied with an accurate wave packet method. A recently developed $3\times$3 diabatic potential model is used, which is based on very accurate {\it ab initio} calculations and includes the long-range interactions for ground and excited states. It is found that for initial H$_2^+$(v=0), the quasi-degenerate H$_2$(v'=4) non-reactive charge transfer product is enhanced, producing an increase of the reaction probability and cross section. It becomes the dominant channel from collision energies above 0.2 eV, producing a ratio, between v'=4 and the rest of v's, that increases up to 1 eV. H+H$_2^+$ $\rightarrow$ H$_2^+$+ H exchange reaction channel is nearly negligible, while the reactive and non-reactive charge transfer reaction channels are of the same order, except that corresponding to H$_2$(v'=4), and the two charge transfer processes compete below 0.2 eV. This enhancement is expected to play an important vibrational and isotopic effect that need to be evaluated. For the three proton case, the problem of the permutation symmetry is discussed when using reactant Jacobi coordinates.
△ Less
Submitted 22 February, 2021;
originally announced February 2021.
-
Neural Network Potential Energy Surface for the low temperature Ring Polymer Molecular Dynamics of the H2CO + OH reaction
Authors:
Pablo del Mazo-Sevillano,
Alfredo Aguado,
Octavio Roncero
Abstract:
A new potential energy surface (PES) and dynamical study are presented of the reactive process between H2CO + OH towards the formation of HCO + H2O and HCOOH + H. In this work a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, what prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PE…
▽ More
A new potential energy surface (PES) and dynamical study are presented of the reactive process between H2CO + OH towards the formation of HCO + H2O and HCOOH + H. In this work a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, what prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many body term has been used fitted with a novel artificial neural networks scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing a good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated to excitations of the normal modes of the ring polymer.
△ Less
Submitted 4 February, 2021;
originally announced February 2021.
-
Bottlenecks to interstellar sulfur chemistry: Sulfur-bearing hydrides in UV-illuminated gas and grains
Authors:
J. R. Goicoechea,
A. Aguado,
S. Cuadrado,
O. Roncero,
J. Pety,
E. Bron,
A. Fuente,
D. Riquelme,
E. Chapillon,
C. Herrera,
C. A. Duran
Abstract:
Hydride molecules lie at the base of interstellar chemistry, but the synthesis of sulfuretted hydrides is poorly understood. Motivated by new observations of the Orion Bar PDR - 1'' resolution ALMA images of SH+; IRAM 30m detections of H2S, H2S34, and H2S33; H3S+ (upper limits); and SOFIA observations of SH - we perform a systematic study of the chemistry of S-bearing hydrides. We determine their…
▽ More
Hydride molecules lie at the base of interstellar chemistry, but the synthesis of sulfuretted hydrides is poorly understood. Motivated by new observations of the Orion Bar PDR - 1'' resolution ALMA images of SH+; IRAM 30m detections of H2S, H2S34, and H2S33; H3S+ (upper limits); and SOFIA observations of SH - we perform a systematic study of the chemistry of S-bearing hydrides. We determine their column densities using coupled excitation, radiative transfer as well as chemical formation and destruction models. We revise some of the key gas-phase reactions that lead to their chemical synthesis. This includes ab initio quantum calculations of the vibrational-state-dependent reactions SH+ + H2 <-> H2S+ + H and S + H2 <-> SH + H. We find that reactions of UV-pumped H2 (v>1) with S+ explain the presence of SH+ in a high thermal-pressure gas component, P_th~10^8 cm^-3 K, close to the H2 dissociation front. However, subsequent hydrogen abstraction reactions of SH+, H2S+, and S with vibrationally excited H2, fail to ultimately explain the observed H2S column density (~2.5x10^14 cm^-2, with an ortho-to-para ratio of 2.9+/-0.3). To overcome these bottlenecks, we build PDR models that include a simple network of grain surface reactions leading to the formation of solid H2S (s-H2S). The higher adsorption binding energies of S and SH suggested by recent studies imply that S atoms adsorb on grains (and form s-H2S) at warmer dust temperatures and closer to the UV-illuminated edges of molecular clouds. Photodesorption and, to a lesser extent, chemical desorption, produce roughly the same H2S column density (a few 10^14 cm-^2) and abundance peak (a few 10^-8) nearly independently of n_H and G_0. This agrees with the observed H2S column density in the Orion Bar as well as at the edges of dark clouds without invoking substantial depletion of elemental sulfur abundances.
△ Less
Submitted 11 February, 2021; v1 submitted 4 January, 2021;
originally announced January 2021.
-
Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction
Authors:
N. Bulut,
O. Roncero,
A. Aguado,
J. -C. Loison,
D. Navarro-Almaida,
V. Wakelam,
A. Fuente,
E. Roueff,
R. Le Gal,
P. Caselli,
M. Gerinm K. M. Hickson,
S. Spezzano,
P. Riviere-Marichalar,
T. Alonso-Albi,
R. Bachiller,
I. Jimenez-Serra,
C. Kramer,
B. Tercero,
M. Rodriguez-Baras,
S. Garcia-Burillo,
J. R. Goicoechea,
S. P. Treviño-Morales,
G. Esplugues,
S. Cazaux,
B. Commercon
, et al. (14 additional authors not shown)
Abstract:
CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the el…
▽ More
CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance.
△ Less
Submitted 18 December, 2020;
originally announced December 2020.
-
Hyperfine excitation of SH$^+$ by H
Authors:
François Lique,
Alexandre Zanchet,
Niyazi Bulut,
Javier R. Goicoechea,
Octavio Roncero
Abstract:
SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of…
▽ More
SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of the physical conditions prevailing in these energetic environments relies on knowing the rate coefficients of inelastic collisions between SH$^+$ molecules and hydrogen atoms, hydrogen molecules, and electrons. In this paper, we derive SH$^+$--H fine and hyperfine-resolved rate coefficients from the recent quantum calculations for the SH$^+$--H collisions, including inelastic, exchange and reactive processes. The method used is based on the infinite order sudden approach. State-to-state rate coefficients between the first 31 fine levels and 61 hyperfine levels of SH$^+$ were obtained for temperatures ranging from 10 to 1000 K. Fine structure-resolved rate coefficients present a strong propensity rule in favour of $Δj = ΔN$ transitions. The $Δj = ΔF$ propensity rule is observed for the hyperfine transitions. {The new rate coefficients will help significantly in the interpretation of SH$^+$ spectra from PDRs and UV-irradiated shocks where the abundance of hydrogen atoms with respect to hydrogen molecules can be significant.
△ Less
Submitted 26 April, 2020;
originally announced April 2020.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS). II. On the quest for the sulphur reservoir in molecular clouds: the $H_{2}S$ case
Authors:
D. Navarro-Almaida,
R. Le Gal,
A. Fuente,
P. Rivière-Marichalar,
V. Wakelam,
S. Cazaux,
P. Caselli,
Jacob C. Laas,
T. Alonso-Albi,
J. C. Loison,
M. Gerin,
C. Kramer,
E. Roueff,
R. Bachiller,
B. Commerçon,
R. Friesen,
S. García-Burillo,
J. R. Goicoechea,
B. M. Giuliano,
I. Jiménez-Serra,
J. M. Kirk,
V. Lattanzi,
J. Malinen,
N. Marcelino,
R. Martín-Domènech
, et al. (8 additional authors not shown)
Abstract:
Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H…
▽ More
Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.
△ Less
Submitted 7 April, 2020;
originally announced April 2020.
-
Formation of interstellar SH$^+$ from vibrationally excited H$_2$: Quantum study of S$^+$ + H$_2$ $\rightleftarrows$ SH$^+$ + H reactions and inelastic collisions
Authors:
Alexandre Zanchet,
Francois Lique,
Octavio Roncero,
Javier R. Goicoechea,
Niyazi Bulut
Abstract:
The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3Σ^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H…
▽ More
The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3Σ^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H$_2$ reaction is endothermic by $\approx$ 9860 K, and requires at least two vibrational quanta on the H$_2$ molecule to yield significant reactivity. Quasi-classical calculations of the total formation rate constant for H$_2$($v=2$) are in very good agreement with the quantum results above 100K. Further quasi-classical calculations are then performed for $v=3$, 4, and 5 to cover all vibrationally excited H$_2$ levels significantly populated in dense photodissociation regions (PDR). The new calculated formation and destruction rate constants are two to six times larger than the previous ones and have been introduced in the Meudon PDR code to simulate the physical and illuminating conditions in the Orion bar prototypical PDR. New astrochemical models based on the new molecular data produce four times larger SH$^+$ column densities, in agreement with those inferred from recent ALMA observations of the Orion bar.
△ Less
Submitted 7 May, 2019;
originally announced May 2019.
-
Gas phase Elemental abundances in Molecular cloudS (GEMS) I. The prototypical dark cloud TMC 1
Authors:
A. Fuente,
D. Navarro,
P. Caselli,
M. Gerin,
C. Kramer,
E. Roueff,
T. Alonso-Albi,
R. Bachiller,
S. Cazaux,
B. Commercon,
R. Friesen,
S. García-Burillo,
B. M. Giuliano,
J. R. Goicoechea,
P. Gratier,
A. Hacar,
I. Jiménez-Serra,
J. Kirk,
V. Lattanzi,
J. C. Loison,
P. J. Malinen,
N. Marcelino,
R. Martín-Doménech,
G. Muñoz-Caro,
J. Pineda
, et al. (8 additional authors not shown)
Abstract:
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 (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the frac…
▽ More
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 (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) 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 Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~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) x10$^{-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~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $ζ_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud.
△ Less
Submitted 7 March, 2019; v1 submitted 13 September, 2018;
originally announced September 2018.
-
Formation and destruction of SiS in space
Authors:
Alexandre Zanchet,
Octavio Roncero,
Marcelino Agúndez,
José Cernicharo
Abstract:
The presence of SiS in space seems to be restricted to a few selected types of astronomical environments.
It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr\,B2, Orion\,KL and more recently L1157-B1.
The kinetics of reactions involving SiS is very poorly known and…
▽ More
The presence of SiS in space seems to be restricted to a few selected types of astronomical environments.
It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr\,B2, Orion\,KL and more recently L1157-B1.
The kinetics of reactions involving SiS is very poorly known and here we revisit the chemistry of SiS in space by studying some potentially important reactions of formation and destruction of this molecule. We calculated {\it ab initio} potential energy surfaces of the SiOS system and computed rate coefficients in the temperature range 50-2500 K for the reaction of destruction of SiS, in collisions with atomic O, and of its formation, through the reaction between Si and SO. We find that both reactions are rapid, with rate coefficients of a few times 10$^{-10}$ cm$^3$ s$^{-1}$, almost independent of temperature.
In the reaction between Si and SO, SiO production is 5-7 times more efficient than SiS formation.
The reaction of SiS with O atoms can play an important role in destroying SiS in envelopes around evolved stars. We built a simple chemical model of a postshock gas to study the chemistry of SiS in protostellar outflows and we found that SiS forms with a lower abundance and later than SiO, that SiS is efficiently destroyed through reaction with O, and that the main SiS-forming reactions are Si + SO and Si + SO$_2$.
△ Less
Submitted 4 July, 2018;
originally announced July 2018.
-
State-to-state chemistry and rotational excitation of CH$^+$ in photon-dominated regions
Authors:
A. Faure,
P. Halvick,
T. Stoecklin,
P. Honvault,
M. D. Epée Epée,
J. Zs. Mezei,
O. Motapon,
I. F. Schneider,
J. Tennyson,
O. Roncero,
N. Bulut,
A. Zanchet
Abstract:
We present a detailed theoretical study of the rotational excitation of CH$^+$ due to reactive and nonreactive collisions involving C$^+(^2P)$, H$_2$, CH$^+$, H and free electrons. Specifically, the formation of CH$^+$ proceeds through the reaction between C$^+(^2P)$ and H$_2(ν_{\rm H_2}=1, 2)$, while the collisional (de)excitation and destruction of CH$^+$ is due to collisions with hydrogen atoms…
▽ More
We present a detailed theoretical study of the rotational excitation of CH$^+$ due to reactive and nonreactive collisions involving C$^+(^2P)$, H$_2$, CH$^+$, H and free electrons. Specifically, the formation of CH$^+$ proceeds through the reaction between C$^+(^2P)$ and H$_2(ν_{\rm H_2}=1, 2)$, while the collisional (de)excitation and destruction of CH$^+$ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000~K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH$^+$ with H atoms at kinetic temperatures below 50~K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H$_2$ has a substantial effect on the excitation of CH$^+$ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH$^+$ toward the Orion Bar and the planetary nebula NGC~7027. Our results further suggest that the population of $ν_{\rm H_2}=2$ might be significant in the photon-dominated region of NGC~7027.
△ Less
Submitted 7 April, 2017;
originally announced April 2017.
-
Trans-cis molecular photoswitching in interstellar Space
Authors:
S. Cuadrado,
J. R. Goicoechea,
O. Roncero,
A. Aguado,
B. Tercero,
J. Cernicharo
Abstract:
As many organic molecules, formic acid (HCOOH) has two conformers (trans and cis). The energy barrier to internal conversion from trans to cis is much higher than the thermal energy available in molecular clouds. Thus, only the most stable conformer (trans) is expected to exist in detectable amounts. We report the first interstellar detection of cis-HCOOH. Its presence in ultraviolet (UV) irradiat…
▽ More
As many organic molecules, formic acid (HCOOH) has two conformers (trans and cis). The energy barrier to internal conversion from trans to cis is much higher than the thermal energy available in molecular clouds. Thus, only the most stable conformer (trans) is expected to exist in detectable amounts. We report the first interstellar detection of cis-HCOOH. Its presence in ultraviolet (UV) irradiated gas exclusively (the Orion Bar photodissociation region), with a low trans-to-cis abundance ratio of 2.8+-1.0, supports a photoswitching mechanism: a given conformer absorbs a stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we specifically study with ab initio quantum calculations, was not considered in Space before but likely induces structural changes of a variety of interstellar molecules submitted to UV radiation.
△ Less
Submitted 3 November, 2016;
originally announced November 2016.
-
Theory of Gas Phase Scattering and Reactivity for Astrochemistry
Authors:
Laurent Wiesenfeld,
Wing-Fai Thi,
Paola Caselli,
Alexandre Faure,
Luca Bizzocchi,
João Brandão,
Denis Duflot,
Eric Herbst,
Stephen J. Klippenstein,
Tamiki Komatsuzaki,
Cristina Puzzarini,
Octavio Roncero,
Hiroshi Teramoto,
Mikito Toda,
Ad van der Avoird,
Holger Waalkens
Abstract:
Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong.
This rep…
▽ More
Because of the very peculiar conditions of chemistry in many astrophysical gases (low densities, mostly low temperatures, kinetics-dominated chemical evolution), great efforts have been devoted to study molecular signatures and chemical evolution. While experiments are being performed in many laboratories, it appears that the efforts directed towards theoretical works are not as strong.
This report deals with the present status of chemical physics/physical chemistry theory, for the qualitative and quantitative understanding of kinetics of molecular scattering, being it reactive or inelastic. By gathering several types of expertise, from applied mathematics to physical chemistry, dialog is made possible, as a step towards new and more adapted theoretical frameworks, capable of meeting the theoretical, methodological and numerical challenges of kinetics-dominated gas phase chemistry in astrophysical environments.
A state of the art panorama is presented, alongside present-day strengths and shortcomings. However, coverage is not complete, being limited in this report to actual attendance of the workshop. Some paths towards relevant progress are proposed.
△ Less
Submitted 3 October, 2016;
originally announced October 2016.
-
Ionization fraction and the enhanced sulfur chemistry in Barnard 1
Authors:
A. Fuente,
J. Cernicharo,
E. Roueff,
M. Gerin,
J. Pety,
N. Marcelino,
R. Bachiller,
B. Lefloch,
O. Roncero,
A. Aguado
Abstract:
Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC). Our aim was to determine the cosmic ray ionization rate and the deplet…
▽ More
Barnard B1b has revealed as one of the most interesting globules from the chemical and dynamical point of view. It presents a rich molecular chemistry characterized by large abundances of deuterated and complex molecules. Furthermore, it hosts an extremely young Class 0 object and one candidate to First Hydrostatic Core (FHSC). Our aim was to determine the cosmic ray ionization rate and the depletion factors in this extremely young star forming region. We carried out a spectral survey towards Barnard 1b as part of the IRAM Large program ASAI using the IRAM 30-m telescope at Pico Veleta (Spain). This provided a very complete inventory of neutral and ionic C-, N- and S- bearing species with, up to our knowledge, the first secure detections of the deuterated ions DCS+ and DOCO+. We used a state-of-the-art pseudo-time-dependent gas-phase chemical model to determine the value of the cosmic ray ionization rate and the depletion factors. The observational data were well fitted with $ζ_{H_2}$ between 3E-17 s$^{-1}$ and 1E-16 s$^{-1}$. Elemental depletions were estimated to be ~10 for C and O, ~1 for N and ~25 for S. Barnard B1b presents similar depletions of C and O than those measured in pre-stellar cores. The depletion of sulfur is higher than that of C and O but not as extreme as in cold cores. In fact, it is similar to the values found in some bipolar outflows, hot cores and photon-dominated regions. Several scenarios are discussed to account for these peculiar abundances. We propose that it is the consequence of the initial conditions (important outflows and enhanced UV fields in the surroundings) and a rapid collapse (~0.1 Myr) that permits to maintain most S- and N-bearing species in gas phase to great optical depths. The interaction of the compact outflow associated with B1b-S with the surrounding material could enhance the abundances of S-bearing molecules, as well.
△ Less
Submitted 27 May, 2016; v1 submitted 16 May, 2016;
originally announced May 2016.
-
OH+ in astrophysical media: state-to-state formation rates, Einstein coefficients and inelastic collision rates with He
Authors:
S. Gomez-Carrasco,
B. Godard,
F. Lique,
N. Bulut,
J. Klos,
O. Roncero,
A. Aguado,
F. J. Aoiz,
J. F. Castillo,
J. R. Goicoechea,
M. Etxaluze,
J. Cernicharo
Abstract:
The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods.
First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3Σ^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-i…
▽ More
The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods.
First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3Σ^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-independent results to asses the accuracy of reaction probabilities at collision energies of about 1 meV. The good agreement between the simulations and the existing experimental cross sections in the $0.01-$1 eV energy range shows the quality of the results.
The calculated state-to-state rate constants have been fitted to an analytical form. Second, the Einstein coefficients of OH$^+$ have been obtained for all astronomically significant ro-vibrational bands involving the $X^3Σ^-$ and/or $A^3Π$ electronic states.
For this purpose the potential energy curves and electric dipole transition moments for seven electronic states of OH$^+$ are calculated with {\it ab initio} methods at the highest level and including spin-orbit terms, and the rovibrational levels have been calculated including the empirical spin-rotation and spin-spin terms. Third, the state-to-state rate constants for inelastic collisions between He and OH$^+(X ^3Σ^-)$ have been calculated using a time-independent close coupling method on a new potential energy surface. All these rates have been implemented in detailed chemical and radiative transfer models. Applications of these models to various astronomical sources show that inelastic collisions dominate the excitation of the rotational levels of OH$^+$. In the models considered the excitation resulting from the chemical formation of OH$^+$ increases the line fluxes by about 10 % or less depending on the density of the gas.
△ Less
Submitted 16 May, 2014;
originally announced May 2014.
-
Depolarizing collisions with hydrogen: neutral and singly ionized alkaline earths
Authors:
Rafael Manso Sainz,
Octavio Roncero,
Cristina Sanz-Sanz,
Alfredo Aguado,
Andres Asensio Ramos,
Javier Trujillo Bueno
Abstract:
Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkaly-earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisio…
▽ More
Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkaly-earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisions with H. We compute ab initio potential curves of the atom-H system and solve the quantum mechanical dynamics. From the scattering amplitudes we calculate the depolarizing rates for Maxwellian distributions of colliders at temperatures T <10000 K. A comparative analysis of our results and previous calculations in the literature is done. We discuss the effect of these rates on the formation of scattering polarization patterns of resonant lines of alkali-earths in the solar atmosphere, and their effect on Hanle effect diagnostics of solar magnetic fields.
△ Less
Submitted 29 April, 2014; v1 submitted 25 April, 2014;
originally announced April 2014.
-
Electron transport signature of H$_2$ dissociation on atomic gold wires
Authors:
Alexandre Zanchet,
Anaís Dorta-Urra,
Octavio Roncero,
Alfredo Aguado,
José Ignacio Martínez,
Fernando Flores,
Nicolás Lorente
Abstract:
Non-equilibrium Green's functions calculations based on density functional theory show a direct link between the initial stages of H$_2$ dissociation on a gold atomic wire and the electronic current supported by the gold wire. The simulations reveal that for biases below the stability threshold of the wire, the minimum-energy path for H$_2$ dissociation is not affected. However, the electronic cur…
▽ More
Non-equilibrium Green's functions calculations based on density functional theory show a direct link between the initial stages of H$_2$ dissociation on a gold atomic wire and the electronic current supported by the gold wire. The simulations reveal that for biases below the stability threshold of the wire, the minimum-energy path for H$_2$ dissociation is not affected. However, the electronic current presents a dramatic drop when the molecule initiates its dissociation. This current drop is traced back to quantum interference between electron paths when the molecule starts interacting with the gold wire.
△ Less
Submitted 26 December, 2013;
originally announced December 2013.
-
Quantum Zeno effect: Quantum shuffling and Markovianity
Authors:
A. S. Sanz,
C. Sanz-Sanz,
T. Gonzalez-Lezana,
O. Roncero,
S. Miret-Artes
Abstract:
The behavior displayed by a quantum system when it is perturbed by a series of von Neumann measurements along time is analyzed. Because of the similarity between this general process with giving a deck of playing cards a shuffle, here it is referred to as quantum shuffling, showing that the quantum Zeno and anti-Zeno effects emerge naturally as two time limits. Within this framework, a connection…
▽ More
The behavior displayed by a quantum system when it is perturbed by a series of von Neumann measurements along time is analyzed. Because of the similarity between this general process with giving a deck of playing cards a shuffle, here it is referred to as quantum shuffling, showing that the quantum Zeno and anti-Zeno effects emerge naturally as two time limits. Within this framework, a connection between the gradual transition from anti-Zeno to Zeno behavior and the appearance of an underlying Markovian dynamics is found. Accordingly, although a priori it might result counterintuitive, the quantum Zeno effect corresponds to a dynamical regime where any trace of knowledge on how the unperturbed system should evolve initially is wiped out (very rapid shuffling). This would explain why the system apparently does not evolve or decay for a relatively long time, although it eventually undergoes an exponential decay. By means of a simple working model, conditions characterizing the shuffling dynamics have been determined, which can be of help to understand and to devise quantum control mechanisms in a number of processes from the atomic, molecular and optical physics.
△ Less
Submitted 3 February, 2012; v1 submitted 16 December, 2011;
originally announced December 2011.
-
Revisiting Zeno's paradox: Flying arrows for atom-diatom reactions
Authors:
C. Sanz-Sanz,
A. S. Sanz,
T. Gonzalez-Lezana,
O. Roncero,
S. Miret-Artes
Abstract:
The possibility to observe quantum Zeno and anti-Zeno scenarios for atom-diatom reactive collisions is investigated for two diferent processes (F+HD and H+O_2) by means of time-dependent wave packet propagations. A novel approach is proposed in which the survival probabilities investigated are those obtained when the initial state is redefined after time intervals tau at which measurements are per…
▽ More
The possibility to observe quantum Zeno and anti-Zeno scenarios for atom-diatom reactive collisions is investigated for two diferent processes (F+HD and H+O_2) by means of time-dependent wave packet propagations. A novel approach is proposed in which the survival probabilities investigated are those obtained when the initial state is redefined after time intervals tau at which measurements are performed on the system. The comparison with the actual probability for the unperturbed system reveals the existence of a regime in which the decay from the initial state appears to be inhibited (the quantum Zeno effect) or accelerated (the anti-Zeno effect). In contrast with preceding interpretations, given that the time-evolving wave packet is not affected at any time ("flying arrow"), the present procedure does not require to invoke to counterintuitive hindered evolutions ("stopping arrow") to explain the results. On the contrary the consecutive update of the reference state of the system solves the apparent modern version of Zeno's paradox.
△ Less
Submitted 7 November, 2011;
originally announced November 2011.
-
Quantum Zeno-based control mechanism for molecular fragmentation
Authors:
C. Sanz-Sanz,
A. S. Sanz,
T. Gonzalez-Lezana,
O. Roncero,
S. Miret-Artes
Abstract:
A quantum control mechanism is proposed for molecular fragmentation processes within a scenario grounded on the quantum Zeno effect. In particular, we focus on the van der Waals Ne-Br$_2$ complex, which displays two competing dissociation channels via vibrational and electronic predissociation. Accordingly, realistic three dimensional wave packet simulations are carried out by using ab initio inte…
▽ More
A quantum control mechanism is proposed for molecular fragmentation processes within a scenario grounded on the quantum Zeno effect. In particular, we focus on the van der Waals Ne-Br$_2$ complex, which displays two competing dissociation channels via vibrational and electronic predissociation. Accordingly, realistic three dimensional wave packet simulations are carried out by using ab initio interaction potentials recently obtained to reproduce available experimental data. Two numerical models to simulate the repeated measurements are reported and analyzed. It is found that the otherwise fast vibrational predissociation is slowed down in favor of the slow electronic (double fragmentation) predissociation, which is enhanced by several orders of magnitude. Based on these theoretical predictions, some hints to experimentalists to confirm their validity are also proposed.
△ Less
Submitted 22 March, 2012; v1 submitted 21 June, 2011;
originally announced June 2011.
