Showing 1–45 of 45 results for author: Lique, F

Rovibrational (de-)excitation of H$_{2}$ by He revisited

Authors: Hubert Jóźwiak, Franck Thibault, Alexandra Viel, Piotr Wcisło, François Lique

Abstract: Collisional (de-)excitation of H$_{2}$ by helium plays an important role in the thermal balance and chemistry of various astrophysical environments, making accurate rate coefficients essential for the interpretation of observations of the interstellar medium. Our goal is to utilize a state-of-the-art potential energy surface (PES) to provide comprehensive state-to-state rate coefficients for He-in… ▽ More Collisional (de-)excitation of H$_{2}$ by helium plays an important role in the thermal balance and chemistry of various astrophysical environments, making accurate rate coefficients essential for the interpretation of observations of the interstellar medium. Our goal is to utilize a state-of-the-art potential energy surface (PES) to provide comprehensive state-to-state rate coefficients for He-induced transitions among rovibrational levels of H$_{2}$. We perform quantum scattering calculations for the H$_{2}$-He system and provide state-to-state rate coefficients for 1 089 transitions between rovibrational levels of H$_{2}$ with internal energies up to 15 000 cm$^{-1}$ for temperatures ranging from 20 to 8 000 K. Our results show good agreement with previous calculations for pure rotational transitions between low-lying rotational levels, but we find significant discrepancies for rovibrational processes involving highly-excited rotational and vibrational states. We attribute these differences to two key factors: the broader range of intramolecular distances covered by ab initio points, and the superior accuracy of the PES, resulting from the utilization of the state-of-the-art quantum chemistry methods, compared to the previous lower-level calculations. Radiative transfer calculations performed with the new collisional data indicate that the population of rotational levels in excited vibrational states experiences significant modifications, highlighting the critical need for this updated dataset in models of high-temperature astrophysical environments. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Journal ref: A&A 685, A113 (2024)

Collisional excitation of PO$^+$ by $para$-H$_2$

Authors: F. Tonolo, L. Bizzocchi, V. M. Rivilla, F. Lique, M. Melosso, C. Puzzarini

Abstract: We report the derivation of rate coefficients for the rotational (de-)excitation of PO$^+$ induced by collisions with H$_2$. The calculations were performed on a four-dimensional potential energy surface, obtained on top of highly accurate $ab$ $initio$ energy points. Preliminary tests pointed out the low influence of the coupling between $j=0$ and the higher rotational levels of H$_2$ on the cros… ▽ More We report the derivation of rate coefficients for the rotational (de-)excitation of PO$^+$ induced by collisions with H$_2$. The calculations were performed on a four-dimensional potential energy surface, obtained on top of highly accurate $ab$ $initio$ energy points. Preliminary tests pointed out the low influence of the coupling between $j=0$ and the higher rotational levels of H$_2$ on the cross sections values, thus allowing to neglect the rotational structure of H$_2$. On this basis, state-to-state collisional rate coefficients were derived for temperatures ranging from 5 to 200 K. Radiative transfer calculations have been used to model the recent observation of PO$^+$ in the G+0.693-0.027 molecular cloud, in order to evaluate the possible impact of non-LTE models on the determination of its physical conditions. The derived column density was found to be approximately $\sim 3.7\times10^{11}$ cm$^{-2}$, which is 60\% (a factor of $\sim 1.7$) smaller than the previously LTE-derived value. Extensive simulations show that PO$^+$ low-$j$ rotational lines exhibit maser behavior at densities between $10^4$ and $10^6$ cm$^{-3}$, thus highlighting the importance of a proper treatment of the molecular collisions to accurately model PO$^+$ emissions in the interstellar medium. △ Less

Submitted 12 October, 2023; originally announced October 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.

Similar levels of deuteration in the pre-stellar core L1544 and the protostellar core HH211

Authors: K. Giers, S. Spezzano, P. Caselli, E. Wirström, O. Sipilä, J. E. Pineda, E. Redaelli, C. T. Bop, F. Lique

Abstract: In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostella… ▽ More In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostellar core HH211 in Perseus. We used single-dish observations of CCH, HCN, HNC, HCO$^+$, and their $^{13}$C-, $^{18}$O- and D-bearing isotopologues, detected with the Onsala 20m telescope. We derived the column densities and the deuterium fractions of the molecules. Additionally, we used radiative transfer simulations and results from chemical modelling to reproduce the observed molecular lines. We used new collisional rate coefficients for HNC, HN$^{13}$C, DNC, and DCN that consider the hyperfine structure of these molecules. We find high levels of deuteration for CCH (10%) in both sources, consistent with other carbon chains, and moderate levels for HCN (5-7%) and HNC (8%). The deuterium fraction of HCO$^+$ is enhanced towards HH211, most likely caused by isotope-selective photodissociation of C$^{18}$O. Similar levels of deuteration show that the process is likely equally efficient towards both cores, suggesting that the protostellar envelope still retains the chemical composition of the original pre-stellar core. The fact that the two cores are embedded in different molecular clouds also suggests that environmental conditions do not have a significant effect on the deuteration within dense cores. Radiative transfer modelling shows that it is necessary to include the outer layers of the cores to consider the effects of extended structures. Besides HCO$^+$ observations, HCN observations towards L1544 also require the presence of an outer diffuse layer where the molecules are relatively abundant. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 27 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 676, A78 (2023)

Cyanopolyyne chemistry in the L1544 prestellar core: new insights from GBT observations

Authors: Eleonora Bianchi, Anthony Remijan, Claudio Codella, Cecilia Ceccarelli, Francois Lique, Silvia Spezzano, Nadia Balucani, Paola Caselli, Eric Herbst, Linda Podio, Charlotte Vastel, Brett McGuire

