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 $\Delta j = \Delta N$ transitions. The $\Delta j = \Delta 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.

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.

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.

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.

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.

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\Sigma$ state and the hydrogen in its ground state $^2S$. The PES are included in the Schr\"oodinger 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\Sigma$ state of CN is partially or completely depolarized by isotropic collisions.

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 $\zeta \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 $\zeta \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.

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.

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{\Sigma}-)-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{\Sigma}-)-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.

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\Sigma^-$) + 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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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\Sigma^+$) 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.

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\Sigma^-, 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\Sigma^-$ and/or $A^3\Pi$ 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\Sigma^-)$ 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.

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.

We review the different excitation processes in the interstellar medium

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.

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.

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 $\Delta j$ transitions is observed in the case of HCN--He collisions whereas a propensity for odd $\Delta 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.