-
Rotational transitions induced by collisions of HD$^{+}$ ions with low energy electrons
Authors:
O. Motapon,
N. Pop,
F. Argoubi,
J. Zs. Mezei,
M. D. Epée Epée,
A. Faure,
M. Telmini,
J. Tennyson,
I. F. Schneider
Abstract:
A series of Multichannel Quantum Defect Theory-based computations have been performed, in order to produce the cross sections of rotational transitions (excitations $N_{i}^{+}-2 \rightarrow$ $N_{i}^{+}$, de-excitations $N_{i}^{+}$ $\rightarrow$ $N_{i}^{+}-2$, with $N_{i}^{+}=2$ to $10$) and of their competitive process, the dissociative recombination, induced by collisions of HD$^+$ ions with elec…
▽ More
A series of Multichannel Quantum Defect Theory-based computations have been performed, in order to produce the cross sections of rotational transitions (excitations $N_{i}^{+}-2 \rightarrow$ $N_{i}^{+}$, de-excitations $N_{i}^{+}$ $\rightarrow$ $N_{i}^{+}-2$, with $N_{i}^{+}=2$ to $10$) and of their competitive process, the dissociative recombination, induced by collisions of HD$^+$ ions with electrons in the energy range $10^{-5}$ to 0.3 eV. Maxwell anisotropic rate coefficients, obtained from these cross sections in the conditions of the Heidelberg Test Storage Ring (TSR) experiments ($k_{B}T_{t}=2.8$ meV and $k_{B}T_{l}=45$ $μ$eV), have been reported for those processes in the same electronic energy range. Maxwell isotropic rate coefficients have been as well presented for electronic temperatures up to a few hundreds of Kelvins. Very good overall agreement is found between our results for rotational transitions and the former theoretical computations as well as with experiment. Furthermore, owing to the full rotational computations performed, the accuracy of the resulting dissociative recombination cross sections is considerably improved.
△ Less
Submitted 10 May, 2024;
originally announced May 2024.
-
Dissociative recombination of the CH$^+$ molecular ion at low energy
Authors:
K. Chakrabarti,
J. Zs Mezei,
O. Motapon,
A. Faure,
O. Dulieu,
K. Hassouni,
I. F. Schneider
Abstract:
The reactive collisions of the CH$^+$ molecular ion with electrons is studied in the framework of the multichannel quantum defect theory, taking into account the contribution of the core-excited Rydberg states. In addition to the $X^1Σ^+$ ground state of the ion, we also consider the contribution to the dynamics of the $a^3Π$ and $A^1Π$ excited states of CH$^+$. Our results - in the case of the di…
▽ More
The reactive collisions of the CH$^+$ molecular ion with electrons is studied in the framework of the multichannel quantum defect theory, taking into account the contribution of the core-excited Rydberg states. In addition to the $X^1Σ^+$ ground state of the ion, we also consider the contribution to the dynamics of the $a^3Π$ and $A^1Π$ excited states of CH$^+$. Our results - in the case of the dissociative recombination in good agreement with the storage ring measurements - rely on decisive improvements - complete account of the ionisation channels and accurate evaluation of the reaction matrix - of a previously used model.
△ Less
Submitted 30 April, 2024;
originally announced April 2024.
-
Maser Activity of Organic Molecules toward Sgr B2(N)
Authors:
Ci Xue,
Anthony Remijan,
Alexandre Faure,
Emmanuel Momjian,
Todd R. Hunter,
Ryan A. Loomis,
Eric Herbst,
Brett McGuire
Abstract:
At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH$_2$NH, HNCNH, and HCOOCH$_3$. However, the lack of spatial distribution information of these new maser species has prevented us from assessing the excitation conditions of the maser emission as…
▽ More
At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH$_2$NH, HNCNH, and HCOOCH$_3$. However, the lack of spatial distribution information of these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 North (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH$_3$OH) maser at 84 GHz and the first interferometric map of the methanimine (CH$_2$NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov-Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH$_3$OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH$_2$NH maser has a close spatial relationship with the 84 GHz CH$_3$OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions.
△ Less
Submitted 19 April, 2024;
originally announced April 2024.
-
NLTE modelling of water-rich exoplanet atmospheres. Cooling and heating rates
Authors:
A. García Muñoz,
A. Asensio Ramos,
A. Faure
Abstract:
The hydrogen and water molecules respond very differently to the collisional-radiative processes taking place in planetary atmospheres. Naturally, the question arises whether H2O-rich atmospheres are more (or less) resilient to long-term mass loss than H2-dominated ones if they radiate away the incident stellar energy more (or less) efficiently. If confirmed, the finding would have implications on…
▽ More
The hydrogen and water molecules respond very differently to the collisional-radiative processes taking place in planetary atmospheres. Naturally, the question arises whether H2O-rich atmospheres are more (or less) resilient to long-term mass loss than H2-dominated ones if they radiate away the incident stellar energy more (or less) efficiently. If confirmed, the finding would have implications on our understanding of the evolution of sub-Neptune exoplanets. As a key step towards answering this question, we present a non-local thermodynamic equilibrium (NLTE) model of H2O for the atmospheric region where the gas accelerates to escape the planet and conditions relevant to close-in sub-Neptunes. Our exploratory calculations for isothermal gas composed of H2, H2O and e- reveal that: 1) In the pressure region ~1e-2 - 1e-4 dyn cm-2 where the stellar extreme-ultraviolet (XUV) photons are typically deposited in the atmosphere, H2O is in rotational LTE but vibrational NLTE. Vibrational LTE is facilitated by high H2O abundances and fractional ionizations, and we report critical densities for the LTE-NLTE transition; 2) Vibrational cooling may locally dominate over rotational cooling, partly because of the comparatively small opacities of ro-vibrational lines; 3) Even low H2O abundances notably enhance the cooling, foreseeably offsetting some of the stellar heating; 4) Heating due to the deposition of stellar infrared (IR) photons is significant at pressures >=0.1 dyn cm-2. We estimate the contribution of H2O excitation to the internal energy of the gas and speculate on the photodissociation from the excited vibrational states. Ultimately, our findings motivate the consideration of NLTE in the mass loss rate calculations of H2O-rich atmospheres.
△ Less
Submitted 3 April, 2024;
originally announced April 2024.
-
Anatomy of the Class I protostar L1489 IRS with NOEMA -- I. Disk, streamers, outflow(s) and bubbles at 3mm
Authors:
M. Tanious,
R. Le Gal,
R. Neri,
A. Faure,
A. Gupta,
C. J. Law,
J. Huang,
N. Cuello,
J. P. Williams,
F. Ménard
Abstract:
Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed stru…
▽ More
Over the past few years, chemical studies have revealed multiple structures in the vicinity of young stellar objects (YSOs). It has become evident that specific physical conditions are associated with the emission of particular molecular lines, allowing us to use molecular probes of the YSO physics. Consequently, chemical surveys are now necessary to fully constrain the origin of the observed structures. Several surveys have been conducted to explore the chemistry of YSOs, focusing on Class 0 and Class II objects. However, our knowledge of intermediate objects, that are Class I objects, remains limited. To bridge the gap and establish the relationship between observed structures and molecular line emission at the Class I evolutionary stage, we investigate the spatial distribution of key molecular gas species in the low-mass Class I protostar L1489 IRS (IRAS 04016+2610), a source part of the ChemYSO survey. We performed a 3mm line survey at high spatial and high spectral resolution using the NOEMA interferometer and the IRAM-30m telescope. We present here the ten brightest lines of our survey, in which we identified a new ~ 3 000 au long streamer in HC3N, C2H, and c-C3H2 emission, likely associated with more localized accretion shocks probed in SO. In addition, two ~ 10 000 au bubbles are seen with the dense molecular tracers HCO+, CS, and HCN around the YSO. Additionally, potential indicators of a second outflow appear in CS and HCN emission, but its nature remains to be confirmed. The late infall identified at large scales may originate from the nearby prestellar core L1489 and is likely responsible for the formation of an external warped disk in this system. The detection of a potential second outflow could be the direct evidence of a binary system. Finally, we hypothesize that the bubbles may result from the magnetic pressure as observed in numerical simulations.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
H2CO and CS in diffuse clouds: Excitation and abundance
Authors:
Maryvonne Gerin,
Harvey Liszt,
Jerome Pety,
Alexandre Faure
Abstract:
To provide constraints on the chemical processes responsible for the observed columns of organic species, we used NOEMA to observe the sight line toward NRAO150 in the 2mm spectral window. We targeted the low excitation lines of o-H2CO 2(1,1)-1(1,0) and p-H2CO 2(0,2)-1(0,1) as well as the nearby transitions of CS(3-2) and c-C3H2. We combined these data with previous observations to determine the e…
▽ More
To provide constraints on the chemical processes responsible for the observed columns of organic species, we used NOEMA to observe the sight line toward NRAO150 in the 2mm spectral window. We targeted the low excitation lines of o-H2CO 2(1,1)-1(1,0) and p-H2CO 2(0,2)-1(0,1) as well as the nearby transitions of CS(3-2) and c-C3H2. We combined these data with previous observations to determine the excitation conditions, column densities, and abundances relative to H2 in the different velocity components. We performed non-LTE radiative transfer calculations including collision cross sections with ortho and para H2 and with electrons. New collision cross sections with electrons were computed for ortho and para formaldehyde. The c-C3H2 line profiles are very similar to those of HCO+ and CCH, while the CS absorption features are narrower and mostly concentrated in two main velocity components at V = -17 and -10 km/s. H2CO absorption lines present an intermediate pattern with absorption in all velocity components but larger opacities in the two main velocity components. The ortho-to-para ratios of H2CO and c-C3H2 are consistent with the statistical value of 3. While the excitation temperature of all c-C3H2 velocity components is consistent with the CMB, the two strong components detected in CS show a clear excess over the CMB indicating that CS resides at higher densities than other species along this particular sightline, n(H2) ~ 2500 cm-3 while n(H2) < 500 cm-3 for the other velocity components. We detected faint absorption from o-H213CO and C34S allowing us to derive isotopic ratios: o-H2CO/o-H213CO = 61 and C32S/C34S = 24. The excitation of the 4.8GHz line of formaldehyde is sensitive to the electron fraction and its excitation temperature is predicted to be lower than the CMB at low and moderate electron fractions, x(e)< 6E-5, and to rise above the CMB at high electron fractions, > 1e-4.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Protonated acetylene in the z=0.89 molecular absorber toward PKS1830-211
Authors:
S. Muller,
R. Le Gal,
E. Roueff,
J. H. Black,
A. Faure,
M. Guelin,
A. Omont,
M. Gerin,
F. Combes,
S. Aalto
Abstract:
We report the first interstellar identification of protonated acetylene, C2H3+, a fundamental hydrocarbon, in the z=0.89 molecular absorber toward the gravitationally lensed quasar PKS1830-211. The molecular species is identified from clear absorption features corresponding to the 2_12-1_01 (rest frequency 494.034 GHz) and 1_11-0_00 (431.316 GHz) ground-state transitions of ortho and para forms of…
▽ More
We report the first interstellar identification of protonated acetylene, C2H3+, a fundamental hydrocarbon, in the z=0.89 molecular absorber toward the gravitationally lensed quasar PKS1830-211. The molecular species is identified from clear absorption features corresponding to the 2_12-1_01 (rest frequency 494.034 GHz) and 1_11-0_00 (431.316 GHz) ground-state transitions of ortho and para forms of C2H3+, respectively, in ALMA spectra toward the southwestern image of PKS1830-211, where numerous molecules, including other hydrocarbons, have already been detected. From the simple assumption of local thermodynamic equilibrium (LTE) with cosmic microwave background photons and an ortho-to-para ratio of three, we estimate a total C2H3+ column density of 2 x 10^12 cm^-2 and an abundance of 10^-10 compared to H_2. However, formation pumping could affect the population of metastable states, yielding a C2H3+ column density higher than the LTE value by a factor of a few. We explore possible routes to the formation of C2H3+, mainly connected to acetylene and methane, and find that the methane route is more likely in PDR environment. As one of the initial hydrocarbon building blocks, C2H3+ is thought to play an important role in astrochemistry, in particular in the formation of more complex organic molecules.
△ Less
Submitted 18 January, 2024;
originally announced January 2024.
-
Identification of Ammonium Salts on Comet 67P/C-G Surface from Infrared VIRTIS/Rosetta Data Based on Laboratory Experiments. Implications and Perspectives
Authors:
Olivier Poch,
Istiqomah Istiqomah,
Eric Quirico,
Pierre Beck,
Bernard Schmitt,
Patrice Theulé,
Alexandre Faure,
Pierre Hily-Blant,
Lydie Bonal,
Andrea Raponi,
Mauro Ciarniello,
Batiste Rousseau,
Sandra Potin,
Olivier Brissaud,
Laurène Flandinet,
Gianrico Filacchione,
Antoine Pommerol,
Nicolas Thomas,
David Kappel,
Vito Mennella,
Lyuba Moroz,
Vassilissa Vinogradoff,
Gabriele Arnold,
Stéphane Erard,
Dominique Bockelée-Morvan
, et al. (7 additional authors not shown)
Abstract:
The nucleus of comet 67P/Churyumov-Gerasimenko exhibits a broad spectral reflectance feature around 3.2 $μ$m, which is omnipresent in all spectra of the surface, and whose attribution has remained elusive since its discovery. Based on laboratory experiments, we have shown that most of this absorption feature is due to ammonium (NH4+) salts mixed with the dark surface material. The depth of the ban…
▽ More
The nucleus of comet 67P/Churyumov-Gerasimenko exhibits a broad spectral reflectance feature around 3.2 $μ$m, which is omnipresent in all spectra of the surface, and whose attribution has remained elusive since its discovery. Based on laboratory experiments, we have shown that most of this absorption feature is due to ammonium (NH4+) salts mixed with the dark surface material. The depth of the band is compatible with semi-volatile ammonium salts being a major reservoir of nitrogen in the comet, which could dominate over refractory organic matter and volatile species. These salts may thus represent the long-sought reservoir of nitrogen in comets, possibly bringing their nitrogen-to-carbon ratio in agreement with the solar value. Moreover, the reflectance spectra of several asteroids are compatible with the presence of NH4+ salts at their surfaces. The presence of such salts, and other NH4+-bearing compounds on asteroids, comets, and possibly in proto-stellar environments, suggests that NH4+ may be a tracer of the incorporation and transformation of nitrogen in ices, minerals and organics, at different phases of the formation of the Solar System.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Chasing Gravitational Waves with the Cherenkov Telescope Array
Authors:
Jarred Gershon Green,
Alessandro Carosi,
Lara Nava,
Barbara Patricelli,
Fabian Schüssler,
Monica Seglar-Arroyo,
Cta Consortium,
:,
Kazuki Abe,
Shotaro Abe,
Atreya Acharyya,
Remi Adam,
Arnau Aguasca-Cabot,
Ivan Agudo,
Jorge Alfaro,
Nuria Alvarez-Crespo,
Rafael Alves Batista,
Jean-Philippe Amans,
Elena Amato,
Filippo Ambrosino,
Ekrem Oguzhan Angüner,
Lucio Angelo Antonelli,
Carla Aramo,
Cornelia Arcaro,
Luisa Arrabito
, et al. (545 additional authors not shown)
Abstract:
The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very…
▽ More
The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA.
△ Less
Submitted 5 February, 2024; v1 submitted 11 October, 2023;
originally announced October 2023.
-
Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array
Authors:
The Cherenkov Telescope Array Consortium,
:,
K. Abe,
S. Abe,
F. Acero,
A. Acharyya,
R. Adam,
A. Aguasca-Cabot,
I. Agudo,
A. Aguirre-Santaella,
J. Alfaro,
R. Alfaro,
N. Alvarez-Crespo,
R. Alves Batista,
J. -P. Amans,
E. Amato,
E. O. Angüner,
L. A. Antonelli,
C. Aramo,
M. Araya,
C. Arcaro,
L. Arrabito,
K. Asano,
Y. Ascasíbar,
J. Aschersleben
, et al. (542 additional authors not shown)
Abstract:
Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med…
▽ More
Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
Hyperfine collisional excitation of ammonia by molecular hydrogen
Authors:
J. Loreau,
A. Faure,
F. Lique,
S. Demes,
P. J. Dagdigian
Abstract:
Ammonia is one of the most widely observed molecules in space, and many observations are able to resolve the hyperfine structure due to the electric quadrupole moment of the nitrogen nucleus. The observed spectra often display anomalies in the satellite components of the lines, which indicate substantial deviations from the local thermodynamic equilibrium. The interpretation of the spectra thus re…
▽ More
Ammonia is one of the most widely observed molecules in space, and many observations are able to resolve the hyperfine structure due to the electric quadrupole moment of the nitrogen nucleus. The observed spectra often display anomalies in the satellite components of the lines, which indicate substantial deviations from the local thermodynamic equilibrium. The interpretation of the spectra thus requires the knowledge of the rate coefficients for the hyperfine excitation of NH$_3$ induced by collisions with H$_2$ molecules, the dominant collider in the cold interstellar medium. In this paper we present the first such calculations using a recoupling approach. The rate coefficients are obtained for all hyperfine levels within rotation-inversion levels up to $j=4$ and temperatures up to 100 K by means of quantum scattering close-coupling calculations on an accurate, five-dimensional, potential energy surface. We show that the rate coefficients depart significantly from those obtained with the statistical approach and that they do not conform to any simple propensity rules. Finally, we perform radiative transfer calculations to illustrate the impact of our new rate coefficients by modelling the hyperfine line intensities of the inversion transition in ground state para-NH$_3$ ($j_k=1_1$) and of the rotational transition $1_0\rightarrow 0_0$ in ortho-NH$_3$.
△ Less
Submitted 4 September, 2023;
originally announced September 2023.
-
An accurate set of H$_3$O$^+ -$ H$_2$ collisional rate coefficients for non-LTE modelling of warm interstellar clouds
Authors:
Sándor Demes,
François Lique,
Alexandre Faure,
Floris F. S. van der Tak
Abstract:
Hydronium (H$_3$O$^+$) was first detected in 1986 in interstellar molecular clouds. It was reported in many galactic diffuse and dense regions, as well as in extragalactic sources. H$_3$O$^+$ plays a major role both in interstellar oxygen and water chemistry. However, despite the large number of H$_3$O$^+$ observations, its collisional excitation was investigated only partially. In the present wor…
▽ More
Hydronium (H$_3$O$^+$) was first detected in 1986 in interstellar molecular clouds. It was reported in many galactic diffuse and dense regions, as well as in extragalactic sources. H$_3$O$^+$ plays a major role both in interstellar oxygen and water chemistry. However, despite the large number of H$_3$O$^+$ observations, its collisional excitation was investigated only partially. In the present work we study the state-to-state rotational de-excitation of $ortho$- and $para$-H$_3$O$^+$ in collisions both with $ortho$- and $para$-H$_2$. The cross sections are calculated within the close-coupling formalism using a highly accurate potential energy surface developed for this system. The rate coefficients are computed up to $300$ K kinetic temperature. Transitions between the lowest 21 rotation-inversion states were studied for $para$-H$_3$O$^+$, and the lowest 11 states for $ortho$-H$_3$O$^+$, i.e. all levels with rotational energies below 430 K ($\sim 300$ cm$^{-1}$) are considered (up to $j\leq5$). In order to estimate the impact of the new rate coefficients on the astrophysical models for H$_3$O$^+$, radiative transfer calculations were also carried out. We have examined how the new collisional data affect the line intensities with respect to older data previously used for the interpretation of observations. By analysing all detected transitions we find that our new, accurate rate coefficients have a significant impact (typically within a factor of 2) on radiation temperatures, allowing more accurate estimation of column densities and relative abundances of hydronium, especially in warm molecular clouds, paving the path towards better interpretation of interstellar water and oxygen chemistry.
△ Less
Submitted 4 November, 2022;
originally announced November 2022.
-
The effect of water and electron collisions in the rotational excitation of HF in comets
Authors:
J. Loreau,
A. Faure,
F. Lique
Abstract:
We present the first set of rate coefficients for the rotational excitation of the 7 lowest levels of hydrogen fluoride (HF) induced by collision with water molecules, the dominant collider in cometary comas, in the 5-150 K temperature range. The calculations are performed with a quantum statistical approach from an accurate rigid rotor ab initio interaction potential. Rate coefficients for excita…
▽ More
We present the first set of rate coefficients for the rotational excitation of the 7 lowest levels of hydrogen fluoride (HF) induced by collision with water molecules, the dominant collider in cometary comas, in the 5-150 K temperature range. The calculations are performed with a quantum statistical approach from an accurate rigid rotor ab initio interaction potential. Rate coefficients for excitation of HF by electron-impact are also computed, within the Born approximation, in the 10-10,000 K temperature range. These rate coefficients are then used in a simplified non-local thermodynamic equilibrium (non-LTE) model of a cometary coma that also includes solar radiative pumping and radiative decay. We investigate the range of H2O densities that lead to non-LTE populations of the rotational levels of HF. We show that to describe the excitation of HF in comets, considering collisions with both water molecules and electrons is needed as a result of the large dipole of HF.
△ Less
Submitted 21 October, 2022;
originally announced October 2022.
-
The ortho-to-para ratio of water in interstellar clouds
Authors:
A. Faure,
P. Hily-Blant,
C. Rist,
G. Pineau des Forêts,
A. Matthews,
D. R. Flower
Abstract:
The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involv…
▽ More
The nuclear-spin chemistry of interstellar water is investigated using the University of Grenoble Alpes Astrochemical Network (UGAN). This network includes reactions involving the different nuclear-spin states of the hydrides of carbon, nitrogen, oxygen and sulphur, as well as their deuterated forms. Nuclear-spin selection rules are implemented within the scrambling hypothesis for reactions involving up to seven protons. The abundances and ortho-to-para ratios (OPRs) of gas-phase water and water ions (H$_2$O$^+$ and H$_3$O$^+$) are computed under the steady-state conditions representative of a dark molecular cloud and during the early phase of gravitational collapse of a prestellar core. The model incorporates the freezing of the molecules on to grains, simple grain surface chemistry and cosmic-ray induced and direct desorption of ices. The predicted OPRs are found to deviate significantly from both thermal and statistical values and to be independent of temperature below $\sim $30~K. The OPR of H$_2$O is shown to lie between 1.5 and 2.6, depending on the spin-state of H$_2$, in good agreement with values derived in translucent clouds with relatively high extinction. In the prestellar core collapse calculations, the OPR of H$_2$O is shown to reach the statistical value of 3 in regions with severe depletion ($n_{\rm H}> 10^7$~cm$^{-3}$). We conclude that a low water OPR ($\lesssim 2.5$) is consistent with gas-phase ion-neutral chemistry and reflects a gas with OPR(H$_2)\lesssim 1$. Available OPR measurements in protoplanetary disks and comets are finally discussed.
△ Less
Submitted 6 January, 2022;
originally announced January 2022.
-
Collisional excitation and non-LTE modelling of interstellar chiral propylene oxide
Authors:
K. Dzenis,
A. Faure,
B. A. McGuire,
A. J. Remijan,
P. J. Dagdigian,
C. Rist,
R. Dawes,
E. Quintas-Sanchez,
F. Lique,
M. Hochlaf
Abstract:
The first set of theoretical cross sections for propylene oxide (CH3CHCH2O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H2 are deduced in the temperature range 5-30 K. These collisional coefficients are combined with radiative dat…
▽ More
The first set of theoretical cross sections for propylene oxide (CH3CHCH2O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H2 are deduced in the temperature range 5-30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made towards the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (~ 10 K) and translucent (nH ~ 2000 cm-3) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is Ntot ~ 3e12 cm-2, corresponding to a fractional abundance relative to total hydrogen of ~ 2.5e-11. The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.
△ Less
Submitted 5 February, 2022; v1 submitted 16 December, 2021;
originally announced December 2021.
-
Sulfur gas-phase abundance in dense cores
Authors:
P. Hily-Blant,
G. Pineau des Forêts,
A. Faure,
F. Lique
Abstract:
The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has…
▽ More
The abundance of volatile sulfur in dense clouds is long-standing problem in studies of the physics and chemistry of star-forming regions. Sulfur is an important species because its low ionization potential may possibly make it an important charge carrier. The observed sulfur-bearing species in the gas-phase of dense clouds represent only a minor fraction of the cosmic sulfur abundance, which has been interpreted as a signature of sulfur depletion into ices at the surface of dust grains. However, atomic sulfur, which could be the main gas-phase carrier, cannot be observed directly in cold cores. We present measurements of the nitrogen sulfide (NS) radical toward four dense cores performed with the IRAM-30m telescope. Analytical chemical considerations and chemical models over a wide parameter space show that the NS:N2H+ abundance ratio provides a direct constraint on the abundance of gas-phase atomic sulfur. Toward early-type cores, we find that $n(\rm S)/n_{\rm H}$ is close, or even equal, to the cosmic abundance of sulfur, 14$\times 10^{-6}$, demonstrating that sulfur is not depleted and is atomic, which is in agreement with chemical models. More chemically evolved cores show sulfur depletion by factors up to 100 in their densest parts. In L1544, atomic sulfur depletion is shown to increase with increasing density. Future observations are needed to discover the solid-phase carrier of sulfur. The initial steps of the collapse of pre-stellar cores in the high sulfur abundance regime also need to be explored from their chemical and dynamical perspectives.
△ Less
Submitted 5 December, 2021; v1 submitted 2 December, 2021;
originally announced December 2021.
-
Observations and analysis of CH$^+$ vibrational emissions from the young, carbon-rich planetary nebula NGC 7027: a textbook example of chemical pumping
Authors:
David A. Neufeld,
Benjamin Godard,
P. Bryan Changala,
Alexandre Faure,
T. R. Geballe,
Rolf Güsten,
Karl M. Menten,
Helmut Wiesemeyer
Abstract:
We discuss the detection of 14 rovibrational lines of CH$^+$, obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. Our observations in the 3.49 - 4.13 $μ$m spectral region, obtained with a 0.375" slit width that provided a spectral resolving power $λ/Δλ\sim 80,000$, have resulted in the unequivocal detection of the $R(0) - R(3)$ and $P(1)-P(10)$ transitio…
▽ More
We discuss the detection of 14 rovibrational lines of CH$^+$, obtained with the iSHELL spectrograph on NASA's Infrared Telescope Facility (IRTF) on Maunakea. Our observations in the 3.49 - 4.13 $μ$m spectral region, obtained with a 0.375" slit width that provided a spectral resolving power $λ/Δλ\sim 80,000$, have resulted in the unequivocal detection of the $R(0) - R(3)$ and $P(1)-P(10)$ transitions within the $v=1-0$ band of CH$^+$. The $R$-branch transitions are anomalously weak relative to the $P$-branch transitions, a behavior that is explained accurately by rovibronic calculations of the transition dipole moment reported in a companion paper (Changala et al. 2021). Nine infrared transitions of H$_2$ were also detected in these observations, comprising the $S(8)$, $S(9)$, $S(13)$ and $S(15)$ pure rotational lines; the $v=1-0$ $O(4) - O(7)$ lines, and the $v=2-1$ $O(5)$ line. We present a photodissociation region model, constrained by the CH$^+$ and H$_2$ line fluxes that we measured, that includes a detailed treatment of the excitation of CH$^+$ by inelastic collisions, optical pumping, and chemical ("formation") pumping. The latter process is found to dominate the excitation of the observed rovibrational lines of CH$^+$, and the model is remarkably successful in explaining both the absolute and relative strengths of the CH$^+$ and H$_2$ lines.
△ Less
Submitted 28 May, 2021;
originally announced May 2021.
-
Deuterium fractionation of nitrogen hydrides: detections of NHD and ND$_2$
Authors:
A. Bacmann,
A. Faure,
P. Hily-Blant,
K. Kobayashi,
H. Ozeki,
S. Yamamoto,
L. Pagani,
F. Lique
Abstract:
Although ammonia is an abundant molecule commonly observed towards the dense interstellar medium, it has not yet been established whether its main formation route is from gas-phase ion-molecule reactions or grain-surface hydrogen additions on adsorbed nitrogen atoms. Deuterium fractionation can be used as a tool to constrain formation mechanisms. High abundances of deuterated molecules are routine…
▽ More
Although ammonia is an abundant molecule commonly observed towards the dense interstellar medium, it has not yet been established whether its main formation route is from gas-phase ion-molecule reactions or grain-surface hydrogen additions on adsorbed nitrogen atoms. Deuterium fractionation can be used as a tool to constrain formation mechanisms. High abundances of deuterated molecules are routinely observed in the dense interstellar medium, with the ratio between deuterated molecules and the main isotopologue enhanced by several orders of magnitude with respect to the elemental D/H ratio. In the case of ammonia, the detection of its triply deuterated isotopologue hints at high abundances of the deuterated intermediate nitrogen radicals, ND, NHD and ND$_2$. So far however, only ND has been detected in the interstellar medium. In this paper, to constrain the formation of ammonia, we aim at determining the NHD/NH$_2$ and ND$_2$/NHD abundance ratios, and compare them with the predictions of both pure gas-phase and grain-surface chemical models. We searched for the fundamental rotational transitions of NHD and ND$_2$ towards the class 0 protostar IRAS16293-2422, towards which NH, NH$_2$ and ND had been previously detected. Both NHD and ND$_2$ are detected in absorption towards the source. The relative abundance ratios NH$_2$ : NHD : ND$_2$ are close to 8 : 4 : 1. These ratios can be reproduced by our gas-phase chemical model within a factor of two-three. Statistical ratios as expected from grain-surface chemistry are also consistent with our data. Further investigations of the ortho-to-para ratio in ND$_2$ , both theoretical and observational, could bring new constraints to better understand nitrogen hydride chemistry.
△ Less
Submitted 25 September, 2020;
originally announced September 2020.
-
Depletion and fractionation of nitrogen in collapsing cores
Authors:
P. Hily-Blant,
G. Pineau des Forêts,
A. Faure,
D. R. Flower
Abstract:
Measurements of the nitrogen isotopic ratio in Solar System comets show a constant value, ~140, which is three times lower than the protosolar ratio, a highly significant difference that remains unexplained. Observations of static starless cores at early stages of collapse confirm the theoretical expectation that nitrogen fractionation in interstellar conditions is marginal for most species. Yet,…
▽ More
Measurements of the nitrogen isotopic ratio in Solar System comets show a constant value, ~140, which is three times lower than the protosolar ratio, a highly significant difference that remains unexplained. Observations of static starless cores at early stages of collapse confirm the theoretical expectation that nitrogen fractionation in interstellar conditions is marginal for most species. Yet, observed isotopic ratios in N2H+ are at variance with model predictions. These gaps in our understanding of how the isotopic reservoirs of nitrogen evolve, from interstellar clouds to comets, and, more generally, to protosolar nebulae, may have their origin in missing processes or misconceptions in the chemistry of interstellar nitrogen. So far, theoretical studies of nitrogen fractionation in starless cores have addressed the quasi-static phase of their evolution such that the effect of dynamical collapse on the isotopic ratio is not known. In this paper, we investigate the fractionation of 14N and 15N during the gravitational collapse of a pre-stellar core through gas-phase and grain adsorption and desorption reactions. The initial chemical conditions, which are obtained in steady state after typically a few Myr, show low degrees of fractionation in the gas phase, in agreement with earlier studies. However, during collapse, the differential rate of adsorption of 14N- and 15N-containing species onto grains results in enhanced 15N:14N ratios, in better agreement with the observations. Furthermore, we find differences in the behavior, with increasing density, of the isotopic ratio in different species. We find that the collapse must take place on approximately one free-fall timescale, based on the CO abundance profile in L183 [see the end in the PDF file]
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
The Leiden Atomic and Molecular Database (LAMDA): Current Status, Recent Updates, and Future Plans
Authors:
Floris van der Tak,
François Lique,
Alex Faure,
John Black,
Ewine van Dishoeck
Abstract:
The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of…
▽ More
The Leiden Atomic and Molecular Database (LAMDA) collects spectroscopic information and collisional rate coefficients for molecules, atoms, and ions of astrophysical and astrochemical interest. We describe the developments of the database since its inception in 2005, and outline our plans for the near future. Such a database is constrained both by the nature of its uses and by the availability of accurate data: we suggest ways to improve the synergies among users and suppliers of data. We summarize some recent developments in computation of collisional cross sections and rate coefficients. We consider atomic and molecular data that are needed to support astrophysics and astrochemistry with upcoming instruments that operate in the mid- and far-infrared parts of the spectrum.
△ Less
Submitted 23 April, 2020;
originally announced April 2020.
-
Ammonium salts are a reservoir of nitrogen on a cometary nucleus and possibly on some asteroids
Authors:
O. Poch,
I. Istiqomah,
E. Quirico,
P. Beck,
B. Schmitt,
P. Theulé,
A. Faure,
P. Hily-Blant,
L. Bonal,
A. Raponi,
M. Ciarniello,
B. Rousseau,
S. Potin,
O. Brissaud,
L. Flandinet,
G. Filacchione,
A. Pommerol,
N. Thomas,
D. Kappel,
V. Mennella,
L. Moroz,
V. Vinogradoff,
G. Arnold,
S. Erard,
D. Bockelée-Morvan
, et al. (7 additional authors not shown)
Abstract:
The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this…
▽ More
The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.
△ Less
Submitted 12 March, 2020;
originally announced March 2020.
-
The effect of CO-H2O collisions in the rotational excitation of cometary CO
Authors:
A. Faure,
F. Lique,
J. Loreau
Abstract:
We present the first accurate rate coefficients for the rotational excitation of CO by H2O in the kinetic temperature range 5-100 K. The statistical adiabatic channel method (SACM) is combined with a high-level rigid-rotor CO-H2O intermolecular potential surface. Transitions among the first 11 rotational levels of CO and the first 8 rotational levels of both para-H2O and ortho-H2O are considered.…
▽ More
We present the first accurate rate coefficients for the rotational excitation of CO by H2O in the kinetic temperature range 5-100 K. The statistical adiabatic channel method (SACM) is combined with a high-level rigid-rotor CO-H2O intermolecular potential surface. Transitions among the first 11 rotational levels of CO and the first 8 rotational levels of both para-H2O and ortho-H2O are considered. Our rate coefficients are compared to previous data from the literature and they are also incorporated in a simple non-LTE model of cometary coma including collision-induced transitions, solar radiative pumping and radiative decay. We find that the uncertainties in the collision data have significant influence on the CO population distribution for H2O densities in the range 10^3-10^8 cm^-3. We also show that the rotational distribution of H2O plays an important role in CO excitation (owing to correlated energy transfer in both CO and H2O), while the impact of the ortho-to-para ratio of H2O is found to be negligible.
△ Less
Submitted 12 February, 2020; v1 submitted 24 January, 2020;
originally announced January 2020.
-
Cyanoacetylene in the outflow/hot molecular core G331.512-0.103
Authors:
N. U. Duronea,
L. Bronfman,
E. Mendoza,
M. Merello,
R. Finger,
N. Reyes,
C. Hervías-Caimapo,
A. Faure,
C. E. Cappa,
E. M. Arnal,
J. R. D. Lépine,
I. Kleiner,
L-Ä Nyman
Abstract:
Using APEX-1 and APEX-2 observations, we have detected and studied the rotational lines of the HC$_3$N molecule (cyanoacetylene) in the powerful outflow/hot molecular core G331.512-0.103. We identified thirty-one rotational lines at $J$ levels between 24 and 39; seventeen of them in the ground vibrational state $v$=0 (9 lines corresponding to the main C isotopologue and 8 lines corresponding to th…
▽ More
Using APEX-1 and APEX-2 observations, we have detected and studied the rotational lines of the HC$_3$N molecule (cyanoacetylene) in the powerful outflow/hot molecular core G331.512-0.103. We identified thirty-one rotational lines at $J$ levels between 24 and 39; seventeen of them in the ground vibrational state $v$=0 (9 lines corresponding to the main C isotopologue and 8 lines corresponding to the $^{13}$C isotopologues), and fourteen in the lowest vibrationally excited state $v_7$=1. Using LTE-based population diagrams for the beam-diluted $v$=0 transitions, we determined $T_{\rm exc}$=85$\pm$4 K and $N$(HC$_3$N)=(6.9$\pm$0.8)$\times$10$^{14}$ cm$^{-2}$, while for the beam-diluted $v_7$=1 transitions we obtained $T_{\rm exc}$=89$\pm$10 K and $N$(HC$_3$N)=2$\pm$1$\times$10$^{15}$ cm$^{-2}$. Non-LTE calculations using H$_2$ collision rates indicate that the HC$_3$N emission is in good agreement with LTE-based results. From the non-LTE method we estimated $T_{\rm kin}$ $\simeq$90~K, $n$(H$_2$)$\simeq$2$\times$10$^7$~cm$^{-3}$ for a central core of 6 arcsec in size. A vibrational temperature in the range from 130~K to 145~K was also determined, values which are very likely lower limits. Our results suggest that rotational transitions are thermalized, while IR radiative pumping processes are probably more efficient than collisions in exciting the molecule to the vibrationally excited state $v_7$=1. Abundance ratios derived under LTE conditions for the $^{13}$C isotopologues suggest that the main formation pathway of HC$_3$N is ${\rm C}_2{\rm H}_2 + {\rm CN} \rightarrow {\rm HC}_3{\rm N} + {\rm H}$.
△ Less
Submitted 22 July, 2019;
originally announced July 2019.
-
Deuterated methanol toward NGC 7538-IRS1
Authors:
Juan Ospina-Zamudio,
Cécile Favre,
Marina Kounkel,
Li-Hong Xu,
Justin Neill,
Bertrand Lefloch,
Alexandre Faure,
Edwin Bergin,
Davide Fedele,
Lee Hartmann
Abstract:
We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E…
▽ More
We investigate the deuteration of methanol towards the high-mass star forming region NGC 7538-IRS1. We have carried out a multi-transition study of CH$_3$OH, $^{13}$CH$_3$OH and of the deuterated fllavors, CH$_2$DOH and CH$_3$OD, between 1.0--1.4 mm with the IRAM-30~m antenna. In total, 34 $^{13}$CH$_3$OH, 13 CH$_2$DOH lines and 20 CH$_3$OD lines spanning a wide range of upper-state energies (E$_{up}$) were detected. From the detected transitions, we estimate that the measured D/H does not exceed 1$\%$, with a measured CH$_2$DOH/CH$_3$OH and CH$_3$OD/CH$_3$OH of about (32$\pm$8)$\times$10$^{-4}$ and (10$\pm$4)$\times$10$^{-4}$, respectively. This finding is consistent with the hypothesis of a short-time scale formation during the pre-stellar phase. We find a relative abundance ratio CH$_2$DOH/CH$_3$OD of 3.2 $\pm$ 1.5. This result is consistent with a statistical deuteration. We cannot exclude H/D exchanges between water and methanol if water deuteration is of the order 0.1$\%$, as suggested by recent Herschel observations.
△ Less
Submitted 23 May, 2019;
originally announced May 2019.
-
Efficient methanol production on the dark side of a prestellar core
Authors:
Jorma Harju,
Jaime E. Pineda,
Anton I. Vasyunin,
Paola Caselli,
Stella S. R. Offner,
Alyssa A. Goodman,
Mika Juvela,
Olli Sipilae,
Alexandre Faure,
Romane Le Gal,
Pierre Hily-Blant,
Joao Alves,
Luca Bizzocchi,
Andreas Burkert,
Hope Chen,
Rachel K. Friesen,
Rolf Guesten,
Philip C. Myers,
Anna Punanova,
Claire Rist,
Erik Rosolowsky,
Stephan Schlemmer,
Yancy Shirley,
Silvia Spezzano,
Charlotte Vastel
, et al. (1 additional authors not shown)
Abstract:
We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particul…
▽ More
We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.
△ Less
Submitted 4 May, 2020; v1 submitted 27 March, 2019;
originally announced March 2019.
-
The Dense Gas Fraction in Galactic Center Clouds
Authors:
E. A. C. Mills,
A. Ginsburg,
K. Immer,
J. M Barnes,
L. Wiesenfeld,
A. Faure,
M. R. Morris,
M. A. Requena-Torres
Abstract:
We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least t…
▽ More
We present an analysis of gas densities in the central R=300 parsecs of the Milky Way, focusing on three clouds: GCM-0.02-0.07 (the 50 km/s cloud), GCM-0.13-0.08 (the 20 km/s cloud), and GCM0.25+0.01 (the "Brick"). Densities are determined using observations of the J=(3-2), (4-3), (5-4), (10-9), (18-17), (19-18), (21-20), and (24-23) transitions of the molecule HC3N. We find evidence of at least two excitation regimes for HC3N and constrain the low-excitation component to have a density less than 10^4 cm^-3 and the high-excitation component to have a density between 10^5 and 10^6 cm^-3. This is much less than densities of 10^7 cm^-3 that are found in Sgr B2, the most actively star-forming cloud in the Galactic center. This is consistent with the requirement of a higher density threshold for star formation in the Galactic center than is typical in the Galactic disk. We are also able to constrain the column density of each component in order to determine the mass fraction of 'dense' (n>10^5 cm^-3) gas for these clouds. We find that this is ~15% for all three clouds. Applying the results of our models to ratios of the (10-9) and (3-2) line across the entire central R=300 pc, we find that the fraction of gas with n>10^4 cm^-3 increases inward of a radius of ~140 pc, consistent with the predictions of recent models for the gas dynamics in this region. Our observations show that HC3N is an excellent molecule for probing the density structure of clouds in the Galactic center.
△ Less
Submitted 29 September, 2018;
originally announced October 2018.
-
The nitrogen isotopic ratio of HC3N towards the L1544 prestellar core
Authors:
P. Hily-Blant,
A. Faure,
C. Vastel,
V. Magalhaes,
B. Lefloch,
R. Bachiller
Abstract:
The origin of the heavily fractionated reservoir of nitrogen in comets remains an issue in the theory of their formation and hence of the solar system. Whether the fractionated reservoir traced by comets is inherited from the interstellar cloud or is the product of processes taking place in the protostar, or in the protoplanetary disk, remains unclear. So far, observations of nitrogen isotopic rat…
▽ More
The origin of the heavily fractionated reservoir of nitrogen in comets remains an issue in the theory of their formation and hence of the solar system. Whether the fractionated reservoir traced by comets is inherited from the interstellar cloud or is the product of processes taking place in the protostar, or in the protoplanetary disk, remains unclear. So far, observations of nitrogen isotopic ratios in protostars or prestellar cores have not securely identified such a fractionated reservoir owing to the intrinsic difficulty of direct isotopic ratios measurements. In this article, we report the detection of 5 rotational lines of \ce{HC3N}, {including the weaker components of the hyperfine multiplets}, and two rotational lines of its $^{15}$N isotopologue, towards the L1544 prestellar core. Based on a MCMC/non-LTE multi-line analysis at the hyperfine level, we derive the column densities of \ce{HC3N} ($8.0\pm0.4\tdix{13}$\cc) and \ce{HC3^{15}N} ($2.0\pm0.4\tdix{11}$\cc) and derive an isotopic ratio of 400$\pm$20(1$σ$). This value suggests that \ce{HC3N} is slightly depleted in $^{15}$N in L1544 with respect to the elemental $^{14}$N/$^{15}$N ratio {of $\approx$330} in the present-day local interstellar medium. Our study also stresses the need for radiative calculations at the hyperfine level. Finally, the comparison of the derived ratio with those obtained in CN and HCN in the same core seems to favor CN+C$_2$H$_2$ as the dominant formation route to HC$_3$N. However, uncertainties in the isotopic ratios preclude definitive conclusions.
△ Less
Submitted 12 June, 2018;
originally announced June 2018.
-
Collisional Excitation and Weak Maser Action of Interstellar Methanimine
Authors:
Alexandre Faure,
François Lique,
Anthony J. Remijan
Abstract:
The inelastic cross sections and rate coefficients for the rotational excitation of methanimine (CH$_2$NH) by cold H$_2$ have been determined quantum mechanically based on a new highly correlated five-dimensional potential energy surface. This surface was fitted to more than 60 000 ab initio points with a root mean square error of $\sim$ 1-2 cm$^{-1}$ in the region of the potential well whose dept…
▽ More
The inelastic cross sections and rate coefficients for the rotational excitation of methanimine (CH$_2$NH) by cold H$_2$ have been determined quantum mechanically based on a new highly correlated five-dimensional potential energy surface. This surface was fitted to more than 60 000 ab initio points with a root mean square error of $\sim$ 1-2 cm$^{-1}$ in the region of the potential well whose depth is -374.0 cm$^{-1}$. The rotationally inelastic rate coefficients have been combined with spectroscopic data to generate the non-equilibrium spectrum of CH$_2$NH towards the interstellar molecular cloud Sgr~B2(N). The transition $1_{10}\to 1_{11}$ at 5.29 GHz is found to be inverted and the predicted weak maser emission spectrum is in excellent agreement with new observations performed with the 100-meter Green Bank Telescope.
△ Less
Submitted 5 June, 2018;
originally announced June 2018.
-
Modelling the molecular composition and nuclear-spin chemistry of collapsing prestellar sources
Authors:
P. Hily-Blant,
A. Faure,
C. Rist,
G. Pineau des Foreêts,
D. R. Flower
Abstract:
We study the gravitational collapse of prestellar sources and the associated evolution of their chemical composition. We use the University of Grenoble Alpes Astrochemical Network (UGAN), which includes reactions involving the different nuclear--spin states of H2, H3+, and of the hydrides of carbon, nitrogen, oxygen, and sulfur, for reactions involving up to seven protons. In addition, species-to-…
▽ More
We study the gravitational collapse of prestellar sources and the associated evolution of their chemical composition. We use the University of Grenoble Alpes Astrochemical Network (UGAN), which includes reactions involving the different nuclear--spin states of H2, H3+, and of the hydrides of carbon, nitrogen, oxygen, and sulfur, for reactions involving up to seven protons. In addition, species-to-species rate coefficients are provided for the ortho/para interconversion of the H3+ + H2 system and isotopic variants. The composition of the medium is followed from an initial steady state through the early phase of isothermal gravitational collapse. Both the freeze--out of the molecules on to grains and the coagulation of the grains were incorporated in the model. The predicted abundances and column densities of the spin isomers of ammonia and its deuterated forms are compared with those measured recently towards the prestellar cores H-MM1, L16293E, and Barnard B1. We find that gas--phase processes alone account satisfactorily for the observations, without recourse to grain-surface reactions. In particular, our model reproduces both the isotopologue abundance ratios and the ortho:para ratios of NH2D and NHD2 within observational uncertainties. More accurate observations are necessary to distinguish between full scrambling processes---as assumed in our gas-phase network---and direct nucleus- or atom-exchange reactions.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
Abundance of HCN and its C and N isotopologues in L1498
Authors:
Victor S. Magalhaes,
Pierre Hily-Blant,
Alexandre Faure,
Mario Hernandez-Vera,
Francois Lique
Abstract:
The isotopic ratio of nitrogen in nearby protoplanetary disks, recently measured in CN and HCN, indicates that a fractionated reservoir of nitrogen is available at the stage of comet formation. This reservoir presents a 3:1 enrichment in 15N relative to the elemental ratio of 330, identical to that between comets and the sun, suggesting that similar processes are responsible for the fractionation…
▽ More
The isotopic ratio of nitrogen in nearby protoplanetary disks, recently measured in CN and HCN, indicates that a fractionated reservoir of nitrogen is available at the stage of comet formation. This reservoir presents a 3:1 enrichment in 15N relative to the elemental ratio of 330, identical to that between comets and the sun, suggesting that similar processes are responsible for the fractionation in the protosolar nebula (PSN) and in PSN analogs. However, how the fractionation of nitrogen takes place is an open question. Previously obtained HCN/HC15N abundance ratios suggest that HCN may be enriched in 15N in prestellar cores, although doubts remain on these measurements, which rely on the double-isotopologue method. Here we present direct measurements of the HCN/H13CN and HCN/HC15N abundance ratios in the L1498 prestellar core based on spatially resolved spectra of HCN(1-0), (3-2), H13CN(1-0), and HC15N(1-0). We use state-of-the-art radiative transfer calculations using ALICO, a 1D radiative transfer code capable of treating hyperfine (hf) overlaps. From an analysis of dust emission maps of L1498, we derive a new physical structure. We also use new, high-accuracy HCN-H2 hf collisional rates, which enable us to quantitatively reproduce all the features of the HCN(1-0) and HCN(3-2) spectra, especially the anomalous hf line ratios. Special attention is devoted to derive meaningful uncertainties on the abundance ratios. The obtained values, HCN/H13CN=45(3) and HCN/HC15N=338(28), indicate that carbon is fractionated in HCN, but nitrogen is not. For the HCN/HC15N abundance ratio, our study validates to some extent analyses based on the single excitation temperature assumption. Comparisons with other measurements suggest significant core-to-core variability. Furthermore, the 13C enrichment we found in HCN could explain the superfractionation of nitrogen measured in solar system chondrites.
△ Less
Submitted 29 March, 2018;
originally announced March 2018.
-
Electron-impact excitation of diatomic hydride cations II: OH$^+$ and SH$^+$
Authors:
James R Hamilton,
Alexandre Faure,
Jonathan Tennyson
Abstract:
R-matrix calculations combined with the adiabatic-nuclei-rotation and Coulomb-Born approximations are used to compute electron-impact rotational rate coefficients for two open-shell diatomic cations of astrophysical interest: the hydoxyl and sulphanyl ions, OH$^+$ and SH$^+$. Hyperfine resolved rate coefficients are deduced using the infinite-order-sudden approximation. The propensity rule…
▽ More
R-matrix calculations combined with the adiabatic-nuclei-rotation and Coulomb-Born approximations are used to compute electron-impact rotational rate coefficients for two open-shell diatomic cations of astrophysical interest: the hydoxyl and sulphanyl ions, OH$^+$ and SH$^+$. Hyperfine resolved rate coefficients are deduced using the infinite-order-sudden approximation. The propensity rule $ΔF=Δj=ΔN=\pm 1$ is observed, as is expected for cations with a large dipole moment. A model for OH$^+$ excitation in the Orion Bar photon-dominated region (PDR) is presented which nicely reproduces Herschel observations for an electron fraction $x_e=10^{-4}$ and an OH$^+$ column density of $3\times 10^{13}$~cm$^{-2}$. Electron impact electronic excitation cross sections and rate coefficients for the ions are also presented.
△ Less
Submitted 21 February, 2018;
originally announced February 2018.
-
An efficient statistical method to compute molecular collisional rate coefficients
Authors:
Jérôme Loreau,
François Lique,
Alexandre Faure
Abstract:
Our knowledge about the "cold" Universe often relies on molecular spectra. A general property of such spectra is that the energy level populations are rarely at local thermodynamic equilibrium. Solving the radiative transfer thus requires the availability of collisional rate coefficients with the main colliding partners over the temperature range 10-1000 K. These rate coefficients are notoriously…
▽ More
Our knowledge about the "cold" Universe often relies on molecular spectra. A general property of such spectra is that the energy level populations are rarely at local thermodynamic equilibrium. Solving the radiative transfer thus requires the availability of collisional rate coefficients with the main colliding partners over the temperature range 10-1000 K. These rate coefficients are notoriously difficult to measure and expensive to compute. In particular, very few reliable collisional data exist for collisions involving reactive radicals or ions. Here we explore the use of a fast quantum statistical method to determine molecular collisional excitation rate coefficients. The method is benchmarked against accurate (but costly) close-coupling calculations. For collisions proceeding through the formation of a strongly-bound complex, the method is found to be highly satisfactory up to room temperature. Its accuracy decreases with the potential well depth and with increasing temperature, as expected. This new method opens the way to the determination of accurate inelastic collisional data involving key reactive species such as H3+, H2O+, and H3O+ for which exact quantum calculations are currently not feasible.
△ Less
Submitted 8 January, 2018;
originally announced January 2018.
-
Direct evidence of multiple reservoirs of volatile nitrogen in a protosolar nebula analogue
Authors:
Pierre Hily-Blant,
Victor Magalhães,
Joel Kastner,
Alexandre Faure,
Thierry Forveille,
Chunhua Qi
Abstract:
Isotopic ratios are keys to understanding the origin and early evolution of the solar system in the context of Galactic nucleosynthesis. The large range of measured $^{14}$N/$^{15}$N isotopic ratios in the solar system reflects distinct reservoirs of nitrogen whose origins remain to be determined. We have directly measured a C$^{14}$N/C$^{15}$N abundance ratio of 323$\pm$30 in the disk orbiting th…
▽ More
Isotopic ratios are keys to understanding the origin and early evolution of the solar system in the context of Galactic nucleosynthesis. The large range of measured $^{14}$N/$^{15}$N isotopic ratios in the solar system reflects distinct reservoirs of nitrogen whose origins remain to be determined. We have directly measured a C$^{14}$N/C$^{15}$N abundance ratio of 323$\pm$30 in the disk orbiting the nearby young star TW Hya. This value, which is in good agreement with nitrogen isotopic ratios measured for prestellar cores, likely reflects the primary present-day reservoir of nitrogen in the solar neighbourhood. These results support models invoking novae as primary $^{15}$N sources as well as outward migration of the Sun over its lifetime, and suggest that comets sampled a secondary, $^{15}$N-rich reservoir during solar system formation.
△ Less
Submitted 30 June, 2017;
originally announced June 2017.
-
NH_3(1_0-0_0) in the pre-stellar core L1544
Authors:
P. Caselli,
L. Bizzocchi,
E. Keto,
O. Sipila,
M. Tafalla,
L. Pagani,
L. E. Kristensen,
F. F. S. van der Tak,
C. M. Walmsley,
C. Codella,
B. Nisini,
Y. Aikawa,
A. Faure,
E. F. van Dishoeck
Abstract:
Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity…
▽ More
Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure.
Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations.
△ Less
Submitted 9 June, 2017;
originally announced June 2017.
-
Chemical complexity induced by efficient ice evaporation in the Barnard 5 molecular cloud
Authors:
Vianney Taquet,
Eva Wirström,
Steven B. Charnley,
Alexandre Faure,
Ana López-Sepulcre,
Carina M. Persson
Abstract:
Cold gas-phase water has recently been detected in a cold dark cloud, Barnard 5 located in the Perseus complex, by targeting methanol peaks as signposts for ice mantle evaporation. Observed morphology and abundances of methanol and water are consistent with a transient non-thermal evaporation process only affecting the outermost ice mantle layers, possibly triggering a more complex chemistry. We p…
▽ More
Cold gas-phase water has recently been detected in a cold dark cloud, Barnard 5 located in the Perseus complex, by targeting methanol peaks as signposts for ice mantle evaporation. Observed morphology and abundances of methanol and water are consistent with a transient non-thermal evaporation process only affecting the outermost ice mantle layers, possibly triggering a more complex chemistry. We present the detection of the Complex Organic Molecules (COMs) acetaldehyde and methyl formate as well as formic acid and ketene, and the tentative detection of di-methyl ether towards the methanol hotspot of Barnard 5 located between two dense cores using the single dish OSO 20m, IRAM 30m, and NRO 45m telescopes. The high energy cis- conformer of formic acid is detected, suggesting that formic acid is mostly formed at the surface of interstellar grains and then evaporated. The detection of multiple transitions for each species allows us to constrain their abundances through LTE and non-LTE methods. All the considered COMs show similar abundances between $\sim 1$ and $\sim 10$ % relative to methanol depending on the assumed excitation temperature. The non-detection of glycolaldehyde, an isomer of methyl formate, with a [glycolaldehyde]/[methyl formate] abundance ratio lower than 6 %, favours gas phase formation pathways triggered by methanol evaporation. According to their excitation temperatures derived in massive hot cores, formic acid, ketene, and acetaldehyde have been designated as "lukewarm" COMs whereas methyl formate and di-methyl ether were defined as "warm" species. Comparison with previous observations of other types of sources confirms that "lukewarm" and "warm" COMs show similar abundances in low-density cold gas whereas the "warm" COMs tend to be more abundant than the "lukewarm" species in warm protostellar cores.
△ Less
Submitted 5 June, 2017;
originally announced June 2017.
-
Collisional excitation of NH3 by atomic and molecular hydrogen
Authors:
N. Bouhafs,
C. Rist,
F. Daniel,
F. Dumouchel,
F. Lique,
L. Wiesenfeld,
A. Faure
Abstract:
We report extensive theoretical calculations on the rotation-inversion excitation of interstellar ammonia (NH3) due to collisions with atomic and molecular hydrogen (both para- and ortho-H2). Close-coupling calculations are performed for total energies in the range 1-2000 cm-1 and rotational cross sections are obtained for all transitions among the lowest 17 and 34 rotation-inversion levels of ort…
▽ More
We report extensive theoretical calculations on the rotation-inversion excitation of interstellar ammonia (NH3) due to collisions with atomic and molecular hydrogen (both para- and ortho-H2). Close-coupling calculations are performed for total energies in the range 1-2000 cm-1 and rotational cross sections are obtained for all transitions among the lowest 17 and 34 rotation-inversion levels of ortho- and para-NH3, respectively. Rate coefficients are deduced for kinetic temperatures up to 200 K. Propensity rules for the three colliding partners are discussed and we also compare the new results to previous calculations for the spherically symmetrical He and para-H2 projectiles. Significant differences are found between the different sets of calculations. Finally, we test the impact of the new rate coefficients on the calibration of the ammonia thermometer. We find that the calibration curve is only weakly sensitive to the colliding partner and we confirm that the ammonia thermometer is robust.
△ Less
Submitted 28 November, 2018; v1 submitted 1 June, 2017;
originally announced June 2017.
-
SOFIA/GREAT Discovery of Terahertz Water Masers
Authors:
David A. Neufeld,
Gary J. Melnick,
Michael J. Kaufman,
Helmut Wiesemeyer,
Rolf Güsten,
Alex Kraus,
Karl M. Menten,
Oliver Ricken,
Alexandre Faure
Abstract:
We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser…
▽ More
We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser transition were carried out towards all four sources using the Effelsberg 100m telescope. The measured line fluxes imply 22 GHz / 1.296 THz photon luminosity ratios of 0.012, 0.12, and 0.83 respectively for W Hya, U Her, and VY CMa, values that confirm the 22 GHz maser transition to be unsaturated in W Hya and U Her. We also detected the 1.884888 THz 8(45)-7(53) transition toward W Hya and VY CMa, and the 1.278266 THz 7(43)-6(52) transition toward VY CMa. Like the 22 GHz maser transition, all three of the THz emission lines detected here originate from the ortho-H2O spin isomer. Based upon a model for the circumstellar envelope of W Hya, we estimate that stimulated emission is responsible for ~ 85% of the observed 1.296 THz line emission, and thus that this transition may be properly described as a terahertz-frequency maser. In the case of the 1.885 THz transition, by contrast, our W Hya model indicates that the observed emission is dominated by spontaneous radiative decay, even though a population inversion exists.
△ Less
Submitted 26 May, 2017;
originally announced May 2017.
-
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.
-
History of the solar-type protostar IRAS16293-2422 as told by the cyanopolyynes
Authors:
A. A. Jaber,
C. Ceccarelli,
C. Kahane,
S. Viti,
N. Balucani,
E. Caux,
A. Faure,
B. Lefloch,
F. Lique,
E. Mendoza,
D. Quenard,
L. Wiesenfeld
Abstract:
Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type pro…
▽ More
Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K.
△ Less
Submitted 4 October, 2016;
originally announced October 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.
-
Collisional excitation of HC3N by para- and ortho-H2
Authors:
A. Faure,
F. Lique,
L. Wiesenfeld
Abstract:
New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-chann…
▽ More
New calculations for rotational excitation of cyanoacetylene by collisions with hydrogen molecules are performed to include the lowest 38 rotational levels of HC3N and kinetic temperatures to 300 K. Calculations are based on the interaction potential of Wernli et al. A&A, 464, 1147 (2007) whose accuracy is checked against spectroscopic measurements of the HC3N-H2 complex. The quantum coupled-channel approach is employed and complemented by quasi-classical trajectory calculations. Rate coefficients for ortho-H2 are provided for the first time. Hyperfine resolved rate coefficients are also deduced. Collisional propensity rules are discussed and comparisons between quantum and classical rate coefficients are presented. This collisional data should prove useful in interpreting HC3N observations in the cold and warm ISM, as well as in protoplanetary disks.
△ Less
Submitted 12 May, 2016;
originally announced May 2016.
-
Deuteration of ammonia in the starless core Ophiuchus/H-MM1
Authors:
Jorma Harju,
Fabien Daniel,
Olli Sipilä,
Paola Caselli,
Jaime E. Pineda,
Rachel K. Friesen,
Anna Punanova,
Rolf Güsten,
Laurent Wiesenfeld,
Philip C. Myers,
Alexandre Faure,
Pierre Hily-Blant,
Claire Rist,
Erik Rosolowsky,
Stephan Schlemmer,
Yancy L. Shirley
Abstract:
Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain…
▽ More
Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.
△ Less
Submitted 22 December, 2016; v1 submitted 19 April, 2016;
originally announced April 2016.
-
N$_2$H$^+$ and N$^{15}$NH$^+$ towards the prestellar core 16293E in L1689N
Authors:
F. Daniel,
A. Faure,
L. Pagani,
F. Lique,
M. Gérin,
D. C. Lis,
P. Hily-Blant,
A. Bacmann,
E. Roueff
Abstract:
Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statist…
▽ More
Understanding the processes that could lead to enrichment of molecules in $^{15}$N atoms is of particular interest in order to shed light on the relatively large variations observed in the $^{14}$N/$^{15}$N ratio in various solar system environments. Currently, the sample of molecular clouds where $^{14}$N/$^{15}$N ratios have been measured is small and has to be enlarged in order to allow statistically significant studies. In particular, the N$_2$H$^+$ molecule currently shows the largest spread of $^{14}$N/$^{15}$N ratios in high-mass star forming regions. However, the $^{14}$N/$^{15}$N ratio in N$_2$H$^+$ was obtained in only two low-mass star forming regions (L1544 and B1b). The current work extends this sample to a third dark cloud. We targeted the 16293E prestellar core, where the N$^{15}$NH$^+$ $J$=1-0 line was detected. Using a model previously developed for the physical structure of the source, we solved the molecular excitation with a non-local radiative transfer code. For that purpose, we computed specific collisional rate coefficients for the N$^{15}$NH$^+$-H$_2$ collisional system. As a first step of the analysis, the N$_2$H$^+$ abundance profile was constrained by reproducing the N$_2$H$^+$ $J$=1-0 and 3-2 maps. A scaling factor was then applied to this profile to match the N$^{15}$NH$^+$ $J$=1-0 spectrum. We derive a column density ratio N$_2$H$^+$ / N$^{15}$NH$^+$ = $330^{+170}_{-100}$. The current estimate $\sim$330 agrees with the value typical of the elemental isotopic ratio in the local ISM. It is however lower than in some other cores, where values as high as 1300 have been reported.
△ Less
Submitted 23 March, 2016;
originally announced March 2016.
-
Detection of protonated formaldehyde in the prestellar core L1689B
Authors:
A. Bacmann,
E. García-García,
A. Faure
Abstract:
Complex organic molecules (COMs) are detected in many regions of the interstellar medium, including prestellar cores. However, their formation mechanisms in cold (~10 K) cores remain to this date poorly understood. The formyl radical HCO is an important candidate precursor for several O-bearing terrestrial COMs in cores, as an abundant building block of many of these molecules. Several chemical ro…
▽ More
Complex organic molecules (COMs) are detected in many regions of the interstellar medium, including prestellar cores. However, their formation mechanisms in cold (~10 K) cores remain to this date poorly understood. The formyl radical HCO is an important candidate precursor for several O-bearing terrestrial COMs in cores, as an abundant building block of many of these molecules. Several chemical routes have been proposed to account for its formation, both on grain surfaces, as an incompletely hydrogenated product of H addition to frozen-out CO molecules, or in the gas phase, either the product of the reaction between H2CO and a radical, or as a product of dissociative recombination of protonated formaldehyde H2COH+. The detection and abundance determination of H2COH+, if present, could provide clues as to whether this latter scenario might apply. We searched for protonated formaldehyde H2COH+ in the prestellar core L1689B using the IRAM 30m telescope. The H2COH+ ion is unambiguously detected, for the first time in a cold (~10 K) source. The derived abundance agrees with a scenario in which the formation of H2COH+ results from the protonation of formaldehyde. We use this abundance value to constrain the branching ratio of the dissociative recombination of H2COH+ towards the HCO channel to ~10-30%. This value could however be smaller if HCO can be efficiently formed from gas-phase neutral-neutral reactions, and we stress the need for laboratory measurements of the rate constants of these reactions at 10 K. Given the experimental difficulties in measuring branching ratios experimentally, observations can bring valuable constraints on these values, and provide a useful input for chemical networks.
△ Less
Submitted 1 March, 2016;
originally announced March 2016.
-
Collisional excitation of doubly and triply deuterated ammonia ND$_2$H and ND$_3$ by H$_2$
Authors:
F. Daniel,
C. Rist,
A. Faure,
E. Roueff,
M. Gérin,
D. C. Lis,
P. Hily-Blant,
A. Bacmann,
L. Wiesenfeld
Abstract:
The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H…
▽ More
The availability of collisional rate coefficients is a prerequisite for an accurate interpretation of astrophysical observations, since the observed media often harbour densities where molecules are populated under non--LTE conditions. In the current study, we present calculations of rate coefficients suitable to describe the various spin isomers of multiply deuterated ammonia, namely the ND$_2$H and ND$_3$ isotopologues. These calculations are based on the most accurate NH$_3$--H$_2$ potential energy surface available, which has been modified to describe the geometrical changes induced by the nuclear substitutions. The dynamical calculations are performed within the close--coupling formalism and are carried out in order to provide rate coefficients up to a temperature of $T$ = 50K. For the various isotopologues/symmetries, we provide rate coefficients for the energy levels below $\sim$ 100 cm$^{-1}$. Subsequently, these new rate coefficients are used in astrophysical models aimed at reproducing the NH$_2$D, ND$_2$H and ND$_3$ observations previously reported towards the prestellar cores B1b and 16293E. We thus update the estimates of the corresponding column densities and find a reasonable agreement with the previous models. In particular, the ortho--to--para ratios of NH$_2$D and NHD$_2$ are found to be consistent with the statistical ratios.
△ Less
Submitted 11 January, 2016;
originally announced January 2016.
-
The origin of gas-phase HCO and CH3O radicals in prestellar cores
Authors:
Aurore Bacmann,
Alexandre Faure
Abstract:
The recent unexpected detection of terrestrial complex organic molecules in the cold (~ 10 K) gas has cast doubts on the commonly accepted formation mechanisms of these species. Standard gas-phase mechanisms are inefficient and tend to underproduce these molecules, and many of the key reactions involved are unconstrained. Grain-surface mechanisms, which were presented as a viable alternative, suff…
▽ More
The recent unexpected detection of terrestrial complex organic molecules in the cold (~ 10 K) gas has cast doubts on the commonly accepted formation mechanisms of these species. Standard gas-phase mechanisms are inefficient and tend to underproduce these molecules, and many of the key reactions involved are unconstrained. Grain-surface mechanisms, which were presented as a viable alternative, suffer from the fact that they rely on grain surface diffusion of heavy radicals, which is not possible thermally at very low temperatures. One of the simplest terrestrial complex organic molecules, methanol is believed to form on cold grain surfaces following from successive H atom additions on CO. Unlike heavier species, H atoms are very mobile on grain surfaces even at 10 K. Intermediate species involved in grain surface methanol formation by CO hydrogenation are the radicals HCO and CH3O, as well as the stable species formaldehyde H2CO. These radicals are thought to be precursors of complex organic molecules on grain surfaces. We present new observations of the HCO and CH3O radicals in a sample of prestellar cores and carry out an analysis of the abundances of the species HCO, H2CO, CH3O, and CH3OH, which represent the various stages of grain-surface hydrogenation of CO to CH3OH. The abundance ratios between the various intermediate species in the hydrogenation reaction of CO on grains are similar in all sources of our sample, HCO:H2CO:CH3O:CH3OH ~ 10:100:1:100. We argue that these ratios may not be representative of the primordial abundances on the grains but, rather, suggest that the radicals HCO and CH3O are gas-phase products of the precursors H2CO and CH3OH, respectively. Gas-phase pathways are considered and simple estimates of HCO and CH3O abundances are compared to the observations. Critical reaction rate constants, branching ratios, and intermediate species are finally identified.
△ Less
Submitted 7 January, 2016;
originally announced January 2016.
-
The NH$_2$D hyperfine structure revealed by astrophysical observations
Authors:
F. Daniel,
L. H. Coudert,
A. Punanova,
J. Harju,
A. Faure,
E. Roueff,
O. Sipilä,
P. Caselli,
R. Güsten,
A. Pon,
J. E. Pineda
Abstract:
The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density est…
▽ More
The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.
△ Less
Submitted 2 January, 2016;
originally announced January 2016.
-
Ortho-to-para ratio of NH2. Herschel-HIFI observations of ortho- and para-NH2 rotational transitions towards W31C, W49N, W51 and G34.3+0.1
Authors:
C. M. Persson,
A. O. H. Olofsson,
R. Le Gal,
E. S. Wirström,
G. E. Hassel,
E. Herbst,
M. Olberg,
A. Faure,
P. Hily-Blant,
J. H. Black,
M. Gerin,
D. Lis,
F. Wyrowski
Abstract:
We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores…
▽ More
We have used the Herschel-HIFI instrument to observe both nuclear spin symmetries of amidogen (NH2) towards the high-mass star-forming regions W31C (G10.6-0.4), W49N (G43.2-0.1), W51 (G49.5-0.4) and G34.3+0.1. The aim is to investigate the ratio of nuclear spin types, the ortho-to-para ratio (OPR), of NH2. The excited NH2 transitions are used to construct radiative transfer models of the hot cores and surrounding envelopes in order to investigate the excitation and possible emission of the ground state rotational transitions of ortho-NH2 N_(K_a,K_c} J=1_(1,1) 3/2 - 0_(0,0) 1/2 and para-NH2 2_(1,2) 5/2 - 1_(0,1) 3/2$ used in the OPR calculations. Our best estimate of the average OPR in the envelopes lie above the high temperature limit of three for W49N, specifically 3.5 with formal errors of \pm0.1, but for W31C, W51, and G34.3+0.1 we find lower values of 2.5\pm0.1, 2.7\pm0.1, and 2.3\pm0.1, respectively. Such low values are strictly forbidden in thermodynamical equilibrium since the OPR is expected to increase above three at low temperatures. In the translucent interstellar gas towards W31C, where the excitation effects are low, we find similar values between 2.2\pm0.2 and 2.9\pm0.2. In contrast, we find an OPR of 3.4\pm0.1 in the dense and cold filament connected to W51, and also two lower limits of >4.2 and >5.0 in two other translucent gas components towards W31C and W49N. At low temperatures (T \lesssim 50 K) the OPR of H2 is <10^-1, far lower than the terrestrial laboratory normal value of three. In such a "para-enriched H2" gas, our astrochemical models can reproduce the variations of the observed OPR, both below and above the thermodynamical equilibrium value, by considering nuclear-spin gas-phase chemistry. The models suggest that values below three arise in regions with temperatures >20-25 K, depending on time, and values above three at lower temperatures.
△ Less
Submitted 14 December, 2015; v1 submitted 17 November, 2015;
originally announced November 2015.
-
Hydrogen isotope exchanges between water and methanol in interstellar ices
Authors:
A. Faure,
M. Faure,
P. Theulé,
E. Quirico,
B. Schmitt
Abstract:
The deuterium fractionation of gas-phase molecules in hot cores is believed to reflect the composition of interstellar ices. The deuteration of methanol is a major puzzle, however, because the isotopologue ratio [CH2DOH]/[CH3OD], which is predicted to be equal to 3 by standard grain chemistry models, is much larger (~20) in low-mass hot corinos and significantly lower (~1) in high-mass hot cores.…
▽ More
The deuterium fractionation of gas-phase molecules in hot cores is believed to reflect the composition of interstellar ices. The deuteration of methanol is a major puzzle, however, because the isotopologue ratio [CH2DOH]/[CH3OD], which is predicted to be equal to 3 by standard grain chemistry models, is much larger (~20) in low-mass hot corinos and significantly lower (~1) in high-mass hot cores. This dichotomy in methanol deuteration between low-mass and massive protostars is currently not understood. In this study, we report a simplified rate equation model of the deuterium chemistry occurring in the icy mantles of interstellar grains. We apply this model to the chemistry of hot corinos and hot cores, with IRAS 16293-2422 and the Orion~KL Compact Ridge as prototypes, respectively. The chemistry is based on a statistical initial deuteration at low temperature followed by a warm-up phase during which thermal hydrogen/deuterium (H/D) exchanges occur between water and methanol. The exchange kinetics is incorporated using laboratory data. The [CH2DOH]/[CH3OD] ratio is found to scale inversely with the D/H ratio of water, owing to the H/D exchange equilibrium between the hydroxyl (-OH) functional groups of methanol and water. Our model is able to reproduce the observed [CH2DOH]/[CH3OD] ratios provided that the primitive fractionation of water ice [HDO]/[H2O] is ~ 2% in IRAS 16293-2422 and ~0.6% in Orion~KL. We conclude that the molecular D/H ratios measured in hot cores may not be representative of the original mantles because molecules with exchangeable deuterium atoms can equilibrate with water ice during the warm-up phase.
△ Less
Submitted 3 November, 2015;
originally announced November 2015.
-
Electron-impact excitation of diatomic hydride cations I: HeH$^+$, CH$^+$, ArH$^+$
Authors:
James R. Hamilton,
Alexandre Faure,
Jonathan Tennyson
Abstract:
{\bf R}-matrix calculations combined with the adiabatic nuclei approximation are used to compute electron-impact rotiational excitation rates for three closed-shell diatomic cations, HeH$^+$, CH$^+$, ArH$^+$. Comparisons with previous studies show that an improved treatment of threshold effects leads to significant changes in the low temperature rates, furthermore the new calculations suggest that…
▽ More
{\bf R}-matrix calculations combined with the adiabatic nuclei approximation are used to compute electron-impact rotiational excitation rates for three closed-shell diatomic cations, HeH$^+$, CH$^+$, ArH$^+$. Comparisons with previous studies show that an improved treatment of threshold effects leads to significant changes in the low temperature rates, furthermore the new calculations suggest that excitation of CH$^+$ is dominated by $ΔJ =1$ transitions as is expected for cations with a large dipole moment. A model for ArH$^+$ excitation in the Crab Nebula is presented which gives results consistent with the observations for electron densities in the range $2-3\times 10^3$~cm$^{-3}$.
△ Less
Submitted 17 October, 2015;
originally announced October 2015.