Abstract: We report a comprehensive study of the cyanopolyyne chemistry in the prototypical prestellar core L1544. Using the 100m Robert C. Byrd Green Bank Telescope (GBT) we observe 3 emission lines of HC$_3$N, 9 lines of HC$_5$N, 5 lines of HC$_7$N, and 9 lines of HC$_9$N. HC$_9$N is detected for the first time towards the source. The high spectral resolution ($\sim$ 0.05 km s$^{-1}$) reveals double-peak… ▽ More We report a comprehensive study of the cyanopolyyne chemistry in the prototypical prestellar core L1544. Using the 100m Robert C. Byrd Green Bank Telescope (GBT) we observe 3 emission lines of HC$_3$N, 9 lines of HC$_5$N, 5 lines of HC$_7$N, and 9 lines of HC$_9$N. HC$_9$N is detected for the first time towards the source. The high spectral resolution ($\sim$ 0.05 km s$^{-1}$) reveals double-peak spectral line profiles with the redshifted peak a factor 3-5 brighter. Resolved maps of the core in other molecular tracers indicates that the southern region is redshifted. Therefore, the bulk of the cyanopolyyne emission is likely associated with the southern region of the core, where free carbon atoms are available to form long chains, thanks to the more efficient illumination of the interstellar field radiation. We perform a simultaneous modelling of the HC$_5$N, HC$_7$N, and HC$_9$N lines, to investigate the origin of the emission. To enable this analysis, we performed new calculation of the collisional coefficients. The simultaneous fitting indicates a gas kinetic temperature of 5--12 K, a source size of 80$\arcsec$, and a gas density larger than 100 cm$^{-3}$. The HC$_5$N:HC$_7$N:HC$_9$N abundance ratios measured in L1544 are about 1:6:4. We compare our observations with those towards the the well-studied starless core TMC-1 and with the available measurements in different star-forming regions. The comparison suggests that a complex carbon chain chemistry is active in other sources and it is related to the presence of free gaseous carbon. Finally, we discuss the possible formation and destruction routes in the light of the new observations. △ Less

Submitted 2 February, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: 16 pages, 7 figure, 3 tables

Linking the dust and chemical evolution: Taurus and Perseus -- New collisional rates for HCN, HNC, and their C, N, and H isotopologues

Authors: D. Navarro-Almaida, C. T. Bop, F. Lique, G. Esplugues, M. Rodríguez-Baras, C. Kramer, C. E. Romero, A. Fuente, P. Caselli, P. Riviére-Marichalar, J. M. Kirk, A. Chacón-Tanarro, E. Roueff, T. Mroczkowski, T. Bhandarkar, M. Devlin, S. Dicker, I. Lowe, B. Mason, C. L. Sarazin, J. Sievers

Abstract: HCN, HNC, and their isotopologues are ubiquitous molecules that can serve as chemical thermometers and evolutionary tracers to characterize star-forming regions. Despite their importance in carrying information that is vital to studies of the chemistry and evolution of star-forming regions, the collision rates of some of these molecules have not been available for rigorous studies in the past. We… ▽ More HCN, HNC, and their isotopologues are ubiquitous molecules that can serve as chemical thermometers and evolutionary tracers to characterize star-forming regions. Despite their importance in carrying information that is vital to studies of the chemistry and evolution of star-forming regions, the collision rates of some of these molecules have not been available for rigorous studies in the past. We perform an up-to-date gas and dust chemical characterization of two different star-forming regions, TMC 1-C and NGC 1333-C7, using new collisional rates of HCN, HNC, and their isotopologues. We investigated the possible effects of the environment and stellar feedback in their chemistry and their evolution. With millimeter observations, we derived their column densities, the C and N isotopic fractions, the isomeric ratios, and the deuterium fractionation. The continuum data at 3 mm and 850 $μ$m allowed us to compute the emissivity spectral index and look for grain growth as an evolutionary tracer. The H$^{13}$CN/HN$^{13}$C ratio is anticorrelated with the deuterium fraction of HCN, thus it can readily serve as a proxy for the temperature. The spectral index $(β\sim 1.34-2.09)$ shows a tentative anticorrelation with the H$^{13}$CN/HN$^{13}$C ratio, suggesting grain growth in the evolved, hotter, and less deuterated sources. Unlike TMC 1-C, the south-to-north gradient in dust temperature and spectral index observed in NGC 1333-C7 suggests feedback from the main NGC 1333 cloud. With this up-to-date characterization of two star-forming regions, we found that the chemistry and the physical properties are tightly related. The dust temperature, deuterium fraction, and the spectral index are complementary evolutionary tracers. The large-scale environmental factors may dominate the chemistry and evolution in clustered star-forming regions. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 25 pages, 20 figures

Journal ref: A&A 670, A110 (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.

Comments: 13 pages, 7 figures

MSC Class: 85-10; 85-08

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.

Comments: 8 pages

Journal ref: MNRAS 516 5964 (2022)

Tracing the contraction of the pre-stellar core L1544 with HC$^{17}$O$^+$ $J$ = 1-0 emission

Authors: J. Ferrer Asensio, S. Spezzano, P. Caselli, F. O. Alves, O. Sipilä, E. Redaelli, L. Bizzocchi, F. Lique, A. Mullins

Abstract: Spectral line profiles of several molecules observed towards the pre-stellar core L1544 appear double-peaked. For abundant molecular species this line morphology has been linked to self-absorption. However, the physical process behind the double-peaked morphology for less abundant species is still under debate. In order to understand the cause behind the double-peaked spectra of optically thin tra… ▽ More Spectral line profiles of several molecules observed towards the pre-stellar core L1544 appear double-peaked. For abundant molecular species this line morphology has been linked to self-absorption. However, the physical process behind the double-peaked morphology for less abundant species is still under debate. In order to understand the cause behind the double-peaked spectra of optically thin transitions and their link to the physical structure of pre-stellar cores, we present high-sensitivity and high-spectral resolution HC$^{17}$O$^+$ $J =$1-0 observations towards the dust peak in L1544. We observed the HC$^{17}$O$^+$ (1-0) spectrum with the Institut de Radioastronomie Millimétrique (IRAM) 30m telescope. By using new state-of-the-art collisional rate coefficients, a physical model for the core and the fractional abundance profile of HC$^{17}$O$^+$, the hyperfine structure of this molecular ion is modelled for the first time with the radiative transfer code LOC applied to the predicted chemical structure of a contracting pre-stellar core. We applied the same analysis to the chemically related C$^{17}$O molecule. The observed HC$^{17}$O$^+$(1-0) and C$^{17}$O(1-0) lines have been successfully reproduced with a non-local thermal equilibrium (LTE) radiative transfer model applied to chemical model predictions for a contracting pre-stellar core. An upscaled velocity profile (by 30%) is needed to reproduce the HC$^{17}$O$^+$(1-0) observations. The double peaks observed in the HC$^{17}$O$^+$(1-0) hyperfine components are due to the contraction motions at densities close to the critical density of the transition ($\sim$10$^{5}$ cm$^{-3}$) and to the fact that the HCO$^{+}$ fractional abundance decreases toward the centre. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Hyperfine resolved rate coefficients of HC17O+ with H2 (j = 0)

Authors: Francesca Tonolo, François Lique, Mattia Melosso, Cristina Puzzarini, Luca Bizzocchi

Abstract: The formyl cation (HCO+) is one of the most abundant ions in molecular clouds and plays a major role in the interstellar chemistry. For this reason, accurate collisional rate coefficients for the rotational excitation of HCO+ and its isotopes due to the most abundant perturbing species in interstellar environments are crucial for non-local thermal equilibrium models and deserve special attention.… ▽ More The formyl cation (HCO+) is one of the most abundant ions in molecular clouds and plays a major role in the interstellar chemistry. For this reason, accurate collisional rate coefficients for the rotational excitation of HCO+ and its isotopes due to the most abundant perturbing species in interstellar environments are crucial for non-local thermal equilibrium models and deserve special attention. In this work, we determined the first hyperfine resolved rate coefficients of HC17O+ in collision with H2 (j=0). Indeed, despite no scattering calculations on its collisional parameters have been performed so far, the HC17O+ isotope assumes a prominent role for astrophysical modelling applications. Computations are based on a new four dimensional (4D) potential energy surface, obtained at the CCSD(T)-F12a/aug-cc-pVQZ level of theory. A test on the corresponding cross section values pointed out that, to a good approximation, the influence of the coupling between rotational levels of H2 can be ignored. For this reason, the H2 collider has been treated as a spherical body and an average of the potential based on five orientations of H2 has been employed for scattering calculations. State-to-state rate coefficients resolved for the HC17O+ hyperfine structure for temperature ranging from 5 to 100 K have been computed using recoupling techniques. This study provides the first determination of HC17O+/H2 inelastic rate coefficients directly computed from full quantum close-coupling equations, thus supporting the reliability of future radiative transfer modellings of HC17O+ in interstellar environments. △ Less

Submitted 21 August, 2022; originally announced August 2022.

Comments: 9 pages, 6 figures

A SUBLIME 3D Model for Cometary Coma Emission: the Hypervolatile-Rich Comet C/2016 R2 (PanSTARRS)

Authors: M. A. Cordiner, I. M. Coulson, E. Garcia-Berrios, C. Qi, F. Lique, M. Zoltowski, M. de Val-Borro, Y. -J. Kuan, W. -H. Ip, S. Mairs, N. X. Roth, S. B. Charnley, S. N. Milam, W. -L Tseng, Y. -L Chuang

Abstract: The coma of comet C/2016 R2 (PanSTARRS) is one of the most chemically peculiar ever observed, in particular due to its extremely high CO/H2O and N2+/H2O ratios}, and unusual trace volatile abundances. However, the complex shape of its CO emission lines, as well as uncertainties in the coma structure and excitation, has lead to ambiguities in the total CO production rate. We performed high resoluti… ▽ More The coma of comet C/2016 R2 (PanSTARRS) is one of the most chemically peculiar ever observed, in particular due to its extremely high CO/H2O and N2+/H2O ratios}, and unusual trace volatile abundances. However, the complex shape of its CO emission lines, as well as uncertainties in the coma structure and excitation, has lead to ambiguities in the total CO production rate. We performed high resolution, spatially, spectrally and temporally resolved CO observations using the James Clerk Maxwell Telescope (JCMT) and Submillimeter Array (SMA) to elucidate the outgassing behaviour of C/2016 R2. Results are analyzed using a new, time-dependent, three dimensional radiative transfer code (SUBLIME), incorporating for the first time, accurate state-to-state collisional rate coefficients for the CO--CO system. The total CO production rate was found to be in the range $(3.8-7.6)\times10^{28}$ s$^{-1}$ between 2018-01-13 and 2018-02-01, with a mean value of $(5.3\pm0.6)\times10^{28}$ s$^{-1}$ at r_H = 2.8-2.9 au. The emission is concentrated in a near-sunward jet, with an outflow velocity $0.51\pm0.01$ km/s, compared to $0.25\pm0.01$ km/s in the ambient (and night-side) coma. Evidence was also found for an extended source of CO emission, possibly due to icy grain sublimation around $1.2\times10^5$ km from the nucleus. Based on the coma molecular abundances, we propose that the nucleus ices of C/2016 R2 can be divided into a rapidly sublimating apolar phase, rich in CO, CO2, N2 and CH3OH, and a predominantly frozen (or less abundant), polar phase containing more H2O, CH4, H2CO and HCN. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: Accepted for publication in ApJ

An improved study of HCO+ and He system: interaction potential, collisional relaxation and pressure broadening

Authors: F. Tonolo, L. Bizzocchi, M. Melosso, F. Lique, L. Dore, V. Barone, C. Puzzarini

Abstract: In light of its ubiquitous presence in the interstellar gas, the chemistry and reactivity of the HCO+ ion requires special attention. The availability of up-to-date collisional data between this ion and the most abundant perturbing species in the interstellar medium is a critical resource in order to derive reliable values of its molecular abundance from astronomical observations. This work intend… ▽ More In light of its ubiquitous presence in the interstellar gas, the chemistry and reactivity of the HCO+ ion requires special attention. The availability of up-to-date collisional data between this ion and the most abundant perturbing species in the interstellar medium is a critical resource in order to derive reliable values of its molecular abundance from astronomical observations. This work intends to provide improved scattering parameters for the HCO+ and He collisional system. We have tested the accuracy of explicitly correlated coupled-cluster methods for mapping the short- and long-range multi-dimensional potential energy surface of atom-ion systems. A validation of the methodology employed for the calculation of the potential well has been obtained from the comparison with experimentally derived bound-state spectroscopic parameters. Finally, by solving the close-coupling scattering equations, we have derived the pressure broadening and shift coefficients for the first six rotational transitions of HCO+ as well as inelastic state-to-state transition rates up to j = 5 in the 5-100 K temperature interval. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 8 pages, 5 figures

Journal ref: J. Chem. Phys. 155, 234306 (2021)

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.

Comments: 9 pages, 5 figures, accepted for publication in ApJ (2021/12/14)

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.

Comments: Accepted for publication in A&A

Journal ref: A&A 658, A168 (2022)

Anomalous HCN emission from warm giant molecular clouds

Authors: J. R. Goicoechea, F. Lique, M. G. Santa-Maria

Abstract: HCN is considered a good tracer of the dense molecular gas that serves as fuel for star formation. However, recent large-scale surveys of giant molecular clouds (GMCs) have detected extended HCN line emission. Such observations often resolve the HCN J=1-0 hyperfine structure (HFS). A precise determination of the physical conditions of the gas requires treating the HFS line overlap effects. Here, w… ▽ More HCN is considered a good tracer of the dense molecular gas that serves as fuel for star formation. However, recent large-scale surveys of giant molecular clouds (GMCs) have detected extended HCN line emission. Such observations often resolve the HCN J=1-0 hyperfine structure (HFS). A precise determination of the physical conditions of the gas requires treating the HFS line overlap effects. Here, we study the HCN HFS excitation and line emission using nonlocal radiative transfer models that include line overlaps and new HFS-resolved collisional rate coefficients for inelastic collisions of HCN with both para-H2 and ortho-H2 (computed via the scaled-IOS approximation up to Tk=500 K). In addition, we account for the role of electron collisions in the HFS level excitation. We find that line overlap and opacity effects frequently produce anomalous HCN J=1-0 HFS line intensity ratios (inconsistent with the common assumption of the same Tex for all HFS lines) as well as anomalous HFS line width ratios. Line overlap and electron collisions also enhance the excitation of the higher J rotational lines. Electron excitation becomes important for molecular gas with H2 densities below a few 10^5 cm-3 and electron abundances above ~10^-5. In particular, electron excitation can produce low-surface-brightness HCN emission from very extended but low-density gas in GMCs. The existence of such a widespread HCN emission component may affect the interpretation of the extragalactic relationship HCN luminosity versus star-formation rate. Alternatively, extended HCN emission may arise from dense star-forming cores and become resonantly scattered by large envelopes of lower density gas. There are two scenarios - namely, electron-assisted (weakly) collisionally excited versus scattering - that lead to different HCN J=1-0 HFS intensity ratios, which can be tested on the basis of observations. △ Less

Submitted 18 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Accepted for publication in Astronomy & Astrophysics. Abridged Abstract. This is v2 after minor editorial changes

Report number: AA/2021/42210

Journal ref: A&A 658, A28 (2022)

CF$^+$ excitation in the interstellar medium

Authors: Benjamin Desrousseaux, François Lique, Javier R. Goicoechea, Ernesto Quintas-Sánchez, Richard Dawes

Abstract: The detection of CF$^+$ in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF$^+$ abundance requires one to model its non-LTE excitation vi… ▽ More The detection of CF$^+$ in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF$^+$ abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF$^+$ in collisions with para- and ortho-H2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF$^+$ emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF$^+$ as a powerful probe to study molecular clouds exposed to strong stellar radiation fields. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Journal ref: A&A 645, A8 (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.

Comments: 10 pages, 5 figures, accepted for publication in MNRAS

Hyperfine excitation of SH$^+$ by H

Authors: François Lique, Alexandre Zanchet, Niyazi Bulut, Javier R. Goicoechea, Octavio Roncero

Abstract: SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of… ▽ More SH$^+$ is a surprisingly widespread molecular ion in diffuse interstellar clouds. There, it plays an important role triggering the sulfur chemistry. In addition, SH$^+$ emission lines have been detected at the UV-illuminated edges of dense molecular clouds, \mbox{so-called} photo-dissociation regions (PDRs), and toward high-mass protostars. An accurate determination of the SH$^+$ abundance and of the physical conditions prevailing in these energetic environments relies on knowing the rate coefficients of inelastic collisions between SH$^+$ molecules and hydrogen atoms, hydrogen molecules, and electrons. In this paper, we derive SH$^+$--H fine and hyperfine-resolved rate coefficients from the recent quantum calculations for the SH$^+$--H collisions, including inelastic, exchange and reactive processes. The method used is based on the infinite order sudden approach. State-to-state rate coefficients between the first 31 fine levels and 61 hyperfine levels of SH$^+$ were obtained for temperatures ranging from 10 to 1000 K. Fine structure-resolved rate coefficients present a strong propensity rule in favour of $Δj = ΔN$ transitions. The $Δj = ΔF$ propensity rule is observed for the hyperfine transitions. {The new rate coefficients will help significantly in the interpretation of SH$^+$ spectra from PDRs and UV-irradiated shocks where the abundance of hydrogen atoms with respect to hydrogen molecules can be significant. △ Less

Submitted 26 April, 2020; originally announced April 2020.

Comments: 3 figures

Journal ref: A&A 638, A72 (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.

Comments: Accepted by Atoms for the special issue on Atomic and molecular databases; 26 pages, 6 tables

Physical and chemical modeling of the starless core L1512

Authors: Sheng-Jun Lin, Laurent Pagani, Shih-Ping Lai, Charlène Lefèvre, François Lique

Abstract: The deuterium fractionation in starless cores gives us a clue to estimate their lifetime scales, thus allowing us to distinguish between different dynamical theories of core formation. Cores also seem to be subject to a differential N2 and CO depletion which was not expected from models. We aim to make a survey of 10 cores to estimate their lifetime scales and depletion profiles in detail. After L… ▽ More The deuterium fractionation in starless cores gives us a clue to estimate their lifetime scales, thus allowing us to distinguish between different dynamical theories of core formation. Cores also seem to be subject to a differential N2 and CO depletion which was not expected from models. We aim to make a survey of 10 cores to estimate their lifetime scales and depletion profiles in detail. After L183, in Serpens, we present the second cloud of the series, L1512 in Auriga. To constrain the lifetime scale, we perform chemical modeling of the deuteration profiles across L1512 based on dust extinction measurements from near-infrared observations and non-local thermal equilibrium radiative transfer with multiple line observations of N2H+, N2D+, DCO+, C18O, and 13CO, plus H2D+ (1$_{10}$--1$_{11}$). We find a peak density of 1.1$\times$10$^5$ cm$^{-3}$ and a central temperature of 7.5$\pm$1 K, which are respectively higher and lower compared with previous dust emission studies. The depletion factors of N2H+ and N2D+ are 27$^{+17}_{-13}$ and 4$^{+2}_{-1}$ in L1512, intermediate between the two other more advanced and denser starless core cases, L183 and L1544. These factors also indicate a similar freeze-out of N2 in L1512, compared to the two others despite a peak density one to two orders of magnitude lower. Retrieving CO and N2 abundance profiles with the chemical model, we find that CO has a depletion factor of $\sim$430-870 and the N2 profile is similar to that of CO unlike towards L183. Therefore, L1512 has probably been living long enough so that N2 chemistry has reached steady state. N2H+ modeling remains compulsory to assess the precise physical conditions in the center of cold starless cores, rather than dust emission. L1512 is presumably older than 1.4 Myr. Therefore, the dominating core formation mechanism should be ambipolar diffusion for this source. △ Less

Submitted 18 February, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: 22 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 635, A188 (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.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS (January 20 2020)

ALMA and ROSINA detections of phosphorus-bearing molecules: the interstellar thread between star-forming regions and comets

Authors: V. M. Rivilla, M. N. Drozdovskaya, K. Altwegg, P. Caselli, M. T. Beltrán, F. Fontani, F. F. S. van der Tak, R. Cesaroni, A. Vasyunin, M. Rubin, F. Lique, S. Marinakis, L. Testi, the ROSINA team

Abstract: To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated w… ▽ More To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars. △ Less

Submitted 26 November, 2019; originally announced November 2019.

Comments: Accepted in MNRAS (20 pages, 12 figures)

Effect of the isotropic collisions with neutral hydrogen on the polarization of the CN solar molecule

Authors: S. Qutub, M. Derouich, Y. N. Kalugina, H. Assiri, F. Lique

Abstract: Our work is concerned with the case of the solar molecule CN which presents conspicuous profiles of scattering polarization. We start by calculating accurate PES for the singlet and triplet electronic ground states in order to characterize the collisions between the CN molecule in its $X \; ^2Σ$ state and the hydrogen in its ground state $^2S$. The PES are included in the Schröodinger equation to… ▽ More Our work is concerned with the case of the solar molecule CN which presents conspicuous profiles of scattering polarization. We start by calculating accurate PES for the singlet and triplet electronic ground states in order to characterize the collisions between the CN molecule in its $X \; ^2Σ$ state and the hydrogen in its ground state $^2S$. The PES are included in the Schröodinger equation to obtain the scattering matrix and the probabilities of collisions. Depolarizing collisional rate coefficients are computed in the framework of the infinite order sudden approximation for temperatures ranging from $T= 2000$ K to $T= 15000$ K. Interpretation of the results and comparison between singlet and triplet collisional rate coefficients are detailed. We show that, for typical photospheric hydrogen density ($n_{H} = 10^{15}-10^{16}$ cm$^{-3}$), the $X \; ^2Σ$ state of CN is partially or completely depolarized by isotropic collisions. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS (Accepted October 22, 2019). 23 pages, 3 tables, 13 figures

C-type shock modelling -- the effect of new H$_2$-H collisional rate coefficients

Authors: A. V. Nesterenok, D. Bossion, Y. Scribano, F. Lique

Abstract: We consider collisional excitation of H$_2$ molecules in C-type shocks propagating in dense molecular clouds. New data on collisional rate coefficients for (de-)excitation of H$_2$ molecule in collisions with H atoms and new H$_2$ dissociation rates are used. The new H$_2$-H collisional data are state of the art and are based on the most accurate H$_3$ potential energy surface. We re-examine the e… ▽ More We consider collisional excitation of H$_2$ molecules in C-type shocks propagating in dense molecular clouds. New data on collisional rate coefficients for (de-)excitation of H$_2$ molecule in collisions with H atoms and new H$_2$ dissociation rates are used. The new H$_2$-H collisional data are state of the art and are based on the most accurate H$_3$ potential energy surface. We re-examine the excitation of rotational levels of H$_2$ molecule, the para-to-ortho-H$_2$ conversion, and H$_2$ dissociation by H$_2$-H collisions. At cosmic ray ionization rates $ζ\geq 10^{-16}$ s$^{-1}$ and at moderate shock speeds, the H/H$_2$ ratio at the shock front is mainly determined by the cosmic ray ionization rate. The H$_2$-H collisions play the main role in the para-to-ortho-H$_2$ conversion and, at $ζ\geq 10^{-15}$ s$^{-1}$, in the excitation of vibrationally excited states of H$_2$ molecule in the shock. The H$_2$ ortho-to-para ratio (OPR) of the shocked gas and column densities of rotational levels of vibrationally excited states of H$_2$ are found to depend strongly on the cosmic ray ionization rate. We discuss the applicability of the presented results to interpretation of observations of H$_2$ emission in supernova remnants. △ Less

Submitted 17 September, 2019; originally announced September 2019.

Comments: 10 pages, 4 figures, accepted for publication in MNRAS

Temperature and density dependent cooling function for H$_2$ with updated H$_2$/H collisional rates

Authors: Carla Maria Coppola, François Lique, Francesca Mazzia, Fabrizio Esposito, Mher V. Kazandjian

Abstract: The energy transfer among the components in a gas determines its fate. Especially at low temperatures, inelastic collisions drive the cooling and the heating mechanisms. In the early Universe as well as in zero- or low- metallicity environments the major contribution comes from the collisions among atomic and molecular hydrogen, also in its deuterated version. The present work shows some updated c… ▽ More The energy transfer among the components in a gas determines its fate. Especially at low temperatures, inelastic collisions drive the cooling and the heating mechanisms. In the early Universe as well as in zero- or low- metallicity environments the major contribution comes from the collisions among atomic and molecular hydrogen, also in its deuterated version. The present work shows some updated calculations of the H$_2$ cooling function based on novel collisional data which explicitely take into account the reactive pathway at low temperatures. Deviations from previous calculations are discussed and a multivariate data analysis is performed to provide a fit depending on both the gas temperature and the density of the gas. △ Less

Submitted 2 July, 2019; originally announced July 2019.

Comments: 5 pages, 2 figures, 1 table

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 486, Issue 2, June 2019, Pages 1590-1593

Collisional excitation of NH(3Σ-) by Ar: A new ab initio 3D potential energy surface and scattering calculations

Authors: D. Prudenzano, F. Lique, R. Ramachandran, L. Bizzocchi, P. Caselli

Abstract: Collisional excitation of light hydrides is important to fully understand the complex chemical and physical processes of atmospheric and astrophysical environments. Here, we focus on the NH(X3Σ-)-Ar van der Waals system. First, we have calculated a new three-dimensional Potential Energy Surface (PES), which explicitly includes the NH bond vibration. We have carried out the ab initio calculations o… ▽ More Collisional excitation of light hydrides is important to fully understand the complex chemical and physical processes of atmospheric and astrophysical environments. Here, we focus on the NH(X3Σ-)-Ar van der Waals system. First, we have calculated a new three-dimensional Potential Energy Surface (PES), which explicitly includes the NH bond vibration. We have carried out the ab initio calculations of the PES employing the open-shell single- and double-excitation couple cluster method with noniterative perturbational treatment of the triple excitations. To achieve a better accuracy, we have first obtained the energies using the augmented correlation-consistent aug-cc-pVXZ (X = T, Q, 5) basis sets and then we have extrapolated the final values to the complete basis set limit. We have also studied the collisional excitation of NH(X3Σ-)-Ar at the close-coupling level, employing our new PES. We calculated collisional excitation cross sections of the fine-structure levels of NH by Ar for energies up to 3000 cm-1 . After thermal average of the cross sections, we have then obtained the rate coefficients for temperatures up to 350 K. The propensity rules between the fine-structure levels are in good agreement with those of similar collisional systems, even though they are not as strong and pronounced as for lighter systems, such as NH-He. The final theoretical values are also compared with the few available experimental data. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Journal ref: The Journal of Chemical Physics 150.21 (2019): 214302

Formation of interstellar SH$^+$ from vibrationally excited H$_2$: Quantum study of S$^+$ + H$_2$ $\rightleftarrows$ SH$^+$ + H reactions and inelastic collisions

Authors: Alexandre Zanchet, Francois Lique, Octavio Roncero, Javier R. Goicoechea, Niyazi Bulut

Abstract: The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3Σ^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H… ▽ More The rate constants for the formation, destruction, and collisional excitation of SH$^+$ are calculated from quantum mechanical approaches using two new SH$_2^+$ potential energy surfaces (PESs) of $^4A''$ and $^2A''$ electronic symmetry. The PESs were developed to describe all adiabatic states correlating to the SH$^+$ ($^3Σ^-$) + H($^2S$) channel. The formation of SH$^+$ through the S$^+$ + H$_2$ reaction is endothermic by $\approx$ 9860 K, and requires at least two vibrational quanta on the H$_2$ molecule to yield significant reactivity. Quasi-classical calculations of the total formation rate constant for H$_2$($v=2$) are in very good agreement with the quantum results above 100K. Further quasi-classical calculations are then performed for $v=3$, 4, and 5 to cover all vibrationally excited H$_2$ levels significantly populated in dense photodissociation regions (PDR). The new calculated formation and destruction rate constants are two to six times larger than the previous ones and have been introduced in the Meudon PDR code to simulate the physical and illuminating conditions in the Orion bar prototypical PDR. New astrochemical models based on the new molecular data produce four times larger SH$^+$ column densities, in agreement with those inferred from recent ALMA observations of the Orion bar. △ Less

Submitted 7 May, 2019; originally announced May 2019.

Comments: 8 pages, 7 figures

Report number: Accepted in Astronomy and Astrophysics 2019

Journal ref: A&A 626, A103 (2019)

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.

Comments: 16 pages, 4 figures

Journal ref: The Journal of Physical Chemistry Letters 2018 9 3199-3204

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.

Comments: 16 pages, 12 figures

Journal ref: A&A 615, A52 (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.

Journal ref: ApJL 853 L5 (2018)

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.

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.

Comments: Accepted for publication in A&A

Journal ref: A&A 597, A40 (2017)

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.

Comments: 8 pages, 2 tables, 4 figures, accepted for publication in MNRAS

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.

Stratified NH and ND emission in the prestellar core 16293E in L1689N

Authors: A. Bacmann, F. Daniel, P. Caselli, C. Ceccarelli, D. Lis, C. Vastel, F. Dumouchel, F. Lique, E. Caux

Abstract: High degrees of deuterium fractionation are commonly found in cold prestellar cores and in the envelopes around young protostars. As it brings strong constraints to chemical models, deuterium chemistry is often used to infer core history or molecule formation pathways. Whereas a large number of observations is available regarding interstellar deuterated stable molecules, relatively little is known… ▽ More High degrees of deuterium fractionation are commonly found in cold prestellar cores and in the envelopes around young protostars. As it brings strong constraints to chemical models, deuterium chemistry is often used to infer core history or molecule formation pathways. Whereas a large number of observations is available regarding interstellar deuterated stable molecules, relatively little is known about the deuteration of hydride radicals, as their fundamental rotational transitions are at high frequencies where the atmosphere is mostly opaque. Nitrogen hydride radicals are important species in nitrogen chemistry, as they are thought to be related to ammonia formation. Observations have shown that ammonia is strongly deuterated, with [NH2D]/[NH3] ~ 10%. Models predict similarly high [ND]/[NH] ratios, but so far only one observational determination of this ratio is available, towards the envelope of the protostar IRAS16293-2422. In order to test model predictions, we aim here at determining [ND]/[NH] in a dense, starless core. We observed NH and ND in 16293E with the HIFI spectrometer on board the Herschel Space Observatory as part of the CHESS guaranteed time key programme, and derived the abundances of these two species using a non-LTE non-local radiative transfer model. Both NH and ND are detected in the source, with ND in emission and NH in absorption against the continuum arising from the cold dust emission. Our model shows however that the ND emission and the NH absorption originate from different layers in the cloud, as further evidenced by their different velocities. In the central region of the core, we can set a lower limit to the [ND]/[NH] ratio of ~2%. This estimate is consistent with recent pure gas-phase models of nitrogen chemistry △ Less

Submitted 8 December, 2015; originally announced December 2015.

Comments: accepted by A&A

Journal ref: A&A 587, A26 (2016)

Si-bearing molecules towards IRC+10216: ALMA unveils the molecular envelope of CWLeo

Authors: L. Velilla Prieto, J. Cernicharo, G. Quintana-Lacaci, M. Agúndez, A. Castro-Carrizo, J. P. Fonfría, N. Marcelino, J. Zúñiga, A. Requena, A. Bastida, F. Lique, M. Guélin

Abstract: We report the detection during the ALMA Cycle 0 of SiS rotational lines in high-vibrational states as well as SiO and SiC$_2$ lines in their ground vibrational state, towards IRC+10216. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC$_2$, and also proves the existence of an increase in the SiC$_2$ emission at the outer shells… ▽ More We report the detection during the ALMA Cycle 0 of SiS rotational lines in high-vibrational states as well as SiO and SiC$_2$ lines in their ground vibrational state, towards IRC+10216. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC$_2$, and also proves the existence of an increase in the SiC$_2$ emission at the outer shells of the circumstellar envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram technique and we used a large velocity gradient model to compare with the observations. We found moderate discrepancies between the observations and the models that could be explained if SiS lines detected are optically thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show a decreasing linewidth with increasing energy levels. This may evidence that these lines could be excited only in the inner shells, i.e. the densest and hottest, of the circumstellar envelope of IRC+10216. △ Less

Submitted 11 June, 2015; v1 submitted 22 April, 2015; originally announced April 2015.

Comments: 20 pages, 3 figures. Published in ApJL

Journal ref: ApJ 805 L13, 2015

Depletion of chlorine into HCl ice in a protostellar core

Authors: M. Kama, E. Caux, A. Lopez-Sepulcre, V. Wakelam, C. Dominik, C. Ceccarelli, M. Lanza, F. Lique, B. B. Ochsendorf, D. C. Lis, R. N. Caballero, A. G. G. M. Tielens

Abstract: The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-… ▽ More The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value. △ Less

Submitted 24 November, 2014; originally announced November 2014.

Comments: 13 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 574, A107 (2015)

Collisional excitation of water by hydrogen atoms

Authors: F. Daniel, A. Faure, P. J. Dagdigian, M. -L. Dubernet, F. lique, G. Pineau des forêts

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

Submitted 29 October, 2014; originally announced October 2014.

Hyperfine excitation of N$_2$H$^+$ by H$_2$: Toward a revision of N$_2$H$^+$ abundance in cold molecular clouds

Authors: François Lique, Fabien Daniel, Laurent Pagani, Nicole Feautrier

Abstract: The modelling of emission spectra of molecules seen in interstellar clouds requires the knowledge of collisional rate coefficients. Among the commonly observed species, N$_2$H$^+$ is of particular interest since it was shown to be a good probe of the physical conditions of cold molecular clouds. Thus, we have calculated hyperfine-structure resolved excitation rate coefficients of N$_2$H$^+$(X… ▽ More The modelling of emission spectra of molecules seen in interstellar clouds requires the knowledge of collisional rate coefficients. Among the commonly observed species, N$_2$H$^+$ is of particular interest since it was shown to be a good probe of the physical conditions of cold molecular clouds. Thus, we have calculated hyperfine-structure resolved excitation rate coefficients of N$_2$H$^+$(X$^1Σ^+$) by H$_2(j=0)$, the most abundant collisional partner in the cold interstellar medium. The calculations are based on a new potential energy surface, obtained from highly correlated {\it ab initio} calculations. State-to-state rate coefficients between the first hyperfine levels were calculated, for temperatures ranging from 5 K to 70 K. By comparison with previously published N$_2$H$^+$-He rate coefficients, we found significant differences which cannot be reproduced by a simple scaling relationship. As a first application, we also performed radiative transfer calculations, for physical conditions typical of cold molecular clouds. We found that the simulated line intensities significantly increase when using the new H$_2$ rate coefficients, by comparison with the predictions based on the He rate coefficients. In particular, we revisited the modelling of the N$_2$H$^+$ emission in the LDN 183 core, using the new collisional data, and found that all three of the density, gas kinetic temperature and N$_2$H$^+$ abundance had to be revised. △ Less

Submitted 20 October, 2014; originally announced October 2014.

OH+ in astrophysical media: state-to-state formation rates, Einstein coefficients and inelastic collision rates with He

Authors: S. Gomez-Carrasco, B. Godard, F. Lique, N. Bulut, J. Klos, O. Roncero, A. Aguado, F. J. Aoiz, J. F. Castillo, J. R. Goicoechea, M. Etxaluze, J. Cernicharo

Abstract: The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods. First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3Σ^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-i… ▽ More The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods. First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3Σ^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-independent results to asses the accuracy of reaction probabilities at collision energies of about 1 meV. The good agreement between the simulations and the existing experimental cross sections in the $0.01-$1 eV energy range shows the quality of the results. The calculated state-to-state rate constants have been fitted to an analytical form. Second, the Einstein coefficients of OH$^+$ have been obtained for all astronomically significant ro-vibrational bands involving the $X^3Σ^-$ and/or $A^3Π$ electronic states. For this purpose the potential energy curves and electric dipole transition moments for seven electronic states of OH$^+$ are calculated with {\it ab initio} methods at the highest level and including spin-orbit terms, and the rovibrational levels have been calculated including the empirical spin-rotation and spin-spin terms. Third, the state-to-state rate constants for inelastic collisions between He and OH$^+(X ^3Σ^-)$ have been calculated using a time-independent close coupling method on a new potential energy surface. All these rates have been implemented in detailed chemical and radiative transfer models. Applications of these models to various astronomical sources show that inelastic collisions dominate the excitation of the rotational levels of OH$^+$. In the models considered the excitation resulting from the chemical formation of OH$^+$ increases the line fluxes by about 10 % or less depending on the density of the gas. △ Less

Submitted 16 May, 2014; originally announced May 2014.

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

Authors: M. Etxaluze, J. Cernicharo, J. R. Goicoechea, P. A. M. van Hoof, B. M. Swinyard, M. J. Barlow, G. C. van de Steene, M. A. T. Groenewegen, F. Kerschbaum, T. L. Lim, F. Lique, M. Matsuura, C. Pearson, E. T. Polehampton, P. Royer, T. Ueta

Abstract: The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 S… ▽ More The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium. △ Less

Submitted 8 April, 2014; originally announced April 2014.

Comments: Accepted for publication in Astronomy and Astrophysics

Molecular excitation in the Interstellar Medium: recent advances in collisional, radiative and chemical processes

Authors: Evelyne Roueff, Francois Lique

Abstract: We review the different excitation processes in the interstellar medium We review the different excitation processes in the interstellar medium △ Less

Submitted 30 October, 2013; originally announced October 2013.

Comments: Accepted in Chem. Rev

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

Authors: F. Daniel, M. Gerin, E. Roueff, J. Cernicharo, N. Marcelino, F. Lique, D. C. Lis, D. Teyssier, N. Biver, D. Bockelee--Morvan

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

Submitted 23 September, 2013; originally announced September 2013.

Comments: accepted in A&A

Astronomical identification of CN-, the smallest observed molecular anion

Authors: M. Agundez, J. Cernicharo, M. Guelin, C. Kahane, E. Roueff, J. Klos, F. J. Aoiz, F. Lique, N. Marcelino, J. R. Goicoechea, M. Gonzalez Garcia, C. A. Gottlieb, M. C. McCarthy, P. Thaddeus

Abstract: We present the first astronomical detection of a diatomic negative ion, the cyanide anion CN-, as well as quantum mechanical calculations of the excitation of this anion through collisions with para-H2. CN- is identified through the observation of the J = 2-1 and J = 3-2 rotational transitions in the C-star envelope IRC +10216 with the IRAM 30-m telescope. The U-shaped line profiles indicate that… ▽ More We present the first astronomical detection of a diatomic negative ion, the cyanide anion CN-, as well as quantum mechanical calculations of the excitation of this anion through collisions with para-H2. CN- is identified through the observation of the J = 2-1 and J = 3-2 rotational transitions in the C-star envelope IRC +10216 with the IRAM 30-m telescope. The U-shaped line profiles indicate that CN-, like the large anion C6H-, is formed in the outer regions of the envelope. Chemical and excitation model calculations suggest that this species forms from the reaction of large carbon anions with N atoms, rather than from the radiative attachment of an electron to CN, as is the case for large molecular anions. The unexpectedly large abundance derived for CN-, 0.25 % relative to CN, makes likely its detection in other astronomical sources. A parallel search for the small anion C2H- remains so far unconclusive, despite the previous tentative identification of the J = 1-0 rotational transition. The abundance of C2H- in IRC +10216 is found to be vanishingly small, < 0.0014 % relative to C2H. △ Less

Submitted 5 July, 2010; originally announced July 2010.

Comments: 5 pages, 4 figures; accepted for publication in A&A Letters

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

Authors: E. Sarrasin, D. Ben Abdallah, M. Wernli, A. Faure, J. Cernicharo, F. Lique

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

Submitted 26 February, 2010; originally announced March 2010.

Comments: 8 figues